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1.0  Introduction 

The  ICPIC  workstation  is  an  extensible  simulation  testbed,  with  a  graphical,  menu-driven 
interface,  that  facilitates  plug  and  play  evaluation  and  design  of  optimal  image  compression  and 
communication  system  configurations.  The  ICPIC  testbed  can  be  used  to  simulate  the  entire 
communication  process  from  original  image  retrieval,  through  image  compression,  error  coding, 
transmission,  propagation  over  a  channel,  error  detection  and  correction  techniques  (EDAC), 
image  decompression,  and  final  image  display  and  evaluation. 

In  addition,  image  compression  algorithms  can  be  evaluated  with  respect  to  their 
performance  when  used  for  image  transmission  over  a  variety  of  communication  channels.  The 
results  of  these  evaluations  can  then  be  used  to  suggest  algorithm  modifications  which  will 
improve  performance  over  a  communication  channel. 

The  ICPIC  testbed  is  modular  in  design  and  is  built  on  the  Signal  Processing  Worksystem 
(SPW)  software  developed  by  the  Alta  Group  of  Cadence  Design  Systems,  Inc.  and  the  Spread 
Spectrum  Vulnerability  Metric  (SSVM)  software  developed  by  TASC.  Use  of  the  ICPIC  testbed 
requires  familiarity  only  with  the  material  covered  in  the  SPW  tutorial  naanual.  ICPIC  expands  the 
SPW  analysis  capabilities  by  providing  the  user  with  one  standardized  graphical  interface  with 
which  to  gain  parallel  access  to  all  of  the  capabilities  of  both  the  SPW  and  SSVM  module  libraries. 

ICPIC  facilitates  the  construction  of  a  complete  simulation  from  analysis  “primitives”  by 
drag  and  drop  operations.  These  primitives  range  in  capability  from  low-level  modules  such  as 
signal  sources  and  simple  scaler  multipliers,  to  high-level  modules  that  perform  such  capabilities  as 
PSK  (Phase  Shift  Keying)  and  EDAC  (Error  Detection  and  Correction).  By  utilizing  the  graphical 
interface  along  with  the  SPW  and  SSVM  module  libraries,  ICPIC  gives  the  user  the  capability  to 
quickly  construct  an  end-to-end  simulation,  allowing  for  rapid  system  design  assessment  and 
optimization. 

The  remainder  of  this  section  describes  the  image  compression,  EDAC,  and  channel 
propagation  models  available  in  ICPIC  as  well  as  the  relationship  between  these  components. 
Section  2  provides  a  detailed  discussion  of  a  complete  end-to-end  BPSK  simulation,  including 
both  the  image  analysis  and  data  preparation  software  tools  which  were  developed  for  this  effort 
and  are  available  as  part  of  the  ICPIC  testbed.  Section  3  presents  sample  data  outputs  from  the 
BPSK  system.  Section  4  presents  a  summary  and  recommendations  for  future  analysis. 
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1.1  Image  Compression 


Image  compression  algorithms  reduce  storage  and  bandwidth  requirements  by  reducing 
information  redundancy  in  the  transmitted  image.  The  smaller  the  image,  the  less  rime  the 
propagation  channel  has  to  degrade  it  by  the  addition  of  noise.  Individual  image  compression 
algorithms  are  grouped  under  the  higher  level  classifications  of  either  Lossless  or  Lossy 
compression  algorithms.  Lossless  means  that  the  image  data  can  be  compressed  and 
decompressed  infinitely  many  times  without  suffering  any  loss  in  the  information  content  With 
lossy  compression  algorithms,  there  is  some  information  loss  every  time  the  image  is  compressed 
and  decompressed.  The  amount  of  loss  depends  on  the  compression  factor,  which  is  an  arbitrary 
user-defined  integer.  Depending  on  the  application,  it  may  be  acceptable  to  lose  some  information 
if  it  significantly  decreases  the  transmission  time  while  maintaining  sufficient  information  content 

Depending  on  the  specific  image  compression  technique,  certain  portions  of  an  image  may 
be  more  sensitive  to  errors  than  other  portions.  The  ED  AC  techiuques  can  therefore  be  optimized 
for  the  more  sensitive  portions  of  the  transmitted  image.  In  general,  errors  in  the  image  header 
need  to  be  eliminated  or  at  least  minimized  as  the  information  contained  in  the  header  portion  of  the 
image  contains  the  prescription  (such  as  the  horizontal  and  vertical  resolution)  for  how  to 
reconstruct  the  transmitted  image  fiom  the  received  data 

To  emphasize  the  need  for  image  compression  consider  the  following  example.  A  digitized 
grey  scale  image  consists  of  an  array  of  picture  elements  (pixels)  represented  by  numbers  that 
correspond  to  the  brightness  (grey  scale)  of  each  pixel,  with  0  =  black  and  255  =  white.  A  grey¬ 
scale  image  has  1024  x  1024  such  numbers  taking  integer  values  between  0  and  255.  Thus  the 

image  is  given  by  a  matrix  (Pj)je{0,...,1023}2  with  Pj  e  {0,...,255}.  To  encode  this  sample 
image  would  require  a  data  file  approximately  one  million  bytes  in  size.  Compression  of  this 
sample  image  with  the  appropriate  image  compression  algorithm  will  facilitate  the  efficient  and 
effective  transmission  of  the  image  across  a  given  communication  system  configuration. 

For  the  ICPIC  effort  the  following  three  image  compression  techniques  were  utilized:  the 
Graphics  Interchange  Format  (GIF),  the  Joint  Photographic  Experts  Group  (JPEG),  and  Wavelet. 
GIF  is  a  popular  commercial  standard  8  bil/256  color  protocol  for  transmission  and  interchange  of 
raster  graphic  data  in  a  way  that  is  independent  of  the  hardware  used  in  it’s  creation  or  display. 
GIF  is  comprised  of  a  lossless  compression  strategy  based  on  LZW  compression  (Lempel-Ziv- 
Welch  algorithm).  The  LZW  technique  was  originally  developed  for  compression  of  textual 
material,  and  compresses  files  by  substituting  commonly  occurring  character  sequences  with  a 
reference  to  the  first  occurrence  of  the  sequence.  For  fullscreen,  8  bit  images  of  moderate 
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complexity,  4:1  compression  is  the  average.  Since  this  is  a  lossless  compression  scheme,  the 
compression  ratio  can  not  be  increased  by  trading  off  image  quality. 

GIF  compression  offers  the  following  advantages  and  strengths: 

•  it  is  a  lossless  compression  scheme  for  8  bit  images  (since  no  image  degradation 
occurs,  repeated  compression  and  decompression  cycles  are  possible,  also 
suitable  for  images  where  information  loss  can  not  be  tolerated). 

•  it  is  ideally  suited  to  styli2ed  images  such  as  line  drawings,  or  those  which  contain 
only  a  limited  number  of  colors  (for  these  images  it  can  produce  good 
compression  ratios). 

•  it  is  widely  used  and  supported,  with  no  runtime  license  required. 

The  weaknesses  and  disadvantages  of  GIF  compression  are  as  follows: 

•  it  is  not  suitable  for  24  bit  images.  When  compressing  such  images,  much  of  the 
color  information  is  lost  during  the  quantization  process  which  reduces  this 

to  8  bits.  Good  algorithms  can  however  optimize  this  process  so  that  the 
resultant  image  will  still  have  good  quality  from  a  human  point  of  view. 

•  the  compression  ratios  are  low.  These  can  not  be  traded  off  against  compression 
times  or  degree  of  loss  of  quality. 

•  it  is  not  intended  for  moving  images/video. 

•  it  is  not  resolution  independent. 

JPEG  is  a  lossy  image  compression  scheme  aimed  at  still  images,  in  color  or  grey  scale, 
and  having  some  degree  of  complexity.  The  greater  the  compression  ratio,  the  greater  the  degree 
of  information  loss.  This  can  be  user  determined  to  optimize  the  trade-off  between  resultant  image 
size  and  image  quality.  The  JPEG  compression  algorithm  exploits  some  of  the  ways  in  which  the 
human  eye  perceives  and  analyzes  images,  so  that  compressed  images  still  appear  to  be  of  high 
quality  when  looked  at  by  human  eyes.  The  total  amount  of  time  it  takes  to  compress  and 
decompress  an  image  are  symmetric,  meaning  that  they  take  roughly  the  same  amount  of  time. 

The  JPEG  algorithm  is  based  on  the  forward  DCT  (Discrete  Cosine  Transform),  as  applied 
to  a  block  breakdown  of  the  original  image  into  8  by  8  blocks,  quantized  down  to  a  finite  set  of 
possible  values,  and  then  further  transformed  (run  length  encoding)  and  finally  entropy  encoded 
using  Huffman  or  arithmetic  coding. 
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The  following  list  compares  the  compression  ratios  with  the  observed  quality.  A  20:1 
compression  ratio  means  that  an  image  originally  900K  will  be  compressed  to  45K,  which  is  one 
twentieth  of  its  original  size. 

•  10:1  to  20:1  -  High  quality,  with  little  or  no  observable  loss  in  image  quality  to  the 
human  viewer. 

•  30: 1  to  50: 1  -  Moderate  quality. 

•  60:1  to  100:1  -  Poor  quality,  suitable  for  thumbnails  and  previews.  Marked 
blockiness  and  “Gibb’s  effect”  occurs.  The  Gibb’s  effect  is  the  name  given 
to  the  phenomenon  where  disturbances  and/or  ripples  can  be  seen  at  the 
margins  of  objects  with  sharp  borders. 

In  comparison  with  GIF  compressed  images,  high  quality  JPEG  image  compression 
produces  an  image  4  to  5  times  smaller. 

JPEG  image  compression  offers  the  following  advantages  and  strengths: 

•  it  provides  support  for  full  24  bit  color  images.  In  contrast,  GIF  compression 
only  supports  8  bit  images. 

•  the  compressed  image  size  and  image  quality  trade-off  can  be  determined  by  the 
user. 

•  it  is  ideally  suited  to  images  of  real  world  scenes,  or  complex  computer  generated 
images. 

•  it  is  platform  independent  for  displaying  24  bit  images.  It  is  currently  the  most 
widely  adhered  to  image  compression  standard,  with  the  algorithm,  source  code 
implementations,  and  public  domain  viewers  readily  available. 

The  weaknesses  and  disadvantages  of  JPEG  are  as  follows: 

•  JPEG  compression  is  a  trade-off  between  the  degree  of  compression,  the  resultant 
image  quality  and  the  amount  of  time  required  to  compress  and  decompress  an 
image.  Also,  blockiness  results  at  the  high  image  compression  ratios. 

•  it  produces  poor  image  quality  when  compressing  text  or  images  containing  sharp 
edges  or  lines  (known  as  the  Gibb’s  effect). 

•  it  is  not  suitable  for  2  bit  black  and  white  images. 

•  the  degree  of  compression  is  greater  for  full  color  images  than  it  is  for  grey  scale 
images. 

•  it  is  not  a  suitable  sffategy  for  images  that  are  still  being  edited,  because  every 
compression  and  decompression  cycle  causes  a  loss  in  information. 
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•  it  is  not  intended  for  moving  images  and/or  video. 

•  it  is  not  resolution  independent.  Does  not  provide  for  scalability,  where  the  image 
is  displayed  optimally  depending  on  the  resolution  of  the  viewing  device. 

Wavelet  transforms  are  high  compression  algorithms  capable  of  analyzing  data  at  multiple 
resolutions  (also  known  as  “scale”).  In  other  words,  they  cut  up  data  into  different  frequency 
components,  and  then  analyze  each  component  with  a  resolution  matched  to  its  scale.  In  addition, 
transient  events  in  the  data  are  preserved  by  the  analysis.  When  the  wavelet  transform  is  applied  to 
a  signal  in  the  time  domain,  the  result  is  a  two  dimensional,  time-scale  analysis  of  the  signal.  The 
transform  of  the  data  exhibits  discrete  steps  in  time  on  one  axis,  and  discrete  steps  of  resolution  on 
the  other.  The  wavelet  transform  has  proven  to  be  a  useful  tool  for  the  compression  and  analysis 
of  both  signals  and  images.  The  discrete  wavelet  transform  (DWT)  is  the  wavelet  transform  as 
applied  to  a  regularly  sampled  data  sequence  and  the  fast  wavelet  transform  (FWT)  is  an  efficient 
implementation  of  the  DWT. 

DWT  has  the  distinct  benefit  of  being  able  to  perform  the  simultaneous  localization  of 
frequency  and  time.  This  gives  the  DWT  the  instantaneous  advantage  of  information  condensing. 
The  wavelet  transform  reduces  information  in  a  signal  roughly  to  averages  and  differences  of 
neighboring  pixels  at  discrete  scale  (frequency)  levels.  In  an  image,  the  information  in  areas  of 
near  constant  color  are  collapsed  into  the  average  component  and  the  difference  components.  Each 
average  and  difference  component  contains  information  about  an  area  4  times  its  size.  If  the 
differences  are  close  to  zero  then  the  components  can  be  discarded  and  all  the  information  about  the 
area  will  be  contained  in  the  average  components. 

The  wavelet  analysis  procedure  adopts  a  wavelet  prototype  function,  called  an  analyzing 
wavelet  or  mother  wavelet.  In  the  FWT  algorithm,  the  sampled  data  set  is  passed  through  scaling 
and  wavelet  filters.  These  filters  are  collectively  known  as  a  quadrature  mirror  filter  (QMF)  pair 
and  are  actually  low-pass  and  high-pass  filters  which  contain  complementary  bandwidths.  The 
outputs  of  both  filters  are  decimated  (desampled)  by  a  factor  of  two.  The  high-pass  filtered  4ata  set 
is  the  wavelet  transform  detail  coefficients  at  the  same  level  of  scale  as  the  transform.  The  low- 
pass  filtered  data  set  is  the  approximation  coefficients  at  the  same  level  of  scale.  Temporal  analysis 
is  performed  by  using  a  high-frequency  version  of  the  prototype  wavelet,  while  frequency  analysis 
is  performed  with  a  dilated,  low-frequency  version  of  the  same  wavelet.  Due  to  the  decimation, 
both  sets  of  coefficients  have  half  as  many  elements  as  the  original  data  set.  The  approximation 
coefficients  can  now  be  used  as  the  sampled  data  input  for  another  pair  of  wavelet  filters,  identical 
to  the  first  pair,  generating  another  set  of  detail  and  approximation  coefficients  at  the  next  lower 
level  of  scale.  This  process  can  continue  until  the  limit  for  the  unit  interval  is  reached.  Because  the 
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original  signal  or  function  can  be  represented  in  terms  of  wavelet  expansion  (using  coefficients  in  a 
linear  combination  of  the  wavelet  functions),  data  operations  can  be  performed  using  just  the 
corresponding  wavelet  coefficients.  Further  compression  is  achieved  by  keeping  only  the  most 
significant  bits  of  the  coefficients  which  is  known  as  scaler  quanti2ation. 

Threshold  coding  is  the  process  whereby  only  those  coefficients  that  exceed  a  specified 
threshold  are  retained.  Threshold  coding  results  in  the  sparse  representation  of  the  data  which 
makes  wavelets  an  excellent  tool  for  data  compression.  Optimally,  the  coefficients  are  transmitted 
in  order  of  decreasing  size.  The  process  of  sending  coefficients  sequentially  across  a 
communication  link  to  allow  for  the  gradual  reconstruction  of  the  image  is  known  as  progressive 
transmission.  Flexibility  is  a  characteristic  that  makes  the  FWT  an  important  tool  for  image 
transmission.  The  FWT  lends  itself  ideally  to  applications  which  require  reordering  of  digital 
information.  Its  adaptability  is  best  seen  in  the  progressive  transmission  of  an  image  over  a  digital 
data  link  where  incoming  information  builds  on  that  already  received,  resulting  in  no  computational 
or  information  overhead.  In  other  words,  the  time  and  information  it  takes  to  decompress  an  image 
piecewise  is  the  same  as  decompressing  the  image  all  at  once. 

As  stated  previously,  the  JPEG  image  compression  method  is  based  on  a  division  of  an 
image  into  blocks  of  8  x  8  pixels,  after  which  each  block  is  transformed  with  a  discrete  cosine 
transform.  While  the  JPEG  method  performs  well  at  high  or  medium  bit  rates,  it  introduces 
perceptibly  annoying  blocking  artifacts  (due  to  block-wise  transform  coding)  at  low  bit  rates.  The 
performance  of  JPEG  deteriorates  rapidly  at  compression  ratios  above  30:1.  Due  to  good 
localization  in  both  space  and  frequency  domains  of  the  basis  function,  the  wavelet  transform  can 
obtain  greater  compression  ratios  thereby  facilitating  high  data  throughput  which  is  required  for 
real-time/video  and  image  compression.  As  a  result,  the  reconstruction  quality  of  wavelet 
compressed  images  has  moved  well  beyond  the  capabilities  of  the  JPEG  method. 

1.2  Error  Correction  and  Detection 

When  information  is  propagated  between  two  points,  errors  are  introduced  due  to  noise 
along  the  propagation  channel.  Two  things  are  therefore  necessary;  error  detection,  and  (if 
possible)  error  correction. 

The  simplest  form  of  error  correction  is  to  use  an  error-detecting  code  so  conceived  that  any 
false  or  incorrect  signal  initiates  a  repetition  of  the  transmitted  character  that  was  incorrectly 
received.  This  form  increases  not  only  the  transmission  time,  but  also  the  probability  of  signal 
detection  and  interception  when  this  may  be  an  issue  of  concern. 
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A  better  method  is  to  divide  the  data  into  packets  of  fixed  size,  where  each  packet  contains 
image  information  and  a  code  that  allows  a  user  to  detect,  with  high  probability,  if  an  error  has 
been  introduced.  Based  upon  the  known  size  of  the  data  packet  (N),  and  the  number  of  bits  of 
information  (k),  it  is  usually  possible  to  remove  the  errors  and  reconstruct  the  original  image. 

SPW  provides  a  number  of  EDAC  algorithms,  such  as  Reed-Solomon  and  Bose- 
Chaudhuri-Hocquenghem  (BCH)  which  includes  a  vector  BCH  encoder  and  decoder.  BCH 
encoding  and  decoding  processes  binary  data  in  blocks.  The  convolutional  encoder,  which  is  often 
used  with  Viterbi  decoders,  typically  works  with  arbitrarily  long,  continuous  flowing  bit  streams. 
The  Viterbi  algorithms  have  also  provided  competitive  performance  with  blocks  of  data.  SPW 
provides  both  the  hard  and  soft  decision  versions  of  the  Viterbi  decoding  algorithm,  with  the  soft 
decision  version  permitting  substantially  fewer  residual  errors.  Viterbi  decoding  is  the  most 
commonly  used  form  of  trellis  decoding,  and  requires  substantially  less  storage  of  prior  data  than 
sequential  decoding.  Many  satellite  communications  systems  use  Viterbi  decoding. 

Reed-Solomon  encoding  and  decoding  is  well  suited  for  use  with  M-ary  signaling,  and  also 
for  error  distributions  that  are  uniform  in  the  symbol  space.  JTIDS  uses  Reed-Solomon  coding. 
Included  in  the  SPW  software  are  vector  Reed-Solomon  coding  and  decoding  and  included  in  the 
ICPIC  module  libraries  are  M-ary  modulation  and  demodulation  to  complement  the  Reed-Solomon 
coding  modules. 

In  addition,  TCM  and  V.32  encoders  and  decoders  are  available  as  part  of  the  ICPIC 
module  libraries.  V.32  coding  is  an  ISO  standard  for  telephone  line  transirdssion. 

1.3  Propagation  Channel 

The  propagation  channel  describes  the  interaction  between  the  signal  and  the  environment 
between  the  transmitter  and  receiver.  Among  the  channel  models  incorporated  into  ICPIC  are 
noise,  interference,  tapped  delay,  and  path  loss. 

The  term  background  refers  to  additive  contamination  of  the  communication  signal  and 
includes  both  noise  and  interference.  Noise  is  generally  natural  in  origin  and  broadband,  and 
includes  thermal  receiver  noise  and  environmental  or  atmospheric  noise  due  to  natural 
electromagnetic  effects.  Interference  is  generally  man-made  and  often  narrowband.  Table  1-1, 
Independent  Parameters  in  Background  and  Channel  Models,  shows  the  occurrences  and 
properties  of  the  independent  parameters  used  in  background  and  channel  components. 
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Independent  Parameter 

Occurrence 

Relationships 

Conflguradon  center  frequency 

Random  Interference 

HF  Interference 

VHF  Interference 

lypically  inside  signal  band 

Configuration  bandwidth 

Random  Interference 

HF  Interference 

VHF  Interference 

Interferers  outside  detector  input  bandwidth 
have  no  effect  on  vulnerability 

Configuration  input  file 

Name  of  file  used  to  store  random  interference 
configuration 

Custom  Interference 

Format  same  as  configuration  output  file 
read  by  Custom  Interference  nradule 

Configuration  input  file 

Name  of  file  used  to  store  random  interference 
configuration 

Random  Interference 

Format  same  as  configuration  output  file 
read  by  Custom  Interference  module 

Input  SNR 

Gaussian  Noise 

Ganuna  Noise 

Assumes  only  one  noise  source  in  use 

Exponent  describing  power  loss  with 
propagation  distance 

Path  Loss 

Frequency  Resolution 

Spacing  of  grid  of  candidate  interferer  carrios 

Random  Interference 

Less  than  configuration  bandwidth 

Minimum  interferer  width 

Random  Interference 

Less  than  maximum  interferer  width 

Maximum  interferer  width 

Random  Interference 

Greater  than  minimum  interfere  width 

Mean  Power  (dB  above  signal) 

Random  Interference 

Number  of  Interferers 

Random  Interference 

Integer,  at  least  one 

Power  standard  deviation  (dB) 

Random  Int^erence 

Greater  than  one  tenth  of  the  mean  interferer 
power 

Scale  factor 

Controls  strength  of  Cauchy  noise,  which  has 
infinite  variance 

Cauchy  Noise 

Input  SNR  is  not  defined 

Spread 

Time  intoval  between  first  and  last  taps 

Thpped  Delay 

Less  than  interception  collect  time 

Tiransmitter  to  communication  receiver 
distance 

Path  Loss 

Hransmitter  to  intercqttion  detector  distance 

Path  Loss 

Table  1-1  Independent  Parameters  in  Background  and  Channel  Modules 
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ICPIC  includes  three  zero  mean,  white  noise  sources:  Gaussian,  Gamma  and  Cauchy.  The 
variance  of  discrete  samples  of  the  complex  envelope  of  noise  depends  on  the  sampling  frequency 
and  power  spectral  density  of  the  real  noise,  as  shown  in  Appendix  A.  The  Gaussian  and  Gamma 
noise  modules  are  parametrized  by  the  input  SNR  (Signal-To-Noise)  ratio  under  the  assumption 
that  the  two  modules  are  used  as  alternatives,  and  not  together. 

The  Gamma  Noise  module  permits  simulation  of  non-Gaussian  effects  and  produces 
radically  symmetric  complex  noise  whose  magnitude  is  a  gamma  variable.  The  marginal 
distributions  of  the  real  and  imaginary  components  are  not  expressible  as  elementary  functions,  but 
decay  exponentially.  The  only  independent  parameter  is  the  input  SNR  and  dB. 

Cauchy  noise  represents  heavy-tailed  noise  distributions  and  is  an  exceptional  case, 
because  it  has  infinite  variance.  The  Cauchy  Noise  module  produces  radially  symmetric  complex 
noise  whose  real  and  imaginary  components  are  marginally  distributed  as  Cauchy  random 
variables.  The  user  may  set  a  scaling  parameter  which  multiplies  a  standard  unit  variable.  The 
Cauchy  variable  has  an  undefined  mean  and  infinite  average  power,  so  that  signal  to  noise  ratio  has 
no  meting.  Instead,  Cauchy  Noise  is  parameterized  by  a  scaling  factor. 

The  interference  module  provides  the  user  with  three  options  with  which  they  may 
configure  the  module  with  an  arbitrary  number  of  interferers,  with  variable  widths,  powers  and 
frequencies. 

The  Random  option  allows  the  user  to  specify  the  number  of  interferers,  the  mean  and 
variance  of  the  interferer  powers,  and  the  maximum  and  minimum  interferer  widths,  and  the 
maximum  and  minimum  interferer  carrier  frequencies.  A  configuration  of  the  interferers  is  then 
randomly  generated  from  the  specified  distributions.  The  frequency  and  bandwidth  of  each 
interferer  are  generated  from  uniform  distributions,  and  the  power  of  the  interferers  are  generated 
from  a  log-log  distribution  on  analysis  of  empirical  observations  as  given  in  Appendix  B.  The 
interferer  width  is  implemented  by  BPSK  (Binary  Phase  Shift  Keying)  modulating  each 
interference  carrier  at  a  modulation  frequency  equal  to  half  the  interferer  width. 

The  Custom  option  allows  the  user  to  directly  specify  the  frequencies,  widths,  and  powers 
of  the  interferers.  The  user  is  given  the  capability  to  specify  the  number  of  interferers;  the 
configuration  center  frequency;  the  interferer  center  frequencies,  expressed  as  offsets  from  the 
configuration  center  frequency;  the  interferer  powers,  in  dB  above  the  communication  signal 
power,  and  the  widths  of  the  interferers.  The  bandwidth  is  implemented  by  BPSK  modulating 
each  interference  carrier  at  a  modulation  frequency  equal  to  half  the  bandwidth. 
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The  HF  and  VHF  options  allow  the  user  to  select  stored  interference  configurations  typical 
of  the  HF  and  weak,  medium  and  strong  VHF  bands.  An  interference  background  generated  from 
a  stored  configuration  chosen  to  be  characteristic  of  typical  operating  environments  is  added  to  the 
signal.  The  user  selects  the  center  and  width  of  the  operating  band,  and  the  modules  generate  the 
interferers  which  lie  in  the  selected  band. 

The  channel  category  models  phenomena  which  delay  or  distort  the  signal  or  change  its 
power,  such  as  multipath,  propagation  losses,  and  Doppler  shifts. 

A  custom  coded  tapped  delay  line  provides  a  foundation  for  channel  modeling.  The  ou^ut 
is  a  linear  combination  of  delayed  versions  of  the  input  with  constant  weights  specified  by  the  user. 
The  independent  parameter  is  the  total  delay  spread.  The  number  of  taps  is  equal  to  the  product  of 
the  delay  spread  and  the  sampling  frequency,  and  the  user  may  specify  arbitrary  tap  weights. 

If  the  communication  receiver  and  interferer  are  at  different  distances  from  the  transmitter, 
the  interferer  will  have  a  range  advantage  or  disadvantage  relative  to  the  communicator.  The  Path 
Loss  module  models  the  loss  with  a  power  law  decay  of  signal  power  as  a  function  of  propagation 
distance. 

1.4  Interaction  Between  Compression,  EDAC  and  Channel  Model 

EDAC  techniques  introduce  delays  due  to  processing,  retransmission,  and  associated 
buffering,  and  each  of  these  effects  can  have  important  system  design  inqjlications.  The  designer 
must  take  into  account  how  long  an  average  transmission  will  take,  the  requirements  for  real-time 
image  display,  the  tolerance  on  errors,  transmit  protocol,  and  how  much  buffering  memory  is 
available.  EDAC  effectiveness  may  be  sensitive  to  the  portion  of  the  compressed  image  to  which  it 
is  applied  based  upon  the  specific  compression  algorithm  being  used.  For  example,  it  may  be 
necessary  to  completely  eliminate  errors  from  the  image  header.  The  errors  themselves  are  a 
function  of  the  noise  model  and  it’s  magnitude.  An  accurate  assessment  of  these  combined  effects 
is  necessary  in  order  to  characterize  system  performance  and  optimize  the  system  design.  This  is 
the  functionality  that  the  ICPIC  workstation  testbed  provides. 
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2.0  An  Example  Of  An  End>To-End  Simulation 

This  section  describes  how  to  set  up  an  ICPIC  simulation.  The  process  is  divided  into 
SPW  and  non-SPW  steps.  The  SPW  steps  involve  building  the  actual  simulation  using  the  Block 
Diagram  Editor  (BDE)  of  SPW,  and  the  non-SPW  steps  generate  the  input  and  analyze  the  output 
signal.  The  tradeoff  between  using  the  custom  coded  modules  of  the  ICPIC  environment  or  off¬ 
line  analysis  modules  is  discussed  in  Section  2.5. 

2.1  Building  The  Simulation 


Building  a  SPW  simulation  requires  care  in  accoimting  for  module-specific  requirements. 

It  is  often  possible  to  build  the  same  simulation  in  slightly  different  ways  that  differ  drastically  in 
computational  efficiency.  Consider  the  simulation  shown  in  Figure  1,  which  models  the 
transmission  of  an  image  using  BPSK  modulation,  BCH  encoding,  with  the  addition  of  White 
Gaussian  Noise.  This  diagram  is  divided  into  eight  sub-blocks  for  the  purposes  of  this  discussion. 
We  will  discuss  each  block  and  (where  appropriate)  possible  alternatives  and  tradeoffs. 


Libniy/Hle 

V«c_iiglib^»e Jmim  « imig* 


Figure  1:  BPSK  BCH  Communication  Model 
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Rlnck  1:  Signal  Source 

This  defines  the  input  signal,  which  is  in  SPW  fonnat.  This  is  defined  by  a  12-line  header 
and  data,  an  example  of  which  is  shown  in  Figure  2.  The  details  of  this  file  and  how  it  is 
generated  from  an  image  will  be  discussed  in  Section  2.2. 

The  Signal  Source  block,  besides  defining  the  input  signal,  displays  it  in  the  Signal 
Calculator  module  of  the  SPW. 


$SIGNAL_FILE  9 
$USER_COMMENT 
$COMMON_lNFO 
SPW  Version  =  3.10 
System  Type  =  solaris2 
Sampling  Frequency  =  16 
Starting  Time  =  0 
$DATA_INFO 
Number  of  points  =  9960 
Signal  Type  =  Double 
$DATA 
1 
1 
0 
1 


Figure  2:  Sample  SPW  File 
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Block  2:  Vector  BCH  Encoder 

In  this  simulation  we  used  a  word  size  of  N=63,  and  a  data  length  of  k=30,  to  give  k/N 
close  to  1/2.  Contrary  to  the  SPW  documentation,  the  allowable  values  of  k  and  N  are  not 
arbitrary.  Figure  3  gives  the  allowable  (k,  N)  pairs^. 

Since  each  block  of  30  data  points  is  being  encoded  into  a  63-data  length  word,  SPW 
requires  that  the  Vector  BCH  Encoder  block  be  fed  in  30-data  chunks.  This  is  achieved  with  a 
Scaler  To  Vector  block,  which  takes  every  30  data  points  and  transforms  it  into  a  30  x  1  vector. 
The  output  from  the  BCH  encoder  is  a  63  x  1  vector  which  is  converted  back  into  set  of  63  scaler 
points  with  a  Vector  To  Scaler  module. 


N 

7 

15 

31 

63 


K 

4 

5, 7,11 

6,  11,  16,21,22 

7,  10, 16,  18,  24,  30,  36,  39,  45,  51,  57 


Figure  3:  Allowable  (N,  k)  Pairings  For  BCH  Module 
Alternative  Block  2:  Scaler  BCH  Encoder 

This  simulation  originally  attempted  to  use  a  scaler  BCH  encoder.  This  required  that  each 
data  point  be  repeated  63  times.  The  result  was  an  unacceptably  long  run  time  and  storage 
requirements  that  were  quickly  exceeded  by  images  files  larger  than  25  kbytes.  These  severe 
restrictions  were  removed  by  using  the  Vector  BCH  Encoder 

It  should  be  mentioned  that  the  vector  BCH  encoder ! decoder  modules  are  not  mentioned  in 
the  SPW  Communication  Library  R^erence. 


Block  3:  Repeat  Block 

The  Repeat  block  takes  each  scaler  input  and  repeats  it  the  specified  number  of  times  (16  in 
this  case).  That  is,  if  the  sequence  { 1, 0, 1 }  is  input  into  a  Repeat  block  with  a  repeat  factor  of  3, 
the  output  would  be  {1,  1,  1,  0,  0,  0,  1,  1, 1).  Adding  a  Repeat  block  is  equivalent  to  sending 
each  data  point  a  repeat  factor  number  of  times.  This  adds  redundancy  to  the  data  and  decreases 
the  effect  of  channel  noise  by  allowing  it  to  be  averaged  out  by  Block  5.  It  was  found  that  for  S/N 
<  20,  BCH  encoding  was  not  sufficient,  by  itself,  to  remove  the  effects  of  noise. 


1  The  SPW  Communications  Library  Reference  states  that  N  is  2"^  -  1 ,  where  2  <  m  <  9 
and  0  <  k  <=  N  is  arbitrary  within  this  range.  In  fact,  k  and  N  are  restricted  to  the  discrete 
pairings  shown  in  Figure  3. 
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This  brings  up  an  important  point  about  transmission  frequency.  The  transmission 
firequency  is  not  an  explicit  input  parameter,  but  is  effectively  set  to  unity.  As  far  as  this  simulation 
is  concerned,  it  does  not  care  what  the  frequency  is,  as  there  is  no  frequency-dependent  module  in 
this  simulation.  If  we  were  to  model  the  transmission  frequency  (typically  9,600  baud)  explicitly, 
then  each  data  point  would  have  to  be  Repeated  with  a  factor  equal  to  the  frequency.  This  would 
require  prohibitively  long  simulation  times  and  disk  space  requirements.  We  were  therefore  forced 
to  model  equivalent  systems  with  manageable  fiequencies. 

If  the  system  being  modeled  operates  at  a  frequency  of  v_Hz  and  a  repeat  factor  of  R,  then 
the  effective  frequency  is  v/_R  Hz  since  each  data  point  needs  to  be  transmitted  R  times. 

Block  4:  Modulator/Channel/De-Modulator 

The  encoded  signal  (consisting  of  zeros  and  ones)  is  BPSK-modulated  and  propagated 
through  a  channel.  In  this  case,  the  channel  is  modeled  as  additive  white  Gaussian  noise,  with  the 
addition  performed  by  the  T.  module.  Once  propagated,  the  signal  is  then  de-modulated  and  fed 
into  an  averager.  The  BPSK  modulator/demodulator  and  the  White  Gaussian  Noise  blocks  are 
custom-coded  routines.2 

Block  5:  Noise  Removal  Bv  Averaging 

This  block  takes  each  segment  of  R-repeated  data  and  averages  it  to  get  a  single  number. 
Remember  that  each  block  of  R  data  points  ouq)ut  from  the  decoder  represents  a  single  number,  if 
noise  were  not  present,  then  all  these  numbers  would  be  identical.  The  average  is  performed  by 
converting  each  block  of  R  numbers  into  a  vector,  multiplying  the  vector  by  1/R  (with  the  aX 
block),  and  then  summing  the  elements  of  the  vector  (Zx  block).  The  constant  1/R  is  stored  in  the 
Scaler  Constant  block  (C). 

Block  6:  Decision  Block 

Since  the  data  consists  of  zeros  and  ones  (as  dictated  by  BPSK  modulation),  a  check  is 
performed  to  determine  if  the  output  from  the  averager  exceeds  some  threshold.  This  threshold  is 
stored  in  another  constant  block.  If  the  average  value  exceeds  this  threshold  then  a  one  is  output,  if 
not,  a  zero  is  output.  The  threshold  generally  increases  with  decreasing  S/N  since  noise  in  this 
simulation  can  only  increase  the  signal  strength.  The  actual  value  of  the  threshold  is  determined  by 
trial  and  error. 


2  Curiously,  the  SPW  Communications  Library  has  its  own  BPSK  modulator,  but  not  a 
BPSK  demodulator  module. 
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Block  7:  Delay 

The  BPSK  modulator/demodulator  modules  (Block  4)  introduce  a  time  delay  into  the 
signal.  This  must  be  accounted  for  in  order  that  the  BCH  decoder  (Block  8)  be  properly 
synchronized  with  the  signal.  Otherwise,  the  decoder  will  operate  on  segments  from  two  separate 
signal  code  words  (of  size  N),  rather  than  a  single  code  word.  The  size  of  the  delay  is  N-1. 

The  delay  can  be  seen  by  placing  Signal  Sinks  before  the  BPSK  modulator  and  after  the 
BPSK  demodulator  (before  and  after  Block  4)  with  a  large  S/N  value  in  the  White  Noise  module 
(or  equivalently  with  the  noise  module  disconnected).  The  two  signals  will  appear  identical,  but 
phase-shifted.  This  test  should  always  be  done  before  using  any  ICPIC  system  configuration  in 
order  to  insure  an  accurate  baseline. 

Block  8:  Vector  BCH  Decoder 

This  block  is  analogous  to  Block  2,  except  the  numbers  are  reversed  since  the  signal  is 
being  decoded.  The  same  comments  mentioned  in  Alternative  Block  2  apply  here  as  welL 

Block  9:  Signal  Sink 

The  Signal  Sink  is  where  the  received  signal  is  stored.  From  the  Signal  Calculator  module 
of  the  SPW  the  received  and  input  signals  can  be  viewed  and  manipulated.  For  example,  the 
difference  signal,  mean,  error,  and  bit  error  rate  can  be  determined  using  the  functionality  of  the 
Signal  Calculator. 

2.2  Constructing  The  Input  Signal 

The  input  signal  is  constrocted  from  an  image  file  which  can  be  in  any  format  The  basic 
procedure  is  to  read  the  image  file,  one  bit  at  a  time,  and  convert  each  bit  into  a  binary 
representation.  Since  a  bit  is  8  bytes,  each  bit  is  converted  into  an  8-digit  binary  number.  The 
program  that  does  this  is  called  build_spw  and  is  invoked  by  typing: 

build_spw  image_file 

where  imagejile  is  the  name  of  the  image.  The  output  is  an  ASCII  file  called  spw_in.asc, 
which  consists  of  a  12-line  header  followed  by  the  binary-translated  data.  Each  digit  of  each  data 
point  is  placed  in  column  1  of  each  line  after  the  header,  as  illustrated  in  Figure  2. 
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Once  this  file  is  created,  it  needs  to  be  transferred  to  the  appropriate  SPW  directory,  which 
is  called  local _par^altadata/spwdata/signal_lib/,  where  local _path  is  the  location  of  the  top-level 
SPW  directory,  and  signal  Jib  is  the  (arbitrary)  name  of  the  signal  library  where  the  user  SPW 
signal  files  are  stored.  An  example  of  local is  /export/home/delmedico/.  This  transfer  is 
accomplished  by  typing: 

transfer 

from  the  C-shell  command  line,  transfer  is  a  script  that  takes  spw_in.asc  and  copies  it  to  a  file 
called  asc_input.asc  in  the  appropriate  SPW  directory.  This  script  needs  to  be  edited  before 
using  it  the  first  time  in  order  to  modify  local  jjoth  for  your  specific  system. 

2.3  Wavelet  To  SPW  File  Format  Conversion 

Since  wavelet  compression  is  not  a  standardized  technique,  and  there  are  many  different 
wavelet  routines  available,  transforming  a  wavelet  file  to  SPW  format  is  a  special  case  since  we 
have  to  create  the  wavelet  file  in  the  first  place.  After  much  searching,  we  found  the  best  package 
to  use  was  DSW;  this  is  a  descendant  of  the  David  Samoff  Labs  wavelet  package  and  is  used  by 
many  government  and  private  organizations.  A  limitation  of  DSW  is  that  it  only  converts  TGA 
image  files  to  a  wavelet  representation.  TGA  is  one  of  the  earlier  image  storage  formats  that  has 
the  benefit  of  being  easy  to  decode.  To  remove  this  limitation,  we  used  the  publicly-available 
pbmplus  graphics  file  conversion  package  which  allows  conversion  among  most  graphics  file 
formats,  pbmplus  consists  of  a  set  of  routines  to  convert  from  one  format  to  another  using 
intermediary  formats.  Using  bold  to  represent  executable  programs,  and  unbold  to  represent  file 
formats,  the  procedure  for  converting  from  GIF  to  wavelet  and  finally  to  SPW  is: 


GIF  ==>giftoppm  ==>PPM  =>ppmtotga  — >tga  ==>DSW  =>wavelet  file  ==>build  =>SPW  . 


Alternatively,  the  starting  point  could  have  been  TIFF  or  BMP,  or  another  graphics  format 
other  than  GIF.  pbmplus  does  not  have  any  JPEG  conversion  routines,  so  in  this  case  we  would 
use  another  public-domain  program  called  Xview,  which  is  invoked  by  typing  xv  at  the  command 
prompt.  From  here  a  JPEG  file  can  be  opened  and  saved  in  a  GIF  format  and  then  converted  using 
gifZwav.  This  entire  process  can  be  invoked  with  a  script  by  typing: 

gifZwav  gif_filename 

where  gifjilename  is  the  name  of  the  GIF  file  with  the  extension  .gif  assumed. 
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2.4  Data  Analysis 


The  output  signal  can  be  analyzed  within  the  Signal  Calculator  or  from  the  UNIX 
command-line.  From  within  the  Signal  Calculator  you  can; 

1.  Visually  compare  the  input  and  ou^ut  signals 

2.  Calculate  and  view  the  difference  signal 

3.  Calculate  statistics  on  the  difference  signal,  such  as  mean,  variance, 
and  number  of  errors. 

From  the  command  line,  a  tool  is  provided  to  calculate  all  the  quantities  in  Item  3  above, 
with  the  added  convenience  of  batch  mode  processing.  SPW  provides  a  mechanism  for  batch 
mode  processing  a  simulation  when  only  one  input  parameter  changes.  In  this  way  a  user  can 
parameterize  the  output  statistics  as  a  function  of  S/N,  with  all  other  variables  kept  constant  Each 
output  is  automatically  stored  in  a  separate  file  which  can  be  parsed  and  fed  into  the  command-line 
analysis  code. 

The  code  that  calculates  the  output  signal  statistics  is  called  rms  and  is  invoked  from  the 
command  line  by 

rms  spw _output_filename 

where  spw_output_filencane  is  the  output  from  SPW  that  is  stored  by  Block  8  in  Section  2.1. 

The  output  from  rms  will  be  discussed  in  more  detail  in  Section  3. 

2.5  Reconstructing  And  Viewing  The  Transmitted  Image 

Although  not  a  quantitative  measure,  the  most  intuitive  way  of  judging  how  well  the 
transmitted  signal  was  received  is  to  view  the  transmitted  and  received  images  side  by  side.  This 
requires  that  the  raw  data  from  the  Signal  Sink  (Block  8  in  Section  2.1)  be  reconstructed  to  form  an 
image  file.  This  is  accomplished  by  typing: 

getdata 

from  the  UNIX  command  line;  this  takes  the  ICPIC/SPW  output  file  from  the  SPW  signal  library 
directory  and  places  it  in  the  user’s  working  directory.  As  with  the  transfer  script  (Section  2.2), 
getdata  needs  to  be  edited  before  using  it  the  firet  time  in  order  to  change  the  relevant  path  names, 
getdata  takes  the  outfile  file  called  real.out  and  copies  it  to  the  working  directory  as  spw.out. 


17 


The  data  in  spw.out  is  then  reconstructed  into  an  image  file  by  typing: 


build 

from  the  UNIX  command  line,  build  takes  each  block  of  8  lines  past  the  12-line  header  file, 
which  consists  of  zeros  or  ones3,  treats  it  as  an  8-bit  binary  number,  and  converts  it  back  to  an  8- 
bit  decimal  integer  (0  -  255). 

The  analysis  program  (rms)  uses  the  same  procedure  to  convert  fiwm  binary  to 
numbers.  To  do  otherwise  would  ignore  the  in^licit  weighting  of  the  binary  digits.  The  outputs 
fix)m  rms  are  the  number  of  bit  errors  (this  is  based  on  the  binary  representation),  the  RMS  error 
between  the  original  image  and  the  output  fiom  ICPIC/SPW,  and  a  quantity  called  the  fv^ntrast 
Error.  This  is  basically  a  normalized  RMS,  where  the  normalization  is  the  average  pixel  value  of 
the  image. 

Once  the  image  file  is  reconstructed,  it  can  be  displayed  using  xv.  In  some  visions  of  xv, 
it  is  possible  to  display  multiple  files  side  by  side.  In  other  versions  only  a  single  file  can  be 
displayed;  in  this  case  xv  must  be  invoked  twice,  once  for  the  original  and  once  for  the 
reconstructed  image.  If  the  S/N  ratio  of  the  Gaussian  noise  is  too  low,  then  xv  may  not  recognize 
the  image  as  a  valid  format  (it  can  handle  most  of  them)  and  will  not  be  able  to  display  it 

3.0  Sample  Output 

The  GIF  image  shown  in  Figure  4  was  used  to  study  the  effects  of  varying  the  encoder 
parameters  (k,  N)  and  S/N.  It  has  sufficient  detail  to  show  any  deviations  from  the  original  and  is 
small  enough  (19.2  Kb)  to  allow  quick  turn  around  times  (under  10  minutes/run). 


Figure  4:  Book  graphic  test  image 


3  This  is  guaranteed  by  the  Decision  Block  (Biock  6  of  Section  2.1). 
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CONTRAST  ERROR 


Figure  5  shows  plots  of  the  Contrast  Error  for  three  (k,  N)  values.  For  all  cases  the 
transmitted  image  was  identical  to  the  original  for  S/N  >  15  dB.  For  7/63  coding  (7  information 
bits  in  a  code  word  of  63  bits),  the  transmitted  image  is  identical  to  the  original  for  S/N  as  low  as  5 
dB. 


BCH  Encoding 
BPSK  Modulation 


S/N  (db) 


Figure  5:  Variation  Of  Contrast  Error  With  S/N 


Notice  that  as  S/N  tends  to  zero,  the  error  does  not  tend  to  infinity.  The  reason  is  that  the 
decision  block  (Block  6)  forces  the  input  to  the  decoder  to  be  zero  or  one.  As  such,  as  S/N 
decreases,  the  output  saturates  to  all  ones  (totally  black  image).  This  represents  the  maximum 
deviatim  from  the  original  image,  and  also  points  out  that  this  deviation  is  image-dependent  For 
exan^le,  if  tiie  image  had  been  mostly  black,  then  the  contrast  error  would  have  been  smaller  for 
die  same  system  parameters.  This  same  behavior  was  seen  when  looking  at  the  fraction  of 
inconect  bits  in  die  received  image. 
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Figure  6  shows  the  signal  before  and  after  Block  4,  respectively,  for  the  case  of  no  noise. 
This  illustrates  the  phase  shift  (([))  introduced  by  the  BPSK  modulator/demodulators,  and  illustrates 

the  need  for  the  delay  module  (Block  7).  The  delay  is  set  to  N  -  (()  sample  points.  If  the  delay 

block  is  not  added  to  the  system,  then  the  BCH  decoder  will  receive  segments  from  two 
consecutive  code  word  blocks,  rather  than  a  single  code  block  of  63  bits.  An  indication  of  an 
incorrect  delay  value  is  a  deviation  between  the  transmitted  and  original  image  when  no  noise  is 
present  This  test  case  should  always  be  performed  before  attempting  any  tradeoff  studies. 

Figure  7  shows  the  signal  before  and  after  Block  4  with  S/N  =  5.  The  two  signals  are  now 
clearly  different,  but  all  the  deviations  (except  for  the  phase  shift)  are  accounted  for  by  the  BCH 
decoder  for  k=7  and  N=63. 


lib/BeforeBod.sig  . 

1201.851 

11904 

92 

0.077  sec  \ 
0  . 


Figure:  6:  Signal  before  (top)  and  after  BPSK  Modulator  block  vitbont  noise. 
The  signals  are  l^nticle,  exc^t  for  a  phase  shift. 


Figure:  7:  Same  as  above,  but  with  noise  added  (S/N  ■•5}. 

The  signals  are  now  different,  but  this  is  corrected  by  the  BCH  decoder  for  K-7  and  H-63. 
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4.0  Summary 


ICPIC  provides  a  flexible  environment  that  gives  the  user  the  capability  of  modeling 
various  communication  channels,  error  detection  and  correction  techniques  (EDAC)  and 
modulation  types,  as  well  as  the  ability  to  assess  their  effect  on  imagery  compressed  with  various 
algorithms.  Software  tools  are  provided  to  convert  any  type  of  image  (including  wavelet)  into  a 
SPW  format,  and  to  convert  from  one  image  format  to  another.  Changing  the  channel  and  EDAC 
models  is  accomplished  by  dragging  and  dropping  operations  in  the  BDE  (Block  Diagram  Editor) 
and  then  setting  the  appropriate  parameters  for  each  module. 

The  image  preparation  and  analysis  tools  were  developed  off-line  from  ICPIC/SPW  to  run 
from  the  UNIX  command  line.  This  was  found  to  be  a  much  more  efficient  approach  then  to  force 
everything  within  the  ICPIC/SPW  framework.  Included  in  these  off-line  software  tools  is  the 
capability  to  reconstruct  the  received  image  from  the  ICPIC/SPW  output  Depending  on  the  errors 
in  the  received  image,  graphics  viewers  (such  as  xv)  may  not  be  able  to  display  the  reconstructed 
image.  This  is  especially  true  if  errors  occur  in  the  header  file.  One  approach  to  minimize  header 
errors  is  to  agree  on  a  particular  fixed  header,  and  then  only  send  the  actual  data.  It  is  also  possible 
to  develop  robust  image  viewers  that  attempt  to  correct  for  corrupt  data.  These  issues  should  be 
investigated  in  the  future. 

To  date  ICPIC’s  analysis  capabilities  have  been  applied  to  simplex  communications 
channels.  These  are  one-way  charmels  where  the  receiver  does  not  have  the  capability  to  request 
re-transmission  of  data  missed  due  to  channel  bit-errors.  Simplex  channels  are  primarily  used 
today  in  the  tactical  arena.  These  are  situations  in  which  the  receiver  does  not  want  to  turn  on  a 
transmitter  out  of  fear  of  detection  by  hostile  forces  or  alternatively  for  portability  reasons  the 
receiver  sight  lacks  sufficient  power  for  transmission.  UHF  satellite  channels,  such  as  FLTS  AT, 
are  a  good  example  of  channels  used  primarily  for  simplex  communications.  The  application  of  the 
developed  capability  to  evaluate  the  performance  on  compressed  imagery  of  channels  such  as  these 
could  also  be  an  area  for  future  investigation. 

Another  area  for  investigation  is  how  to  specify  quality  of  service  and  to  best  evaluate  the 
results  of  chaimel  induced  corruption  on  compressed  imagery.  Quantitative  comparison  (RMS 
error,  etc.)  of  the  reconstructed  corrupted  images  and  original  images  can  often  be  misleading. 
Qualitative  con^arison  by  side-by-side  viewing  of  the  reconstructed  corrupted  image  and  original 
image  is  often  more  useful  but  is  subjective.  The  basic  problem  is  that  of  objectifying  a  basically 
subjective  phenomena.  In  addition,  the  performance  depends  on  the  particular  scenario.  For 
example,  is  the  goal  to  pick  out  a  rocket  launcher  in  the  woods,  or  to  distinguish  between  classes 
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of  tanks?  It  may  be  acceptable  to  use  much  higher  S/N  ratios  to  answer  the  first  question  and  this 
may  effect  the  final  choice  of  system  operating  parameters. 

The  fact  that  extremely  corrupted  images  cannot  be  decompressed  at  all  for  viewing  leads  to 
yet  another  area  for  future  investigation.  Since  an  image  that  cannot  be  decon:q)ressed  is  useless,  it 
is  desirable  to  conduct  more  detailed  tradeoff  studies  with  more  complicated  channel  models,  better 
imagery  and  different  compression  algorithms  in  order  to  investigate  the  error  threshold  for 
“decompressibility”  for  different  EDAC  methods  and  compression  algoridims.  This  could  be 
accomplished  for  various  channels  and  modulation  types.  Hie  goal  would  be  to  document  which 
EDAC  methods,  modulation  types  and  compression  algorithms  are  best  suited  to  getting  a 
deconqiressible  (and  hopefully  viewable)  image  through  under  various  channel  conditioils. 

These  are  significant  issues  that  need  to  be  fully  addressed,  and  ICPIC/SPW  provides  the 
fiamework  to  address  them. 
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Appendix  A 

Complex,  Sampled,  White  Noise 


Wavefram  simulations  generally  process  discrete  sanq)les  of  continuous  waveforms.  In 
die  case  of  white  noise  a  continuous  wavefcom  is  not  avail^le  to  sample,  so  random  discrete 
sequences  are  generated  directly,  and  a  method  is  needed  to  set  die  parameters  of  the  ^nerating 
probability  distribution.  Ibe  chosen  mediod  is  to  r^uire  the  discrete  simuladon  of  traimtted 
power  through  any  filter  to  mateh  exacdy  the  transmitted  power  in  the  corresponding  continuous 
process. 

Conqilex  envelopes  are  processed  so  that  real  white  noise,  which  c^ot  be  expressed  in 
tPTms  of  quadrature  components,  are  modeled  as  real  bandpass  white  noise  n  with  a  power 
spectral  density  given  by 

So  ||/|-/c«|<^/2 

-S.(/)  =  (  (A-l) 

0  otherwise 

where/ggg  is  the  signal  center  frequency,  and  B  is  much  larger  than  the  largest  bandwidth  being 
simulated,  but  less  than  2/  ggg  .  The  con^lex  envelope  z(t)  of  n(t)  relative  to/  ggg  has  a  power 
spectral  density  given  as 


\f\<B/2 

otherwise 


(A-2) 


and  the  quadrature  components  X  andy  have  identical  power  spectra: 

S,(f)  =  Sy(f)  =  S.(f)  (A-3) 


Farh  quadrature  cfwnponent  is  modeled  by  a  discrete  random  sequence  with  variance 
de»«rminftd  by  the  criterion  that  Ae  noise  power  transmitted  through  any  r^  lowpass  filter  acting 
as  tiie  noise  baseband  can  be  accuratdy  modeled  by  the  discrete  simulatioiL  die  quadrature 
rfwnpnnftnt  x  is  passed  tiuough  a  filter  with  transfer  function  H(f)  and  real  in^ulse  response 


function  h  (t),  the  transmitted  power  is: 


foo  > 

2' 

jx{t-T)h{T)dT 

► 

V-o  J 

(A-4) 


A-1 


Expanding  the  square  as  the  product  of  two  integrals  and  taking  the  expectation  of  the  integrand 
gives: 

oo  oo 

^=J  !‘^x(-f2-ti)>‘{ti)h(r2)dridT2  (A-5) 

— OO  *-00 

where  (t)  is  the  autocorrelation  function  of  x  Expressing  (t)  as  the  Fourier  transform  of 
and  using  (A-2)  and  (A-3)  gives: 

2 

/>  =  25o||H(/)|  df  (A.6) 

Li  a  discrete  simulation  (A-5)  becomes: 

P^screte  =fs'a'LPx  ^  ^  h  (A-T) 

ijc  V  Jsa  y  \Jsa  y  \Jsa  j 

which  simplifies,  since  the  discrete  random  numbers  are  uncorrelated,  to 

Pdixcre,e=ft(P\Y.hH^  (A-8) 

i  \fsa  j 

where 

<^l=Px(0)  (A-9) 

is  the  variance  of  the  discrete  time  sequence  modeling  the  x-component  For  large  sample 
ftequency  the  sum  in  (A-8)  approximates  an  integral,  giving: 

oo  2 

^discrete  —fsa  (A- 10) 

— OO 

SO  using  the  norm  preserving  property  of  Fourier  transforms  gives: 

2 

Pdiscrete  =  J  W  (t)|  df  (A-11) 
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Comparing  (A-1 1)  to  (A-6)  and  requiring  the  discrete  simulation  to  model  the  continuous  process 
accurately  means  the  variance  of  the  discrete  quadrature  component  must  be  given  by: 


(A-12) 

while  (A-3)  shows  both  quadrature  components  have  the  same  power  so, 

al=al=2f^S^  (A-13) 

The  intended  communicator’s  input  signal-to-noise  ratio  is  a  dimensionless  quantity  that  is 
of  greater  interest  to  the  user  than  noise  variances  of  power  spectral  densities.  For  commuiucation 
signal  power  equal  to  unity  shows  the  input  signal-to-noise  ratio  is  given  by: 


(SIN)i  = 


fda^O 


Combining  (A-13)  with  (A- 14)  gives: 


_^2  _ 
(Jjj  Gy 


2/. 


sa 


fda(S/N), 


(A-14) 


(A-15) 


which  is  built  into  the  white  Gaussian  noise  block  so  the  independent  parameter  set  by  the  user  is 
the  communicator’s  input  signal-to-noise  ratio. 
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Appendix  B 

Interference  Power  Distributions 

Literferaice  environments  change  dramatically  with  geographical  location,  season, 
time  of  day,  and  frequency  band,  and  while  trrads  have  been  established  [1],  many  of  the 
details  of  die  environments  remain  unpredictable.  To  model  die  niolt^licity  of 
unpredictable  environments  a  Random  option  is  available  which  randomly  gm^rates  a 
cc^guration  of  inteiferers,  subject  to  probability  ^tiibudons  whose  parameters  are 
specified  1>y  die  user.  Enseal  observations  repmted  and  analyzed  suggest  the  use  of  a 
specific  distribution  of  inte^erer  powers.  Sin^lifying  die  distribution  and  paraitv»«»riTing 
in  terms  of  die  mean  and  variance  of  the  interferer  powers  enhances  the  Random  option’s 
ease  of  use. 

The  probability  distribution  used  is 


f(p)= 


p<p„M 

Ptom  ^  Ptaax. 


P  ^  Pmax 


(B-1) 


and  gives  the  probability  that  a  randomly  chosen  interference  line  will  have  power  less  than 
p.  Interferer  powers  are  measured  in  units  of  the  signal  power  at  Ae  cnTnmnninntinn 
receiver,  sop  is  the  ratio  of  interferer  power  to  signal  power.  For  simplicity  PmOTwill  be 
taken  to  approach  infinity,  and  will  be  chosen  so  that  F(p)  qiproaches  zero 

contmuously  at  P^,  which  requires  1  -  Cp^  =  0,  so 

^uiin=C‘'‘'  (B-2) 


Difierentiating  (B-1)  gives  the  probability  density  as 


P<Pmin 
P  ^  Pmin 


(B-3) 


The  mean  power  is  given  by 


lpf(p)dp 

Paia 


(B-4) 


and  using  (B-3)  and  evaluating  the  integral  gives: 


Cd 

d-l 


(B-5) 


B-1 


provided 


(B-6) 


otherwise  the  density  in  (B-3)  has  infinite  mean.  The  variance  is  given  by 

<^^  =  ]  P^f  (P)  ^P  " 

Pmin 


and  using  (B-3)  and  evaluating  the  integral  gives: 


^2  Cd  2-d  ,,2 

"d-2^  ^ 

(B-8) 

provided 

d-2>0 

(B-9) 

otherwise  the  density  in  (B-3)  has  infimte  variance.  Using  (B-2)  to  eliminate 
(B-5)  and  (B-8)  gives 

tom 

U=  ‘‘  C''‘‘ 

(B-10) 

d-\ 

and 

d-2 

Solving  for  the  parameters  C  and  d  in  terms  of  the  mean  and  variance  proceeds  by 
using  (B-10)  to  eliminate  C  tom  (B-11)  giving: 

d^-2d-i^  =  0  (B-12) 

(7^ 


and  taking  the  positive  root  gives: 


and  rearranging  (B-10)  gives: 


C  = 


(B-14) 


Q>mparing  (B-13)  with  (B-9)  shows  that  as  the  standard  deviation  of  the  distribution 

becomes  much  larger  than  the  mean  JJ,  f (7  —>0,d  —>2  and  the  distribution 
approaches  a  singular  limit 

The  user  supplies  the  mean  power  Ep  and  the  power  standard  deviation  Jp  in  dB 
so 


=  10 


E,I\Q 


and 


Gr  =  10 


£p/10 


(B-15) 


(B-16) 


(B-15)  and  (B-16)  in  (B-13)  and  (B-14)  are  used  to  supply  the  coefficient  values  to  the 
distribution  in  (B-1),  and  then  randomly  generate  the  interferer  powers. 
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Appendix  C 
ICPIC  Modules 

The  following  list  of  entries  uses  a  format  based  on  the  Unix  man  pages,  this  format 
provides  easily  accessible  documentation  of  the  modules  included  in  the  ICPIC  workstation 
testbed.  The  modules  are  ordered  alphabetically  by  component  categories.  The  entry  for  each 
module  provides  a  functional  description  of  each  module,  the  path  leading  to  the  module  when 
using  the  SPW  module  selection  under  the  “logical  grouping”  option,  the  input  and  output  ports, 
and  the  parameters  in  the  module. 

There  are  two  fundamental  types  of  modules,  known  in  SPW  parlance  as  “custom-coded” 
modules  and  “multi-level”  modules.  The  former  refers  to  modules  developed  from  scratch  in  C 
code.  These  are  designated  in  the  documentation  by  “Coded  in  C”.  Modules  not  indicated  as  such 
are  multi-level.  These  are  typically  subsystems  built  from  lower-level  modules  such  as  SPW 
primitives,  custom-coded  modules,  or  even  lower  level  multi-level  modules. 
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CONTENTS 


BACKGROUND 

Cauchy  Noise 
Custom  Interference 
Gamma 

Random  Interference 


CHANNEL 

Communication  Coupler 
Intercept  Coupler 
Tapped  Delay  Line 


MODULATOR 

8PSK  modulator 
BPSK  modulator 
MSK  modulator 
QPSK  modulator 


DEMODULATOR 

8PSK  demodulator 
BPSK  demodulator 
FSK  demodulator 
MSK  demodulator 
QPSK  demodulator 


Overflow  Protect 
Phasor 

Random  Digits 
Read  From  File 
Square  Wave 
Stepper 
Store  Power 
Store  Signals 

Store  Three  Dimensional  Data 
Two  Dimensional  Stepper 
Vector  Magnitude  Squared 
Vector  Mean 
Vector  Subsampler 


SPREAD  SPECTRUM 


DETECTORS 


PREPROCESSORS 


Bandpass  Cyclic  Feature 
Correlation  Doubler 
Correlation  Radiometer 
Delay  Chip  Rate 
Envelope  Chip  Rate 
Delay,  FFT  Search,  Chip  Rate 
Envelope,  FFT  Search,  Chip  Rate 
Frequency  Doubler 
Lowpass  Cyclic  Feature 
Scanner 
Spectra 
Spectrum 

Etelay,  Stepping  Search,  Ch  Rate 
Envelope,  Stepping  Search,  Ch  Rate 
Store  Detector  Data 
Store  Spectra 
Store  Spectrum 
Switch^  Radiometer 
Total  Power  Radiometer 


Adaptive  Threshold  Excision 
Adaptive  Whitening  Filter 
Communicator’s  Composite 
Density  Based 
Fixed  Excision 
Frequency  Hop  Downconvert 
MMSE  Whitening  Filter 
Spectral  Estimate  and  Match 
\^tening  Excision 

UTILITIES 

Auto  Stepper 
Bandpass  Integrator 
Center  Frequency  Mismatch 
Chip  Release 
Complex  Input  Switch 
Complex  Moving  Average 
Complex  Variance 


8ary  DS  Spreader 
Agility 

Binary  DS  Spreader 
Burst 

Chip  Rate  Agile  8ary  DS 
Chip  Rate  Agile  Binary  DS 
Chip  Rate  Agile  4ary  DS 
Chip  Rate  Agility 
Direct  Sequence  Despreading 
Filtering 

Frequency  Agility 
Frequency  Hopping 
Jitter 

Jittered  8ary  DS  Spreader 
Jittered  Binary  DS  Spreader 
Jittered  4ary  DS  Spreader 
Qary  DS  Spreader 
Triggered  Frequency  Hopping 


Interceptor  >  Metrics 

Exceedance 
HistogramFile  Stepper 
Composite  Output  SNR 
ROC 

Ricean  ROC 
Simulated  Output  SNR 


D/A  Converter 
Deinterleaver 
File  Stepper 
Heterodyne 
Interleaver 

Leading  Edge  Mod  N  Impulse  Train 
Lowpass/Bandpass  Filter 
Magnitude  Squared 
Nyquist  Filtering 
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8Ary  DS  SPREADER 


DESCRIPTION 

Multiplies  the  complex  input  by  an  8-fold  complex  pseudonoise  sequence  with  a  chip  rate 
specified  by  the  user  as  a  multiple  of  the  data  rate.  If  the  input  signal  has  constant  unit 
magnitude,  the  output  power  is  1.  Delays  the  output  signal  by  the  input  delay  value  for 
automatic  synchronization.  Outputs  the  spreading  sequence  as  a  global  output  connector 
called  ‘code’. 

MENU  PATH 

Spread  Spectrum  >  Direct  Sequence  >  8ary 
INPUT  PORT 

IN  Complex 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

nch  Number  of  chips  per  symbol 

fda  Data  rate 

sfreq  Sampling  frequency 

nsamples  Number  of  samples  in  simulation 

LIBRARY  /  FUNCTION 
walib/8ary_ds 


C-3 


8PSK  DEMODULATOR 


DESCRIPTION 

An  8-ary  PSK  demodulator. 

MENU  PATH 

Communicator  >  Demodulator  >  8PSK 
INPUT  PORT 


IN 

Complex 

DELAY_IN 

Real 

OUTPUT  PORT 

OUT 

Real 

DELAY  OUT 

Real 

PARAMETERS 

fda 

sfreq 

LIBRARY  /  FUNCTION 
modIib/8psk-demod 


C-4 


8PSK  MODULATOR 


DESCRIPTION 

An  8-aiy  PSK  modulator.  A  binary  input  stream  is  gray  coded  onto  an  eight  phase  shift 
keyed  signal 

MENU  PATH 

Communicator  >  Modulator  >  8PSK 
INPUT  PORT 


IN 

Real 

DELAYJN 

Real 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

fda 

Bit  rate 

sfreq 

Sampling  rate 

LIBRARY  /  FUNCTION 
modlib/8PSK-mod 


C-5 


8PSK  SOURCE 


DESCRIPTION 

Outputs  a  narrowband  complex  envelope  resulting  from  interleaving  a  pseudorandom  bit 
sovuce,  and  using  the  result  as  the  data  stream  for  an  8-ary  phase  shift  keying  modidation. 
The  user  sets  the  numbers  of  columns  and  rows  to  define  the  interleaving,  and  accepting  the 
default  value  of  1  for  both  produces  a  non-interleaved  data  stream. 

MENU  PATH 

Communicator  >  Modulated  Source  >  8PSK 
INPUT  PORT 
None 

OUTPUT  PORT 

OUT  Complex 

DELAY_OUT  Real 

PARAMETERS 

ncols  Interleaver  columns 

n  rows  Interleaver  rows 

fda  Bit  rate 

sfreq  Sampling  rate 

LIBRARY  /  FUNCTION 
modlib/8psk-source 


C-6 


ADAPTIVE  THRESHOLD  EXCISION  CORE 


DESCRIPTION 

Combines  the  two  input  vectors  to  form  a  complex  vector,  and  removes  the  components 
with  squared  magnitude  greater  that  an  exceedance  factor  times  a  mean  square  magnitude. 
The  mean  square  magnitude  is  of  all  the  bins  except  a  number  of  strongest  bins.  The  user 
specifies  the  length  of  the  vectors,  the  number  of  strongest  bins  to  be  excluded  from  the 
calculation  of  mean  square  magnitude,  and  the  exceedance  factor.  The  spectral  resolution, 
or  bin  size  is  the  sampling  frequency  divided  by  the  dimensionality.  The  algorithm 
recomputes  the  threshold  for  each  successive  block  of  data  and  so  adapts  to  changing 
background  levels. 

MENU  PATH 

None 

INPUT  PORTS 


X_re 

X_im 

OUTPUT  PORT 
Y_re 
Y_im 

PARAMETERS 

blocksize 

nexempt 

exceedfactor 


Real  vector 
Real  vector 

Real  vector 
Real  vector 

Number  of  points  used  in  the  FFTs 

Number  of  points  excluded  in  the  calculation  of  mean  square 
magnitude 

Exceedance  factor 


window_type  Type  of  window  used  in  the  forward  FFT 

sfreq  Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 
prlib/exadapt_cmp 
Codedin  C 


C-7 


ADAPTIVE  WHITENING  FILTER 


DESCRIPTION 

Implements  a  non-overlapped  version  of  the  adaptive  whitening  matched  filter,  which 
serves  as  the  basis  for  the  mmse  adaptive  whitening  filter. 

MENU  PATH 

None 

INPUT  PORT 

IN  Complex 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

blocksize  Length  of  FFT 

average_type  Periodogram  averaging  type,  either  ‘sliding_block’  or  ‘decaying’ 

averaging_block  Length  of  sliding  window  for  ‘sliding_block’  averaging  type 

daniell_param  Width  of  frequency  domain  smoothing  window 

alpha  Decay  parameter  for  ‘decaying’  averaging  type 

window_type  Tapering  window  applied  to  input  blocks  to  provide  additional 

smoothing 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
prlib/adapt_whiten 


C-8 


AGILITY 


DESCRIPTION 

Produces  trigger  pulses  separated  in  time  by  random  intervals  called  blocks.  The  user 
specifies  the  upper  and  lower  bounds  of  a  uniform  distribution;  each  sample  of  the  uniform 
distribution  is  rounded  to  the  nearest  multiple  of  the  data  symlxil  duration;  and  the  result  is 
used  as  the  cxirrent  block  length.  Setting  an  integer  flag  to  1  stores  the  trigger  pulses  and 
block  durations  as  signals  for  subsequent  inspection  in  SDE. 

MENU  PATH 

Spread  Spectrum  >  Agility 

INPUT  PORT 

IN 

OUTPUT  PORT 
TRIGGER 
PARAMETERS 

min_duration 
max_duration 
sfreq 
fda 

store_agility 
LIBRARY  /  FUNCTION 
walib/agility 


None 

Complex 

Minimum  block  length 
Maximum  block  length 
Sampling  frequency 
Data  rate 

1  to  store  trigger  pulses  and  block  lengths,  0  otherwise 


C-9 


ANALOG  TO  VECTOR  CONVERSION 


DESCRIPTION 
Coded  in  C 
MENU  PATH 

Utilities  >  Vector  >  Analog  to  Vector 
INPUT  PORT 

IN  Real 

HOLD  Real 

OUTPUT  PORT 

OUT  Real  Vector 

PARAMETERS 

OUTJOVEC_LEN 

Output  vector  length 

LIBRARY  /  FUNCTION 
utlib/real_to_vector 


C-10 


ANNOUNCE 


DESCRIPTION 

Produces  a  message  in  a  pop-up  window  at  the  completion  of  a  simulation  run. 
MENU  PATH 

Utilities  >  Input/Output  >  Announce 
INPUT  PORT 

HOLD  Real 

OUTPUT  PORT 
None 

PARAMETERS 

control  Print  message  if  >  0 

niessage_l  First  message  line 

message__2  Second  message  line 

LIBRARY  /  FUNCTION 
utlib/announce 
Codedin  C 


C-11 


AUTO  STEPPER 


DESCRIPTION 

Produces  a  stepped  output  which  cycles  through  a  set  of  values  specified  by  the  user.  The 
values  are  re^arly  spaced  and  specified  by  the  starting  value,  the  increment,  and  the 
number  of  vdues.  The  user  controls  the  dwell  time,  which  is  Ae  time  interval  between 
changes  in  the  output. 

MENU  PATH 

Utilities  >  Sources  >  Auto  Stepper 
INPUT  PORT 

HOLD  Real 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 

dweil_time  Time  interval  between  changes  in  output 

start_freq  Initial  output  value 

step_freq  Step  in  output  value 

no_of_values  Number  of  output  values 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/auto_stepper 


C-12 


BANDPASS  CYCLIC  FEATURE  DETECTOR  SYSTEM 


DESCRIPTION 


Uses  a  communication  signal  input  and  a  background  input  to  form  a  signal  present  input  to 
one  detector  and  a  signal  absent  input  to  a  second  detector.  The  detectors  are  each  cyclic 
feature  detectors  for  bandpass  cyclic  frequencies.  The  result  of  processing  the  signal 
present  case  is  the  output  and  the  result  of  processing  the  signal  absent  case  in  ihe 
H 0  output.  Stores  detector  parameter  values  in  files. 

MENU  PATH 


Interceptor  >  Detector  >  Cyclic  Feature  >  Bandpass,  system 
INPUT  PORTS 


SIGNAL 

Complex 

BACKGROUND 

Complex 

OUTPUT  PORT 

HI 

Complex 

HO 

Complex 

PARAMETERS 

cyclic_offset 

fiIter_zone 

lp_edge 

bp_  center 

bp_bw 

delay_time 

colIect_time 

wswgn 

nwssb 


Offset  of  cyclic  frequency  from  twice  signal  center  frequency 
Spectral  frequency  zone:  0  for  lowpass,  1  for  bandpass 
Spectral  lowpass  edge  frequency 
Spectral  bandpass  center  frequency 
Spectral  and  pass  bandwidth 

Delay  in  delay  and  multiply  section  of  detector  systems 
Output  integration  time  in  detector  systems 
1  for  weak  signal,  white  Gaussian  noise  mode,  0  otherwise 
1  for  nonweak  signal,  simulated  background  mode,  0  otherwise 


C-13 


BANDPASS  CYCLIC  FEATURE  DETECTOR  SYSTEM,  continued 


sfreq 

nsamples 

fda 


Sampling  frequency 
Number  of  samples  in  simulation 
Data  rate 


LIBRARY  /  FUNCTION 
delib/cyc_feat2 


BANDPASS  CYCLIC  FEATURE  DETECTOR  CORE 


DESCRIPTION 


Simulates  the  action  of  a  cyclic  feature  detector  for  bandpass  cyclic  frequencies  with 
complex  envelope  input.  Provides  optional  storage  of  intermediate  signals. 

MENU  PATH 

Interceptor  > 

Detector  >  Cyclic  Feature  >  Bandpass,  core 

INPUT  PORT 

IN 

Complex 

HOLD 

Real 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

cyclic_offset 

Offset  of  cyclic  frequency  form  twice  signal  center  frequency 

filter_zone 

Spectral  frequency  zone:  0  for  lowpass,  1  for  bandpass 

Ip_edge 

Spectral  lowpass  edge  frequency 

bp_center 

Spectral  bandpass  center  frequency 

bp_bw 

Spectral  bandpass  bandwidth 

delay_time 

Delay  in  delay,  conjugate  and  multiply  section 

collect_time 

Output  integration  time 

sfreq 

Sampling  frequency 

nsamples 

Number  of  samples  in  simulation 

LIBRARY  /  FUNCTION 


deIib/cyc_feat2_core 


C-15 


BANDPASS  CYCLIC  SPECTRUM 


DESCRIPTION 

Computes  the  bandpass  component  of  the  cyclic  spectmm  of  a  complex  signal  in  SDE.  The 
cyclic  frequency  of  the  bandpass  cyclic  spectrum  computed  is  entered  as  an  offset  from 
twice  the  assumed  carrier  frequency.  The  spectrum  is  computed  for  spectral  frequencies 
between  H  and/^^  H  where is  the  sampling  frequency. 

MENU  PATH 

SDE:  SYSTEM/MACRO/EXECUTE/vmIic/Cyclic_bandpass 
PARAMETERS 

numsig  Display  number  of  signal  to  analyze 

alpha  Cyclic  frequency  desired 

samp_freq  Sampling  frequency 

numpts  Number  of  samples  to  use  in  the  computation 

numffpts  Number  of  points  in  the  FFT  (must  be  a  power  of  two) 

LIBRARY  /  FUNCTION 

spwdata/sdemacro/vmIib/CycIlc_high.mac 


C-16 


BANDPASS  INTEGRATOR 


DESCRIPTION 

Multiplies  the  input  by  a  complex  tone  with  frequency  equal  to  the  bandpass  integrator 
center  frequency,  and  averages  the  result  over  the  integration  time. 

MENU  PATH 

Utilities  >  Signal  Processing  >  Bandpass  Integrator 
INPUT  PORT 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

fcebi  Center  Frequency 

inttime  Integration  time 

sfreq  Sampling  Frequency 

LIBRARY  /  FUNCTION 
utlib/bp_integrator 


C-17 


BINARY  DS  SPREADER 


DESCRIPTION 

Multiplies  the  complex  input  by  a  binary  complex  pseudonoise  sequence  with  a  chip  rate 
specified  by  the  user  as  a  multiple  of  the  data  rate.  If  the  input  signal  has  constant  unit 
magnitude,  the  output  power  is  1. 

MENU  PATH 

Spread  Spectrum  >  Direct  Sequence  >  Binary 
INPUT  PORT 

IN  Complex 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

nch  Number  of  chips  per  symbol 

fda  Data  rate 

sfreq  Sampling  frequency 

nsamples  Number  of  samples  in  simulation 

LIBRARY  /  FUNCTION 
walib/binary_ds 


BIT  ERROR  ESTIMATOR 


DESCRIPTION 

Compares  the  input  bit  stream  to  a  synchronized  reference  data  stream  on  the  global 
connector  called  ‘data’,  and  computes  the  fraction  of  the  total  number  of  input  bits  which  do 
not  match  the  reference.  The  module  writes  the  result  to  a  file  with  a  name  selected  by  the 
user. 

MENU  PATH 

Communicator  >  Metrics  >  Bit  Error  Rate 
INPUT  PORT 

BITS  OUT 
DELAYJN 
OUTPUT  PORT 
None 

PARAMETERS 
filename 
fda 
sfreq 

LIBRARY  /  FUNCTION 
vuIib/bit_error_rate 


Real 

Real 


File  name 
Bit  rate 
Sampling  rate 


C-19 


BPSK  DEMODULATOR 


DESCRIPTION 

Demodulates  a  binary  PSK  waveform,  producing  a  binary  sequence. 
MENU  PATH 

Communicator  >  Demodulator  >  BPSK 
INPUT  PORT 

IN  Complex 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 

fda  Bit  rate 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
modiib/bpsk_demod 


C-20 


BPSK  MODULATOR 


DESCRIPTION 

A  binary  PSK  modulator. 

MENU  PATH 

Communicator  >  Modulator  >  BPSK 
INPUT  PORT 

IN  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

fda  Bit  rate 

sfreq  Sampling  rate 

LIBRARY  /  FUNCTION 
modlib/bpsk_mod 


C-21 


BPSK  SOURCE 


DESCRIPTION 

Outputs  a  narrowband  complex  envelope  resulting  from  interleaving  a  pseudorandom  bit 
source,  and  using  the  result  as  the  data  stream  for  a  binary  phase  shift  keying  modulation. 
The  user  sets  the  numbers  of  columns  and  rows  to  define  the  interleaving,  and  accepting  the 
default  value  of  1  for  both  produces  a  non-interleaved  data  stream. 

MENU  PATH 

Communicator  >  Modulated  Source  >  BPSK 


INPUT  PORT 
None 

OUTPUT  PORT 
OUT 

DELAY_OUT 

PARAMETERS 

ncols 

nrows 

fda 

sfreq 

USAGE  EXAMPLES 


Complex 

Real 

Interleaver  columns 
Interleaver  rows 
Bit  rate 
Sampling  rate 


SSVM>Demonstrations>  Correlation  Radiometer 
SSVM>Demonstrations>  Scanner 


LIBRARY  /  FUNCTION 


modulib/bpsk_source 


BURST 


DESCRIPTION 

Reduces  the  duration  of  the  transmitted  signal  from  the  collect  time  to  the  burst  duration 

MENU  PATH 


Spread  Spectrum  >  Burst 
INPUT  PORT 

IN  Complex 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 


duration 

collect_time 

sfreq 


Scale  parameter 
Collect  time 
Sampling  frequency 


LIBRARY  /  FUNCTION 
walib/burst 


CAUCHY  NOISE 


DESCRIPTION 

Produces  radially  symmetric  complex  noise  whose  real  and  imaginary  components  are 
marginally  distributed  as  Cauchy  random  variables.  The  user  may  set  a  scaling  parameter, 
which  multiplies  a  standard  unit  variable.  The  Cauchy  variable  has  an  undefined  mean  and 
infinite  average  power,  so  that  signal  to  noise  ratio  has  no  meaning.  If  the  scaling  parameter 
is  a,  the  marginal  probability  densities  of  the  real  and  imaginary  parts  are  given  by 

1  1 
na  l  +  (jc/a)^ 


MENU  PATH 

Background  >  Noise  >  Cauchy 
INPUT  PORT 
None 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

mult  Scale  parameter 

LIBRARY  /  FUNCTION 
bklib/cauchy_noise 


C-24 


CENTER  FREQUENCY  MISMATCH 


DESCRIPTION 

Accounts  for  the  difference,  if  any,  between  the  center  frequency  of  the  communication 
signal  and  the  center  frequency  of  the  receiver. 

MENU  PATH 

Utilities  >  Signal  Processing  >  Center  Freq  Mismatch 
INPUT  PORT 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

fees  Signal  center  frequency 

f  ce  d  Receiver  center  frequency 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/cent_freq_mm 


C-25 


CHIP  RATE  AGILITY 


DESCRIPTION 

Outputs  trigger  pulses  at  a  rate  which  varies  pseudorandomly  from  one  block  of  time  to  the 
next,  but  is  constant  within  each  block. 

MENU  PATH 

Spread  Spectrum  >  Chip  Rate  Agility 
INPUT  PORT 
None 

OUTPUT  PORT 


OUT 

Complex 

PARAMETERS 

min_duration 

Minimum  block  length 

max_duration 

Maximum  block  length 

nchav 

Average  number  of  chips  per  data  symbol 

hw 

One  sided  percent  chip  length  fluctuation 

sfreq 

Sampling  frequency 

fda 

Data  rate 

Store_agility 

1  to  store  trigger  pulses  and  block  lengths,  0  otherwise 

LIBRARY  /  FUNCTION 


walib/chip_rate_agty 


C-26 


CHIP  RATE  AGILE  8ary  DS  SPREADER 


DESCRIPTION 

Multiplies  the  complex  input  by  an  8aiy  complex  pseudonoise  sequence,  which  has  a  rate 
that  remains  constant  over  blocks  of  time  and  varies  pseudorandomly  between  blocks.  The 
chip  rate  agile  DS  spreaders  use  the  Chip  Rate  Agility  module  to  trigger  the  generation  of  the 
M-ary  symbols  used  in  spreading.  The  chip  rates  result  from  sampling  a  uniform 
distribution  and  rounding  to  the  nearest  multiple  of  the  data  rate.  The  Chip  Rate  Agility 
module  uses  the  Agility  module  to  trigger  generation  of  new  blocks.  Block  lengths  vary 
pseudorandomly,  and  are  rounded  to  the  nearest  multiple  of  the  data  symbol  duration. 

The  user  specifies  the  minimum  and  maximum  block  lengths;  the  average  chip  rate  as  a 
multiple  of  the  data  rate;  the  maximum  percentage  one  sided  chip  rate  fluctuation;  and 
whether  the  block  lengths,  block  boundary  pulses,  instantaneous  chip  rate,  and  chip 
boundary  pulses  are  stored.  If  the  input  signal  has  constant  unit  magnitude,  the  output 
power  is  1. 

Delays  the  output  signal  by  the  input  delay  value  for  automatic  synchronization.  Outputs  the 
spreading  sequence  as  a  global  output  connector  called  ‘code’. 

MENU  PATH 


Spread  Spectrum  >  Agile  Direct  Sequence  >8ary 


INPUT  PORT 

IN 

Complex 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

min_duration 

Minimum  block  length 

max_du  ration 

Maximum  block  length 

nchav 

Average  number  of  chips  per  data  symbol 

hw 

One  sided  percent  chip  length  fluctuation 

store_agiIity 

1  to  store  block  lengths,  block  boundary  pulses,  chip  rate  and  chip 
boundary  pulses;  0  otherwise 

sfreq 

Sampling  frequency 

fda 

Data  rate 

LIBRARY  /  FUNCTION 
walib/a_8ary_ds 


C-27 


CHIP  RATE  AGILE  BINARY  DS  SPREADER 


DESCRIPTION 

Multiplies  the  complex  input  by  a  binary  complex  pseudonoise  sequence,  which  has  a  rate 
that  remains  constant  over  blocks  of  time  and  varies  pseudorandomly  between  blocks.  The 
chip  rate  agile  DS  spreaders  use  the  CHIP  RATE  AGILITY  module  to  trigger  the  generation 
of  the  M-ary  symbols  used  in  spreading.  The  chip  rates  result  from  sampling  a  uniform 
distribution  and  rounding  to  the  nearest  multiple  of  the  data  rate.  The  CHIP  RATE 
AGILITY  module  uses  the  AGILITY  module  to  trigger  generation  of  new  blocks.  Block 
lengths  vary  pseudorandomly,  and  are  rounded  to  die  nearest  multiple  of  the  data  symbol 
duration. 

The  user  specifies  the  minimum  and  maximum  block  lengths;  the  average  chip  rate  as  a 
multiple  of  the  data  rate;  the  maximum  percentage  one  sided  chip  rate  fluctuation;  and 
whether  the  block  lengths,  block  boundary  pulses,  instantaneous  chip  rate,  and  chip 
boundary  pulses  are  stored.  If  the  input  signal  has  constant  unit  magnitude,  the  ouqiut 
power  is  1. 

Delays  the  output  signal  by  the  input  delay  value  for  automatic  synchronization.  Outputs  the 
spreading  sequence  as  a  global  output  connector  called  ‘code’. 

MENU  PATH 

Spread  Spectrum  >  Agile  Direct  Sequence  >  Binary 
INPUT  PORT 


IN 

Complex 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

min_duration 

Minimum  block  length 

max_duration 

Maximum  block  length 

nchav 

Average  number  of  chips  per  data  symbol 

hw 

One  sided  percent  chip  length  fluctuation 

store_agility 

1  to  store  block  lengths,  block  boundary  pulses,  chip  rate  and  chip 
boundary  pulses;  0  otherwise 

sfreq 

Sampling  frequency 

fda 

Data  rate 

LIBRARY  /  FUNCTION 


walib/a_binary_ds 


C-28 


CHIP  RATE  AGILE  QUATERNARY  DS  SPREADER 


DESCRIPTION 

Multiplies  the  complex  input  by  a  binary  complex  pseudonoise  sequence,  which  has  a  rate 
that  remains  constant  over  blocks  of  time  and  varies  pseudorandomly  between  blocks.  The 
chip  rate  agile  DS  spreaders  use  the  CHIP  RATE  AGILITY  module  to  trigger  the  generation 
of  the  M-ary  symbols  used  in  spreading.  The  chip  rates  result  from  sampling  a  uniform 
distribution  and  rounding  to  the  nearest  multiple  of  the  data  rate.  The  CHIP  RATE 
AGILITY  module  uses  the  AGILITY  module  to  trigger  generation  of  new  blocks.  Block 
lengths  vary  pseudorandomly,  and  are  rounded  to  Ae  nearest  multiple  of  the  data  symbol 
duration. 

The  user  specifies  the  minimum  and  maximum  block  lengths;  the  average  chip  rate  as  a 
multiple  of  the  data  rate;  the  maximum  percentage  one  sided  chip  rate  fluctuation;  and 
whether  the  block  lengths,  block  boundary  pulses,  instantaneous  chip  rate,  and  chip 
bounda^  pulses  are  stored.  If  the  input  signal  has  constant  unit  magnitude,  the  ouq)ut 
power  is  1. 

Delays  the  output  signal  by  the  input  delay  value  for  automatic  synchronization.  Outputs  the 
spreading  sequence  as  a  global  output  connector  called  ‘code’. 

MENU  PATH 

Spread  Spectrum  >  Agile  Direct  Sequence  >  Quaternary 
INPUT  PORT 


IN 

Complex 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

min_duration 

Minimum  block  length 

max_duration 

Maximum  block  length 

nchav 

Average  number  of  chips  per  data  symbol 

hw 

One  sided  percent  chip  length  fluctuation 

store_agility 

1  to  store  block  lengths,  block  boundary  pulses,  chip  rate  and  chip 
boundary  pulses;  0  otherwise 

sfreq 

Sampling  frequency 

fda 

Data  rate 

LIBRARY  /  FUNCTION 


walib/a_quatary_ds 


C-29 


CHIP  RELEASE 


DESCRIPTION 

Outputs  trigger  pulses  at  a  rate  equal  to  the  input  signal.  The  output  pulses  result  from 
differentiation  a  square  wave,  which  in  turn  is  formed  by  thresholding  a  siiiusoid  with 
frequency  equal  to  the  input  signal.  The  frequency  is  not  rounded  to  synchronize  with  the 
interval  tetween  simulation  iterations. 

MENU  PATH 

Utilities  >  Sources  >  Chip  Release 
INPUT  PORT 
FREQ 

OUTPUT  PORT 
PULSES 
PARAMETERS 
sfreq 

LIBRARY  /  FUNCTION 
walib/chip_release 


Real 

Real 

Sampling  frequency 


C-30 


COMMUNICATION  COUPLER 


DESCRIPTION 

Adds  the  transmitted  signal  and  background  noise  to  form  the  signal  received  by  the 
communication  receiver. 

MENU  PATH 

Channel  >  Couplers  >  Communicator 
INPUT 


SIGNAL 

Real 

BACKGROUND 

Real 

OUTPUT 

SIG+NOISE 

Real 

PARAMETERS 

None 

LIBRARY  /  FUNCTION 
chlib/comm_coupler 


C-31 


COMMUNICATORS  COMPOSITE  PREPROCESSOR 


DESCRIPTION 


Lowpass  filters  the  input  signal  then  optionally  applies  a  frequency  based  excision 
technique. 

MENU  PATH 

Communicator 

>  Preprocessor  >  Composite 

INPUT  PORTS 

IN 

Complex 

DELAY_IN 

Real 

OUTPUT  PORTS 

OUT 

Complex 

DELAY_OUT 

Real 

PARAMETERS 

Lc 

Filter  cutoff  frequency 

delta_f 

Filter  transition  region  width 

a_r 

Filter  stopband  attenuation 

sfreq 

Sampling  frequency 

index 

Selection  index  (0-no  excision,  1 -fixed,  2- Adaptive,  3-Whitening) 

store_spectra 

Store  spectra?  (index  >  0) 

blocksize 

Number  of  points  used  in  the  FFT’s  (index  >  0) 

nlargest 

Number  of  bins  to  be  excised  (index  =1) 

exceedfactor 

Exceedance  factor  (index  =  2) 

nexempt 

Number  of  points  excluded  in  the  calculation  of  the  mean  square 
magnitude  (index  =  2) 

LIBRARY  /  FUNCTION 


prlib/communication 


C-32 


COMMUNICATOR’S  ADAPTIVE  THRESHOLD  EXCISION  SYSTEM 


DESCRIPTION 

Performs  a  spectral  analysis  of  the  input  and  removes  the  components  with  squared 
magnitnHft  greater  than  an  exceedance  factor  times  a  mean  square  magnitude.  The  mean 
square  magnitude  is  of  all  the  bins  except  a  number  of  strongest  bins.  TTie  user  specifies  the 
size  of  the  FFT’s,  the  number  of  strongest  bins  to  be  excluded  from  the  calculation  of  mean 
square  magnitude,  the  exceedance  factor,  and  the  type  of  window  used  in  the  forward  FFT. 
The  inverse  FFT  always  uses  a  rectangular  window.  This  algorithm  recomputes  the 
threshold  for  each  successive  block  of  data  and  so  adapts  to  changing  background  levels. 

The  user  has  the  option  of  storing  the  spectra  before  and  after  excision.  The  user  specifies  a 
filename  stem,  and  the  Vulnerability  Metric  appends  “.in”  and  “.out”  to  stem  to  form  the 
names  of  the  files  storing  the  spectra  before  and  after  excision  respectively.  The 
Vulnerability  Metric  writes  all  signd  files  to  the  directory  /spwdataA^nsigs. 

Taking  the  FFT  delays  die  signal  by  a  number  of  iterations  equal  to  the  FFT  size,  so  the 
signal  value  presented  at  the  DELAY  OUT  port  is  equal  to  the  value  at  the  DELAY  IN  port 
plus  the  FFT  size. 

MENU  PATH 

Communicator  >  Preprocessor  >  Adaptive  Threshold  Excision 
INPUT  PORTS 
IN 

HOLD 

OUTPUT  PORTS 
OUT 

PARAMETERS 

exceedfactor 
nexempt 


blocksize 

window_type 


Complex 

Real 

Complex 

Exceedance  factor 

Number  of  points  excluded  in  the  calculation  of  the  mean  square 
magnitude 

Number  of  points  used  in  the  FFTs 
Type  of  window  used  in  the  forward  FFT 
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COMMUNICATOR’S  ADAPTIVE  THRESHOLD  EXCISION  SYSTEM,  continued 


store_spectra  Store  spectra? 

s  f r  e  q  Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 
prlib/exadapt 


C-34 


COMMUNICATOR’S  FIXED  EXCISION  SYSTEM 


DESCRIPTION 

Performs  a  spectral  analysis  of  the  input  and  removes  the  strongest  components.  The  user 
specifies  the  size  of  the  FFT’s,  the  number  of  strongest  bins  to  be  excised,  and  the  ^  of 
window  used  in  the  forward  FFT.  The  inverse  ^  always  us^  a  rectanpl^mdow 
The  frequency  resolution,  or  bin  size,  is  the  sampling  frequency  divided  by  me  FFT  siz^  11 
the  FFT  size  is  N,  simulation  samples  numbered  0  through  N-2  are  in  the  first  block,  N-1 
through  N-2  in  the  second  block,  and  so  on. 

The  user  has  the  option  of  storing  the  spectra  before  and  after  excision.  The  user  specifies  a 
filename  stem,  and  the  Vulnerability  Metric  appends  “.in”  and  “.out”  to  stem  to  form  fte 
names  of  the  files  storing  the  spectra  before  and  after  excision  respectively.  The 
Vulnerability  Metric  writes  all  signi  files  to  the  directory  /spwdata/vmsigs. 

Taking  the  FFT  delays  the  signal  by  a  number  of  iterations  equal  to  the  FFT  size,  so  tiie 
signal  value  presented  at  the  DELAY  OUT  port  is  equal  to  the  value  at  the  DELAY  IN  port 
plus  the  FFT  size. 

MENU  PATH 

Communicator  >  Preprocessor  >  Fixed  Excision 


INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORTS 


OUT  Complex 

PARAMETERS 


nlargest 

blocksize 

window_type 

store_spectra 

sfreq 

LIBRARY  /  FUNCTION 


Number  of  bins  to  be  excised 
Number  of  points  used  in  the  FFTs 
Type  of  window  used  in  the  forward  FFT 
Store  spectra? 

Simulation  sampling  frequency 


prlib/enfixed 


C-35 


COMMUNICATOR’S  FREQUENCY  BASED  EXCISION 


DESCRIPTION 

Optionally  applies  one  of  three  frequency  based  excision  techniques:  fixed,  adaptive 
threshold,  or  whitening. 

MENU  PATH 

Communicator  >  Preprocessor  >  Composite  Freq  Based  Excision 
INPUT  PORTS 


IN 

Complex 

DELAYJN 

Real 

OUTPUT  PORTS 

OUT 

Complex 

DELAY_OUT 

Real 

PARAMETERS 

index 

Selection  index  (0-no  excision,  1 -Fixed,  2- Adaptive,  3-Whitening) 

store_spectra 

Store  spectra?  (index>  0) 

blocksize 

Number  of  points  used  in  the  FFTs 

window_type 

Type  of  window  used  in  the  forward  FFT 

nlargest 

Number  of  bins  to  be  excised 

exceedfactor 

Exceedance  factor  )index  =1) 

nexempt 

Number  of  points  excluded  in  the  calculation  of  mean  square 
magnitude  (index  =  2) 

sfreq 

Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 
prlib/freq_excision 
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COMMUNICATOR’S  WHITENING  EXCISION  SYSTEM 


DESCRIPTION 

Performs  a  spectral  analysis  of  the  input  and  divides  each  spectral  component  by  its 
magnitude  resulting  in  flat  output  spectrum..  The  user  specifies  the  size  of  the  FFT’s,  and 
the  type  of  window  used  in  the  forward  FFT.  The  inverse  FFT  always  uses  a  rectangular 
window.  The  spectral  resolution,  or  bin  size,  is  the  sampling  frequency  divided  by  the  FFT 
size.  If  the  input  spectrum  in  any  bin  is  smaller  than  the  machine’s  arithmetic  resolution,  tfie 
magnitude  is  set  to  1,  and  the  phase  is  arbitrarily  set  to  0.  If  the  FFT  size  is  N,  simulation 
samples  numbered  0  through  N-2  are  in  the  first  block,  N-1  through  N-2  in  the  second 
block,  and  so  on. 

The  user  has  the  option  of  storing  the  spectra  before  and  after  excision.  The  user  specifies  a 
filename  stem,  and  the  Vulnerability  Metric  appends  “.in”  and  “.out”  to  stem  to  form  the 
names  of  the  files  storing  the  spectra  before  and  after  excision  respectively.  The 
Vulnerability  Metric  writes  all  signd  files  to  the  directory  /spwdataA^^isigs. 

Taking  the  FFT  delays  the  signal  by  a  number  of  iterations  equal  to  the  FFT  size,  so  the 
signal  value  presented  at  the  DELAY  OUT  port  is  equal  to  the  value  at  the  DELAY  IN  port 
plus  the  FFT  size. 

MENU  PATH 


Communicator  >  Preprocessor  >  Whitening  Excision 


INPUT  PORTS 
IN 

DELAY_IN 

HOLD 

DELAY_OUT 
OUTPUT  PORTS 
OUT 

PARAMETERS 

blocksize 

window_type 


Complex 

Real 

Real 

Real 

Complex 

Number  of  points  used  in  the  FFTs 
Type  of  window  used  in  the  forward  FFT 
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COMMUNICATOR’S  WHITENING  EXCISION  SYSTEM,  continued 


store_spectra  Store  spectra? 

sfreq  Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 
prlib/exwhite 
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COMPLEX  INPUT  SWITCH 


DESCRIPTION 

Selects  the  output  from  one  of  the  complex  inputs.  Parameters  specify  both  the  number  of 
permitted  inputs,  and  which  one  is  used  for  die  output.  The  naodule  tests  the  value  of  the 
both  parameters,  and  terminates  with  an  error  message  if  either  is  out  of  range.  Setting  the 
parent  m<xlu1e  parameter  permits  the  error  message  to  give  the  containing  module  of  the 
instance  that  has  an  error. 


MENU  PATH 

Utilities  >  Logic  >  Switch,  Complex,  Parameter 
INPUT  PORTS 


IN_0 

Complex 

IN_1 

Complex 

IN_2 

Complex 

IN_3 

Complex 

IN_4 

Complex 

IN_5 

Complex 

OUTPUT  PORTS 

OUT 

Complex 

PARAMETERS 

no_of_inputs 

Number  of  inputs  (1-6) 

index 

Selection  index  (0-no_of_inputs-l) 

LIBRARY  /  FUNCTION 
utlib/switch_in 
Coded  in  C 
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COMPLEX  MOVING  AVERAGE 


DESCRIPTION 

Outputs  the  input  averaged  over  an  immediately  preceding  time  interval. 

MENU  PATH 

Utilities  >  Estimators  >  Complex  Moving  Avge 
INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

sfreq  Simulation  sampling  frequency 

inttime  Length  if  signal  averaging  time  interval 

LIBRARY  /  FUNCTION 

utlib/cmplx_mv_avge 
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COMPLEX  SWITCHED  INPUT 


DESCRIPTION 

Outputs  the  top  signal  input  if  the  control  parameter  is  0  and  outputs  the  bottom  signal 
input  if  the  control  parameter  is  1. 

MENU  PATH 

Utilities  >  Logic  >  Switch,  Complex,  Signal 
INPUT  PORTS 


IN 

Complex 

CONTROL 

Integer 

HOLD 

Real 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

None 

LIBRARY  /  FUNCTION 
utlib/cmplx_sw_in 
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COMPLEX  VARIANCE 


DESCRIPTION 

Outputs  the  mean  and  variance  of  the  input  signal,  computed  from  the  beginning  of  the 
simulation. 

MENU  PATH 

Utilities  >  Estimators  >  Complex  Variance 


INPUT  PORTS 

IN 

Complex 

HOLD 

Real 

OUTPUT  PORTS 

MEAN 

Complex 

VAR 

Real 

PARAMETERS 

None 

LIBRARY  /  FUNCTION 
utlib/cmplx_var 
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CORRELATION  DOUBLER  SYSTEM 


DESCRIPTION 

Uses  a  communication  signal  input  and  a  pair  of  background  inputs  to  form  a  signal  present 
pair  and  signal  absent  pair,  and  passes  the  signal  present  pair  through  one  correlation 
doubler  detector  and  the  signal  absent  pair  through  a  second  correlation  doubler  detector. 
The  result  of  processing  the  signal  present  is  the  Hj  output  and  the  result  of  processing  the 

signal  absent  pair  is  the  Hq  output  Stores  detector  parameter  values  in  files. 

MENU  PATH 

Interceptor  >  Detector  >  Freq  Multr  >  Corr  Doubler,  system 
INPUT  PORTS 


SIGNAL 

Complex 

BACKGROUND!  Complex 

BACKGROUND! 

Complex 

OUTPUT  PORTS 

HI 

Complex 

HO 

Complex 

PARAMETERS 

feed 

Receiver  center  frequency 

bwif 

Receiver  bandwidth 

collect_time 

Output  integration  time  in  detector  systems 

wswgn 

1  for  weak  signal,  white  Gaussian  noise  mode,  0  otherwise 

nwssb 

1  for  nonweak  signal,  simulated  background  mode,  0  otherwise 

fees 

Signal  center  frequency 

sfreq 

Simulation  sampling  frequency 

nsamples 

Number  of  samples  in  simulation 

fda 

Data  rate 

LIBRARY  /  FUNCTION 


delib/cor_doub 
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CORRELATION  DOUBLER  CORE 

DESCRIPTION 

Simulates  the  action  of  a  correlation  doubler  on  complex  envelope  inputs.  Provides  optional 
storage  of  the  outputs  of  one  receiver,  the  multiplier  and  the  detector  system. 

MENU  PATH 

Interceptor  >  Detector  >  Freq  Multr  >  Corr  Doubler,  core 
INPUT  PORTS 

INI  Complex 

IN2  Complex 

HOLD  Complex 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

write_sigs  1  to  store  intermediate  signals,  0  otherwise 

feed  Receiver  center  frequency 

b  w  i  f  Receiver  bandwidth 

int_time  Integration  time  in  output  integrator 

fees  Signal  center  frequency 

sfreq  Sampling  frequency 

nsamples  Numbers  of  samples  in  simulation 

LIBRARY/FUNCTION 

delib/cor  doub  core 
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CORRELATION  RADIOMETER  SYSTEM 


DESCRIPTION 

Uses  a  coi^unication  signal  input  and  a  pair  of  background  inputs  to  form  a  signal  present 
pair  and  signal  absent  pair,  and  passes  the  signal  present  pair  through  one  correlation 
radiometer  detector  and  the  signd  absent  pair  through  a  second  correlation  radiometer 
detector.  The  result  of  processing  the  signal  present  pair  is  the  Hj  output  and  the  result  of 

processing  the  signal  absent  pair  is  the  Hq  output.  Stores  detector  parameter  values  in  files. 
MENU  PATH 


Interceptor  >  Detector  >  Radiometer  >  Correlation,  system 
INPUT  PORTS 


SIGNAL 

Complex 

BACKGROUND!  Complex 

BACKGROUND2 

Complex 

OUTPUT  PORTS 

HI 

Complex 

H2 

Complex 

PARAMETERS 

feed 

Receiver  center  frequency 

bwif 

Receiver  bandwidth 

coIIect_time 

Output  integration  time  in  detector  frequency 

wswgn 

1  for  weak  signal,  white  Gaussian  noise  mode,  0  otherwise 

nwssb 

1  for  non  weak  signal,  simulated  background  mode,  o  otherwise 

fees 

Signal  center  frequency 

sfreq 

Sampling  frequency 

C-45 


CORRELATION  RADIOMETER  SYSTEM,  continued 


nsamples  Number  of  samples  in  simulation 

fda  Data  rate 

USAGE  EXAMPLE 

SSVM>Demonstrates>CorreIation  Radiometer 
LIBRARY  /  FUNCTION 
delib/cor_radiom 
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CORRELATION  RADIOMETER  CORE 

DESCRIPTION 

Simulates  the  action  of  a  correlation  doubler  on  complex  envelope  inputs.  Provides  optional 
storage  of  the  outputs  of  one  receiver,  the  multiplier  and  the  detector  system. 

MENU  PATH 

Interceptor  >  Detector  >  Freq  Multr  >  Corr  Doubler,  core 
INPUT  PORTS 


INI 

Complex 

IN2 

Complex 

HOLD 

Complex 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

write_sigs 

1  to  store  intermediate  signals,  0  otherwise 

feed 

Receiver  center  frequency 

bwif 

Receiver  bandwidth 

int_time 

Integration  time  in  output  integrator 

wswgn 

1  for  weak  signal,  white  Gaussian  noise  mode,  0  otherwise 

nwssb 

1  for  nonweak  signal,  simulated  background  mode,  0  otherwise 

fees 

Signal  center  frequency 

sfreq 

Sampling  frequency 

nsamples 

Number  of  samples  in  simulation 

USAGE  EXAMPLE 

SSyM>Demonstrations>CorreIation  Radiometer 
Double  click  on  “CORRELATION  RADIOMETER”  module 


LIBRARY  /  FUNCTION 

delib/cor  radiom  core 


C-47 


CUSTOM  INTERFERENCE 


DESCRIPTION 

Reads  an  interference  configuration  supplied  by  the  user  and  adds  the  result  to  the  signd. 
The  user  specifies  the  number  of  interferers,  the  configuration  center,  interferer  powers,  in 
dB  above  communication  signal  power,  and  the  widths.  The  bandwidth  is  implemented  by 
BPSK  modulating  each  int^erence  carrier  frequency  equal  to  the  bandwidth.  The  user 
specifies  the  name  of  the  file  configuration,  the  default  name  begins  custom.conf.  The 
format  of  the  file  is  a  format  created  by  the  Random  Interference  Module,  so  that 
configuration  randomly  and  saved  may  be  reused  directly  by  the  Custom  Interference 
Module. 

The  mathematical  form  of  the  interference  is  the  same  as  that  detailed  for  Random 
Interference. 

MENU  PATH 

Background  >  Interference  >  Custom 
INPUT  PORT 
None 

OUTPUT  PORT 

OUT  Complex 

INPUT  FILE 

data/custom.conf 

The  input  file  has  the  following  format: 

Interference  Configuration 

Num_Freqs:  100  Center_Freq:  l.OOOOOOe+06 
-3.363  lOOe+05  +2.2735  lOe+05  ...  +3.575460e+05 

Powers: 

9.535084e+00  9.055280e+00  ...  1.087717e+01 
Widths: 

3.341063e+02  4.591780e+02  ...  8.631209e+02 
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CUSTOM  INTERFERENCE,  continued 


PARAMETERS 

conf_file 

Configuration  input  filename,  defaults  to  custom.conf 

simtime 

Simulation  time 

sfreq 

Sampling  frequency 

fees 

Signal  center  frequency 

LIBRARY  /  FUNCTION 
bklib/ifcus 
Codedin  C 


C-49 


DATA  SOURCE 


DESCRIPTION 

Outputs  a  pseudorandom  binary  sequence. 
MENU  PATH 

Utilities  >  Input  /Output  >  Data  Source 
INPUT  PORTS 
None 

OUTPUT  PORT 

OUT  Complex 

DELAY_OUT  Real 
PARAMETERS 

fda  Data  rate 

sf  r e  q  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/data_source 
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DEINTERLEAVER 


DESCRIPTION 

Implements  a  block  deinterleaver,  which  reverses  the  effect  of  the  block  interleaver.  An 
NxM  deinterleaver  is  the  same  as  an  MxN  interleaver. 

MENU  PATH 

Communicator  >  Demodulator  >  Deinterleaver 
INPUT  PORTS 

INI  Real 

HOLD  Real 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 

n  c  o  1  s  Number  of  columns  in  corresponding  interleaver  matrix 

n  r  o  ws  Number  of  rows  in  corresponding  interleaver  matrix 

fda  Data  rate 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
modlib/deinterleave 
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DELAY  CHIP  RATE  DETECTOR  SYSTEM 

DESCRIPTION 

Transmits  the  signal  present  and  signal  absent  inputs  to  separate  delay  chip  rate  detectors. 
Outputs  the  results,  as  well  as  the  noise  bandwidth. 

MENU  PATH 

Interceptor  >  Detector  >  Chip  Rate  >  Delay,  system 
INPUT  PORTS 


SIGNAL 

Complex 

BACKGROUND 

Complex 

OUTPUT  PORT 

HI 

Complex 

HO 

Complex 

NOISE_BW 

Real 

PARAMETERS 

feed 

Receiver  center  frequency 

bwif 

Receiver  bandwidth 

delay_factor 

Delay,  in  units  of  half  the  estimated  chip  dilation 

wswgn 

1  for  weak  signal,  white  Gaussian  noise  mode,  0  otherwise 

nwssb 

1  for  nonweak  signal,  simulated  background  mode,  0  otherwise 

fees 

Signal  center  frequency 

sfreq 

Sampling  frequency 

nsamples 

Number  of  samples  in  simulation 

fda 

Data  rate 
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DELAY  CHIP  RATE  DETECTOR  SYSTEM,  continued 


collect_time  Output  integration  time  in  detector  systems 

LIBRARY  /  FUNCTION 
delib/delay_chip 
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DELAY  CHIP  RATE  DETECTOR  CORE 


DESCRIPTION 

Simulates  the  action  of  a  delay  chip  rate  detector  with  complex  envelope  input  Provides 
optional  storage  of  intermediate  signals. 

MENU  PATH 

Interceptor  >  Detector  >  Chip  Rate  >  Delay,  core 
INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

write_sigs  1  to  store  intermediate  signals,  0  otherwise 

feed  Receiver  center  frequency 

b  w  i  f  Receiver  bandwidth 

coIIect_time  Output  integration  time 

fees  Signal  center  frequency 

sfreq  Sampling  frequency 

nsamples  Number  of  samples  in  simulation 

LIBRARY  /  FUNCTION 

deHb/delay_chip_cor 


C-54 


DELAY,  FFT  SEARCH,  CHIP  RATE  DETECTOR 


DESCRIPTION 

Implements  a  delay  chip  rate  detector  with  an  FFT  output  stage  to  permit  detection  at 
midtiple  chip  rates.  In  addition  the  output  stage  provides  filtering,  subsampling  and  block 
averaging.  The  user  can  therefore  tailor  the  search  to  cover  efficiently  the  range  of  chip  rate 
uncertainty.  The  documentation  of  the  Envelope  Spectram  Detector  provides  detailed 
information  on  the  filtering,  subsampling  and  block  averaging. 

MENU  PATH 


Interceptor  >  Detector  >  Chip  Rate  >  Delay,  FFT  Search 
INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 
None 


PARAMETERS 

feed 

bwif 

filter_^order 

edge_atten 

store_signals 

sfreq 

fees 

fchd 

delay_factor 
start  time 


Receiver  center  frequency 
Receiver  bandwidth 
Receiver  filter  order 

Receiver  passband  edge  attenuation,  dB  power 

1  to  store  intermediate  signals,  0  otherwise 

Sampling  frequency 

Signal  center  frequency 

Interceptor  chip  rate  estimate 

Delay 

Time  to  begin  collecting  date  to  determine  averaged  spectrum 
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DELAY,  FFT  SEARCH,  CHIP  RATE  DETECTOR,  continued 


spec_bwif 

subsampling 

fft_size 

window_type 

spec_avg_biock 

freq_zoom 


Filter  width  (lowpass,  2-sided) 

Subsampling  order 
FFT  size,  power  of  2 

Window  (blackman,  bartlett,  hamming,  banning,  rectangular) 
Number  of  spectra  to  average 
Frequency  zoom,  power  of  2 


USAGE  EXAMPLE 


SSVM>Demonstrations>Search  Chip  Rate  Detector 
LIBRARY  /  FUNCTION 
delib/delay_chip_fft 
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DELAY,  STEPPING  SEARCH,  CHIP  RATE  DETECTOR 


DESCRIPTION 

Combines  a  delay  chip  rate  detector  with  a  two-dimensional  search  that  scans  receiver  center 
frequency  and  chip  rate.  The  chip  rate  cycles  throu^  a  fixed  number  of  values  specified  by 
the  user,  while  the  receiver  center  frequency  keeps  incrementing  throughout  the  simulation, 
once  after  each  cycle  through  the  chip  rate,  so  the  number  of  center  frequencies  used 
depends  on  the  number  of  collects  simulated  and  the  number  of  chip  rates  used.  The  user 
controls  the  order  and  passband  edge  attenuation  of  the  receiver  filter  and  specifies  the 
name,  filename,  of  the  file  to  which  the  results  are  written. 

After  the  simulation,  the  user  obtains  a  Gnuplot  screen  plot  by  entering 

DelSrch  filename 

in  an  operating  system  window.  The  plot  also  displays  the  name  of  the  underlying  data  file. 
To  obtain  a  hard  copy  in  addition  to  the  screen  plot,  the  user  enters 

DelSrch  -h  filename 

If  the  default  printer  has  not  been  set  during  the  session,  the  user  may  have  to  specify  the 
printer  to  be  used  by  entering,  in  the  operating  system  window,  setenv  PRINTER 
followed  by  the  printer  name.  Quitting  the  graphics  window,  or  typing  a  keystroke  in  the 
operating  system  window  causes  the  graphical  output  to  disappear,  and  die  Gnuplot  session 
to  close. 

DelSrch  is  a  Bourne  shell  script  which  selects  and  edits  the  file 

DelSrch  .hard,  template  if  the  flag  -h  is  present,  Md 

DelSrch .  soft .  template  otherwise,  and  calls  Gnuplot  to  execute  the  resulting 
Gnuplot  script,  DelSrch .  gnu .  The  data  file  and  the  script  files  are  in  the  directory 
/spwdata/vmsign,  but  DelSrch  can  be  invoked  fiom  any  directory  without 
specifying  a  path. 

MENU  PATH 

Interceptor  >  Detector  >  Chip  Rate  >  Delay,  Step  Search 
INPUT  PORTS 

IN  Complex 

HOLD  Complex 

OUTPUT  PORTS 
None 


C-57 


DELAY,  STEPPING  SEARCH,  CHIP  RATE  DETECTOR,  continued 


PARAMETERS 

start_fced 

Starting  receiver  center  frequency 

step_fced 

Receiver  center  frequency  increment 

bwif 

Receiver  bandwidth 

fchd 

Interceptor  chip  rate  estimate 

fiIter_order 

Receiver  filter  order 

edge_atten 

Receiver  passband  edge  attenuation,  dB  power 

store_signals 

1  to  store  intermediate  signals,  0  otherwise 

start_fchd 

Starting  chip  rate 

step_fchd 

Chip  rate  increment 

collect_tinie 

Output  integration  time 

sfreq 

Sampling  frequency 

fees 

Signal  center  frequency 

USAGE  EXAMPLE 

SSyM>Demonstrations>Search  Chip  Rate  Detector 


LIBRARY  /  FUNCTION 
delib/del_chip_srch 


C-58 


DENSITY  BASED  PREPROCESSOR 


DESCRIPTION 

Implements  the  Digital  Density  Detector.  Estimates  the  probability  density  of  the  additive 
noise  by  computing  a  histogram  over  a  fixed  data  block.  The  probability  density  estimate  is 
used  to  compute  a  nonlinearity  which  is  then  applied  to  the  input  signal,  to  produce  a 
nonlinearly  processed  data  signal,  delayed  by  the  length  of  the  block.  Tbe  vector 
corresponding  to  the  accumulated  histogram  vdues  is  also  output,  on  a  block-by-block 
basis. 

MENU  PATH 
None 


INPUT  PORT 
IN 

OUTPUT  PORTS 
SIG 
HIST 

PARAMETERS 

blk_size 

maxval 

nbins 


Complex 

Complex 
Real  Vector 

Processing  block  size 
Clipping  Level 

Number  of  bins  for  histogram  computation 


LIBRARY  /  FUNCTION 


prlib/ddd_sys 


C-59 


DETECTOR  INPUT  SWITCH 


DESCRIPTION 

Switches  the  signal  present,  Hj,  and  signal  absent,  Hq,  detector  inputs  to  the  two  detector 

cores,  as  required  by  the  selected  Operating  Mode.  The  Mode  selection  is  done  for  the  user 
based  on  the  specification  of  global  parameters  in  the  Vulnerability  Metric  Viewpoint 

MENU  PATH 

Interceptor  >  Detector  >  Utilities  >  Detector  Input  Switch 
INPUT  PORTS 

SIGNAL  Complex 

BACKGROUND  Real 
OUTPUT  PORTS 

H 1  Complex 

HO  Complex 

PARAMETERS 

long_stat_con  Long,  stationary  collect? 
quad_weak  Quadratic  processing,  with  weak  signal? 

Bideal_bkgd  Ideal  background? 

LIBRARY  /  FUNCTION 
delib/mode_switch 


C-60 


DIGITAL  TO  ANALOG  CONVERTER 


DESCRIPTION 

Convey  a  binary  vector  to  a  real  number.  The  low  order  bit  is  interpreted  as  the  one’s 
coefficient,  the  next  higher  bit  as  the  two’s  coefficient,  and  so  on.  The  output  is  a  real 
number  between  0  and  2  (number  of  bits  in)  -1. 

MENU  PATH 

Utilities  >  Signal  Processing  >  D/A  Converter 
INPUT  PORTS 


IN 

Real  Vector 

LOAD_IN 

Real 

HOLD 

Real 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

INJOVEC_LEN 

Number  of  bits  in  the  input  vector 

write_sigs 

1  to  store  intermediate  signals,  0  otherwise 

fda 

Data  rate 

s  f  r  e  q  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/d_a 

Coded  in  C 


C-61 


DIRECT  SEQUENCE  DESPREADER 


DESCRIPTION 

Reverses  the  effect  of  direct  sequence  spreading.  Since  the  spreading  sequence  is  conveyed 
by  means  of  a  global  connector,  a  single  block  serves  from  arbitrary  number  of  phases,  and 
for  any  type  of  jittering  of  agility. 

MENU  PATH 

Spread  Spectrum/DS  Despreading 
INPUT  PORT 

IN  Complex 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 
None 

LIBRARY  /  FUNCTION 
walib/ds_despreader 
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ENVELOPE  CHIP  RATE  DETECTOR  SYSTEM 


DESCRIPTION 

Uses  a  communication  signal  input  and  a  background  input  to  form  a  signal  present  input  to 
one  detector  and  a  signal  absent  input  to  another  detector.  The  detectors  are  each  envelope 
chip  rate  detectors.  The  result  of  processing  the  signal  present  case  is  the  Hj  output  and  Ae 

result  of  processing  the  signal  absent  case  is  the  Hq  output.  Stores  detector  parameter 
values  in  files. 

MENU  PATH 

Interceptor  >  Detector  >  Chip  Rate  >  Envelope,  system 
INPUT  PORTS 


SIGNAL 

Complex 

BACKGROUND 

Complex 

OUTPUT  PORT 

HI 

Complex 

H2 

Complex 

PARAMETERS 

feed 

Receiver  center  frequency 

bwif 

Receiver  bandwidth 

fchd 

Interceptor’s  estimate  of  chip  rate 

wswgn 

1  for  weak  signal,  white  Gaussian  noise  mode,  0  otherwise 

nwssb 

1  for  nonweak  signal,  simulated  background  mode,  0  otherwise 

fees 

Signal  center  frequency 

sfreq 

Sampling  frequency 

nsamples 

Number  of  samples  in  simulation 

fda 

data  rate 

conect_time 

Output  integration  time  in  detector  systems 

LIBRARY  /  FUNCTION 


delib/envel_chip 
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ENVELOPE  CHIP  RATE  DETECTOR  CORE 

DESCRIPTION 

Simulates  the  action  of  an  envelope  chip  rate  detector  with  complex  envelope  input 
Provides  optional  storage  of  intermediate  signals. 

MENU  PATH 

Detectors  >  Chip  Rate  >  Envelope,  core 
INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

write_signals  1  to  store  intermediate  signal,  0  otherwise 
f  ce  d  Receiver  center  frequency 

b  w  i  f  Receiver  bandwidth 

coIlect_time  Output  integration  time 

s  f r e  q  Sampling  frequency 

fees  Signal  center  frequency 

nsamples  Number  of  samples  in  simulation 

LIBRARY  /  FUNCTION 

delib/envel_chip_cor 


C-64 


ENVELOPE,  STEPPING  SEARCH  CHIP  RATE  DETECTOR 


DESCRIPTION 

Combines  a  envelope  chip  rate  detector  with  a  two-dimensional  search  that  scans  receiver 
center  frequency  and  chip  rate.  The  chip  rate  cycles  through  a  fixed  number  of  values 
specified  by  the  user,  while  the  receiver  center  frequency  keeps  incrementing  throughout  Ae 
simulation,  once  after  each  cycle  through  the  chip  rate,  so  the  number  of  center  frequencies 
used  depends  on  the  number  of  collects  simulated  and  the  number  of  chip  rates  used.  The 
user  controls  the  order  and  passband  edge  attenuation  of  the  receiver  filter  and  specifies  the 
name,f  ilename,  of  the  file  to  which  the  results  are  written. 

After  the  simulation,  the  user  obtains  a  Gnuplot  screen  plot  by  entering 

EnvSrch  filename 

in  an  operating  system  window.  The  plot  also  displays  the  name  of  the  underlying  data  file. 
To  obtain  a  hard  copy  in  addition  to  the  screen  plot,  the  user  enters 

EnvSrch  -h  filename 

If  the  default  printer  has  not  been  set  during  the  session,  the  user  may  have  to  specify  the 
printer  to  be  used  by  entering,  in  the  operating  system  window,  setenv  PRINTER 
followed  by  the  printer  name.  Quitting  the  graphics  window,  or  typing  a  keystroke  in  tiie 
operating  system  window  causes  the  graphical  output  to  disappear,  and  the  Gnuplot  session 
to  close. 

EnvSrch  is  a  Bourne  shell  script  which  selects  and  edits  the  file 

EnvSrch.  hard,  template  if  the  flag  -h  is  present,  and 

EnvSrch .  soft .  template  otherwise,  and  calls  Gnuplot  to  execute  the  resulting 
Gnuplot  script,  EnvSrch. gnu.  The  data  file  and  the  script  files  are  in  the  directory 
/spwdata/vmsign,  but  EnvSrch  can  be  invoked  from  any  directory  without 
specifying  a  path. 

MENU  PATH 

Interceptor  >  Detector  >  Chip  Rate  >  Envelope,  Step  Search 
INPUT  PORTS 

IN  Complex 

HOLD  Complex 

OUTPUT  PORTS 
None 


C-65 


ENVELOPE,  STEPPING  SEARCH,  CHIP  RATE  DETECTOR,  continued 


PARAMETERS 

start_fced 

Starting  receiver  center  frequency 

step_fced 

Receiver  center  frequency  increment 

bwif 

Receiver  bandwidth 

fchd 

Interceptor  chip  rate  estimate 

fiIter_order 

Receiver  filter  order 

edge_atten 

Receiver  passband  edge  attenuation,  dB  power 

stQre_signals 

1  to  store  intermediate  signals,  0  otherwise 

start_fchd 

Starting  chip  rate 

step_fchd 

Chip  rate  increment 

no_of_chip_rates 

Number  of  chip  rates  Chip  rate  increment 

collect^time 

Output  integration  time 

sfreq 

Sampling  frequency 

fees 

Signal  center  frequency 

USAGE  EXAMPLE 

SSyM>Denionstrations>Search  Chip  Rate  Detector 
LIBRARY  /  FUNCTION 

delib/envl_chip_srch 
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ENVELOPE  SPECTRUM 


DESCRIPTION 

Lowpass  filters  the  incoming  complex  signal;  buffers  the  signal  into  blocks;,  subsamples 
the  blocks;  forms  a  spectrum  by  performing  an  FFT  and  magnitude  squaring  the  output; 
averages  the  results  over  a  number  of  blocks,  and  stores  the  logarithm  of  the  results  for 
subsequent  plotting. 

After  the  simulation,  the  user  enters 

Spectrum  filename 

in  an  operating  system  window,  where  filename  is  the  data  file  name  selected  by  the 
user  in  the  spectrum  module.  TTie  result  is  a  Gnuplot  screen  plot  showing  spectral  strength 
as  a  function  of  frequency.  The  plot  also  displays  the  name  of  the  underlying  data  file.  To 
obtain  a  hard  copy  in  addition  to  the  screen  plot,  the  user  enters 

Spectrum  -h  filename 

If  the  default  printer  has  not  been  set  during  the  session,  the  user  may  have  to  specify  the 
printer  to  be  used  by  entering,  in  the  operating  system  window,  setenv  PRINTER 
followed  by  the  printer  name.  Quitting  the  graphics  window,  or  typing  a  keystroke  in  the 
operating  system  window  causes  the  graphical  output  to  disappear,  and  tiie  Gnuplot  session 
to  close. 

Spectrum  is  a  Bourne  shell  script  which  selects  and  edits  the  file 
Spectrum,  hard,  template  if  the  flag  -h  is  present,  and 

Spectrum .  sof t .  template  otherwise,  and  calls  Gnuplot  to  execute  the  resulting 
Gnuplot  script.  Spectrum .  gnu .  The  data  file  and  the  script  files  are  in  the  directory 
/spwdata/vmsign,  but  Spectrum  can  be  invoked  from  any  directory  without 
specifying  a  path. 

MENU  PATH 

None 

INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORTS 
None 
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ENVELOPE  SPECTRUM,  continued 

PARAMETERS 


start_time 

Starting  receiver  center  frequency 

bwif 

Receiver  bandwidth 

filter_order 

Receiver  filter  order 

edge_atten 

Receiver  passband  edge  attenuation,  dB  power 

subsampling 

Subsampling  order 

fft_si2e 

FFT  size,  power  of  2 

window_type 

Window  (blackman,  bartlett,  hamming,  banning,  rectangular) 

spec_avg_bIock 

Number  of  spectra  to  average 

center_freq 

Center  frequency 

freq_20om 

Frequency  zoom,  power  of  2 

sfreq 

Sampling  frequency 

fees 

Signal  center  frequency 

LIBRARY  /  FUNCTION 
delib/env_spectrum 


C-68 


ERROR  HANDLER 


DESCRIPTION 

Produces  an  error  message  in  a  pop-up  window  (whenever  IN  >  0.0) 
MENU  PATH 

Utilities  >  Logic  >  Error  Handler 
INPUT  PORT 

IN  Real 

OUTPUT  PORT 
None 

PARAMETERS 

parent_niodule  Parent  module  string 

message  Error  message  string 

LIBRARY  /  FUNCTION 
utlib/error_handler 

Codedin  C 


C-69 


EXCEEDANCE 


DESCRIPTION 

Treats  the  inputs  as  the  outputs  of  the  histogram  module,  that  is  as  the  underflow,  binned, 
and  overflow  counts  of  a  real  signal,  and  computes  the  exceedance.  The  exceedance  of  each 
bin  is  the  number  of  counts  larger  than  that  bin’s  lower  boundary,  that  is  the  number  of 
counts  in  that  bin  plus  the  number  of  counts  in  all  higher  bins.  The  overflow  is  treated  as  a 
bin  so  the  first,  or  highest  component  in  the  output  vector  is  equal  to  the  OVER  input.  The 
UNDER  input  is  passed  to  the  output  for  subsequent  processing. 

MENU  PATH 

Interceptor>  Metrics  >  Exceedance 
INPUT  PORTS 

OVER  Real 

BINS  Real  Vector 

UNDERJn  Real 

HOLD  Real 

OUTPUT  PORTS 

OUT  Real  vector,  of  dimension  equal  to  one  plus  the  dimension  of  the 

BINS  input 

UNDER_out  Real 

PARAMETERS 

BINS_IOVEC_LEN 

Number  of  components  in  the  input  vector  BINS 

LIBRARY  /  FUNCTION 
vulib/exceedance 

Codedin  C 


C-70 


ENVELOPE,  EFT  SEARCH  CHH*  RATE  DETECTOR 


DESCRIPTION 

Implements  an  envelope  chip  rate  detector  with  an  FFT  output  stage  to  permit  detection  at 
mdtiple  chip  rates.  In  addition  the  output  stage  provides  filtering,  subsampling  and  block 
averaging.  The  user  can  therefore  tailor  the  search  to  cover  efficiently  the  range  of  chip  rate 
uncertainty.  The  documentation  of  the  Envelope  Spectmm  Detector  provides  detailed 
information  on  the  filtering,  subsampling  and  block  averaging. 

MENU  PATH 

Interceptor  >  Detectors  >  Chip  Rate  >  Envelope,  FFT  Search 
INPUT  PORTS 


IN  Complex 


HOLD 

Real 

OUTPUT  PORT 

None 

PARAMETERS 

write^signals 

1  to  store  intermediate  signal,  0  otherwise 

feed 

Receiver  center  frequency 

bwif 

Receiver  bandwidth 

filter_order 

Receiver  filter  order 

edge_atten 

Receiver  passband  edge  attenuation,  dB  power 

store_signaIs 

1  to  store  intermediate  signals,  0  otherwise 

sfreq 

Sampling  frequency 

fees 

Signal  center  frequency 

start_time 

Time  to  begin  collecting  data  to  determine  averaged  spectmm 

spec_bwif 

Filter  width  (lowpass,  2-sided) 

subsampling 

Subsampling  order 

C-71 


ENVELOPE,  EFT  SEARCH  CHIP  RATE  DETECTOR,  continued 
fft_size  FFT  size,  power  of  2 

window_type  Window  (blackman,  bartlett,  hamming,  banning,  rectangular) 
spec_avg_block  Number  of  spectra  to  average 
freq_zooni  Frequency  zoom,  power  of  2 

USAGE  EXAMPLE 

SSVM>Demonstrations>Search  Chip  Rate  Detector 
LIBRARY  /  FUNCTION 
delib/env_chip_fft 


C-72 


FILE  STEPPER 


DESCRIPTION 

Produces  a  stepped  ou^ut  which  cycles  through  a  set  of  values.  The  user  stores  arbitrary 
values  in  a  file,  and  specifies  a  scaling  factor  and  an  offset.  The  user  also  controls  the  dwell 
time,  which  is  the  time  interval  between  changes  in  the  output 

MENU  PATH 

Utilities  >  Sources  >  File  Stepper 
INPUT 

HOLD 

OUTPUT  PORT 
OUT 

PARAMETERS 
dwell_time 
filename 
scale_factor 
offset 
sfreq 

LIBRARY  /  FUNCTION 
utlib/file_stepper 


Real 

Real 

Time  interval  between  changes  in  output 
Name  of  file  used  for  file  mode 
Scaling  applied  to  values  from  file 
Offset  applied  to  scaled  values  from  file 
Sampling  frequency 
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FILTERING 


DESCRIPTION 

Filters  the  complex  envelope  input  with  a  lowpass  filter  that  may  be  offset  from  the  signal 
center  frequency.  The  user  controls  the  filter  offset  and  bandwidth.  Ibe  filter  is  a 
ButterworA  HR  filter  with  an  adjustable  order  preset  at  20  and  an  adjustable  passband  edge 
attenuation  preset  to  3dB.  The  filtering  module  may  only  be  used  in  the  two  weak  signal 
operational  modes  of  the  Vulnerability  Metric. 

Communicators  adjust  the  signal  power  so  the  transmitted  power  is  constant,  independent 
of  the  filtering.  Direct  sim^ation  of  power  compensation  would  require  an  additional 
simulation  to  measure  the  impact  of  filtering  on  power.  To  maintain  computational 
efficiency,  the  Vulnerability  Metric  allows  the  filter  m^ule  to  reduce  the  transmitted  power, 
but  measures  and  stores  the  reduction  factor  on  disk  at  the  end  of  the  simulation.  In  both  of 
the  weak  signal  modes,  the  vulnerability  assessment  at  the  end  of  the  simulation  reads  the 
power  r^uction  factor  and  rescales  the  signal  strength.  From  the  user’s  perspective,  the 
vulnerability  evaluation  runs  exactly  as  if  transmitted  power  were  maintained  constant 

MENU  PATH 


Spread  Spectrum  >  Filtering 
INPUT  PORTS 


IN 

Complex 

HOLD 

Real 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

offset 

Offset  of  filter  passband  center  form  signal  center  frequency 

bw 

Filter  bandwidth 

fiIter_order 

Receiver  filter  order 

edge_atten 

Receiver  passband  edge  attenuation,  dB  power 

sfreq 

Sampling  frequency 

fees 

Signal  center  frequency 

nsamples 

Number  of  iterations  in  simulation 

LIBRARY  /  FUNCTION 


walib/filter 


C-74 


FIXED  EXCISION  CORE 


DESCRIPTION 

Performs  a  spectral  analysis  of  the  input  and  removes  the  strongest  components.  The  user 
specifies  the  length  of  the  vectors,  and  the  number  of  strongest  bins  to  be  excised. 

MENU  PATH 
None 

INPUT  PORTS 

X-re  Real  Vector 

X-im  Real  Vector 

OUTPUT  PORTS 

Y-re  Real  Vector 

Y-im  Real  Vector 

PARAMETERS 

blocksize  Number  of  points  used  in  the  FFTs 

nlargest  Number  of  bins  to  be  excised 

window  type  Type  of  window  used  in  the  forward  FFT 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
prlib/exfixed_cnip 

Coded  in  C 


C-75 


FREQUENCY  AGILITY 


DESCRIPTION 

Shifts  the  incoming  signal’s  frequency  and  phase  by  random  offsets  which  remain  constant 
over  time  intervals  c^ed  hops,  and  vary  pseudorandomly  between  hops.  Hop  length, 
generated  by  the  AGILITY  module,  also  vary  pseudorandomly,  and  are  rounded  to  die 
nearest  multiple  of  the  data  symbol  duration.  Typically  the  difference  between  the  minimum 
and  maximum  hop  lengths  exceeds  the  data  symbol  duration,  so  the  FREQUENCY 
AGILITY  module  performs  slow  frequency  hopping  with  fluctuation  hop  lengths. 

The  user  specifies  the  minimum  and  maximum  block  lengths;  the  number  of  frequency 
offsets;  the  maximum  one  sided  frequency  deviation;  and  whether  the  block  lengths  and 
block  boundaty  pulses  are  stored.  If  the  input  signal  has  constant  unit  magnitude,  the 
output  power  is  1. 

MENU  PATH 

Spread  Spectrum  >  Frequency  Agility 
INPUT  PORT 


IN 

Complex 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

min_duration 

Minimum  block  length 

max_duration 

Maximum  block  length 

nfreq 

Number  of  available  frequency  offsets 

maxfreqs 

Maximum  frequency  deviation 

store_hops 

1  to  store  block  lengths  and  block  boundary  pulses,  0  otherwise 

sfreq 

Sampling  frequency 

fda 

Data  rate 

LIBRARY  /  FUNCTION 


walib/freq_agility 
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FREQUENCY  DOUBLER  SYSTEM 

DESCRIPTION 

Uses  a  communication  signal  input  and  a  backgroimd  input  to  form  a  signal  present  input  to 
one  detector  and  a  signal  absent  input  to  another  detector.  The  detectors  are  each  fiequency 
doubler  detectors.  The  result  of  processing  the  signal  present  case  is  the  Hj  output  and  tiie 

result  of  processing  the  signal  absent  case  is  the  Hq  output.  Stores  detector  parameter 
values  in  files. 

MENU  PATH 

Interceptor  >  Detector  >  Freq  Multr  >  Doubler,  system 
INPUT  PORTS 


SIGNAL 

Complex 

BACKGROUND 

Complex 

OUTPUT  PORT 

HI 

Complex 

HO 

Complex 

PARAMETERS 

feed 

Receiver  center  frequency 

bwif 

Receiver  bandwidth 

fchd 

Interceptor’s  estimate  of  chip  rate 

wswgn 

1  for  weak  signal,  white  Gaussian  noise  mode,  0  otherwise 

nwssb 

1  for  nonweak  signal,  simulated  background  mode,  0  otherwise 

fees 

Signal  center  frequency 

sfreq 

Sampling  fiequency 

nsamples 

Number  of  samples  in  simulation 

fda 

data  rate 

coIlect_time 

Output  integration  time  in  detector  systems 

LIBRARY  /  FUNCTION 
delib/doubler 
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FREQUENCY  DOUBLER  CORE 

DESCRIPTION 

Simulates  the  action  of  an  frequency  doubler  detector  with  complex  envelope  input 
Provides  optional  storage  of  intermediate  signals. 

MENU  PATH 

Detectors  >  Freq  Multrs  >  Doubler,  core 
INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

write_signals  1  to  store  intermediate  signal,  0  otherwise 

feed  Receiver  center  frequency 

b  w  i  f  Receiver  bandwidth 

collect^time  Output  integration  time 

sfreq  Sampling  frequency 

fees  Signal  center  frequency 

nsamples  Number  of  samples  in  simulation 

LIBRARY  /  FUNCTION 
delib/doubIer_core 


C-78 


FREQUENCY  HOPPING 


DESCRIPTION 

Shifts  the  incoming  signal’s  frequency  and  phase  by  random  offsets.  The  offsets  are 
constant  during  intervals  defined  as  hops,  and  are  reselected  for  each  hop.  If  the  input 
signal  has  constant  unit  magnitude,  the  output  power  is  1. 

MENU  PATH 

Spread  Spectrum  >  Frequency  Hopping 
INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

hop_rate  Hop  rate 

nfreq  Number  of  available  frequency  offsets 

maxfreqs  Maximum  frequency  deviation 

write_freqs  1  to  store  frequency  offsets 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 

walib/incoh_freq_hop 


C-79 


FREQUENCY  HOPPING  DOWNCONVERTER 


DESCRIPTION 

Coherently  reverses  the  effect  of  the  Frequency  Hopping  and  Frequency  Agility  modules. 
MENU  PATH 

Communicator/Preprocessor  >  Freq  Hop  Downconvert 
INPUT  PORTS 


IN 

Complex 

DELAYJN 

Real 

HOLD 

Real 

OUTPUT  PORTS 

OUT 

Complex 

DELAY_OUT 

Real 

PARAMETERS 

sfreq 

Sampling  frequency 

LIBRARY  /  FUNCTION 


prlib/fh_downconvert 


C-80 


FSK  DEMODULATOR 


DESCRIPTION 

A  demodulator  of  M-ary  frequency  shift  keying. 
MENU  PATH 

Communicator  >  Demodulator  >  FSK 
INPUT  PORT 


IN 

Complex 

DELAYJN 

Real 

OUTPUT  PORT 

OUT 

Real 

DELAY_OUT 

Real 

PARAMETERS 

m 

Alphabet  size 

fsep 

Tone  separation 

fda 

Bit  rate 

sfreq 

Sampling  frequency 

LIBRARY  /  FUNCTION 
modlib/fsk_demod 


C-81 


FSK  MODULATOR 


DESCRIPTION 

A  M-ary  frequency  shift  keyed  modulator.  Either  coherent  or  incoherent  modulation  may  be 
chosen. 


MENU  PATH 

Communicator  >  Modulator  >  FSK 
INPUT  PORT 


IN 

Real 

DELAY^IN 

Real 

OUTPUT  PORT 

OUT 

Complex 

DELAY  OUT 

Real 

PARAMETERS 


m 

fsep 

coherence_fIag 

fda 

sfreq 


Alphabet  size 
Tone  separation 

Flag  to  determine  whether  coherent  FSK  is  chosen.  1  of  coherent, 
0  for  incoherent. 

Bit  rate 

Sampling  frequency 


LIBRARY  /  FUNCTION 


modlib/fsk_mod 


C-82 


FSK  SOURCE 


DESCRIPTION 

Outputs  a  narrowband  complex  envelope  resulting  from  interleaving  a  pseudorandom  bit 
source,  and  using  the  result  as  the  data  stream  for  a  frequency  shift  keying  modulation.  The 
user  sets  the  numbers  of  columns  and  rows  to  define  the  interleaving,  and  accepting  the 
default  value  of  1  for  both  produces  a  non-interleaved  data  stream.  The  user  also  sets  the 
alphabet  size,  that  is  the  number  of  frequencies  to  be  used  in  the  modulation;  the  frequency 
separation,  and  whether  the  modulation  is  coherent  or  not 

MENU  PATH 

Communicator  >  Modulated  Source  >  FSK 
INPUT  PORT 
None 

OUTPUT  PORT 

OUT  Complex 

DELAY_OUT  Real 

PARAMETERS 

ncols  Interleaver  columns 

n r o ws  Interleaver  rows 

m  Alphabet  size 

fsep  Tone  Separation 

coherence_flag  Flag  to  determine  whether  coherent  FSK  is  chosen,  1  for  coherent,  0 

for  incoherent. 

fda  Bit  rate 

sfreq  Sampling  rate 

USAGE  EXAMPLES 

SSVM>Demonstrations>  Correlation  Radiometer 
SSVM>Demonstrations>  Scanner 

LIBRARY  /  FUNCTION 

modulib/fsk  source 


ESTIMATED  OUTPUT  SNR 


DESCRIPTION 

Estimates  the  output  signal  to  noise  ratio 


(S/JV)o  = 


Var^z  |//i} 


The  denominator  is  a  variance,  and  the  first  output  of  the  variance  estimators  is  zero,  so  the 
module’s  first  output  sample  should  be  ignored. 

MENU  PATH 

Interceptor  >  Metrics  >  Estimated  Output  SNR 


INPUT  PORTS 

HI 

Complex 

HO 

Complex 

HOLD 

Real 

OUTPUT  PORTS 

OUT 

Complex 

PARAMETERS 

None 


LIBRARY  /  FUNCTION 
vulib/full  snro 


C-84 


GAMMA  NOISE 


DESCRIPTION 

Produces  a  discrete  complex  envelope  model  of  continuous  bandpass  white  noise  whose 
magnitude  is  distributed  as  a  Gamma  random  variable,  and  whose  phase  is  uniformly 
distributed.  The  marginal  distributions  of  the  real  and  imaginaiy  components  are  not 
expressible  as  elementary  fxmctions,  but  decay  exponentially.  The  probability  density  of  the 
magnitude  is  given  by 


xe-^ 


The  noise  power  in  both  the  Gaussian  and  Gamma  Noise  modules  is  parametrized  by  the 
communicator’s  input  signal  to  noise  ratio,  which  is  computed  for  the  case  where  only  one 
noise  module  is  present.  Using  more  than  one  noise  module  at  a  time  requires  manually 
determining  the  input  signal  to  noise  ratio. 

MENU  PATH 


Background  >  Noise  >  Gamma 


INPUT  PORT 
IN 

OUTPUT  PORTS 
OUT 

PARAMETERS 

dB_in 

sfreq 

fda 

nsamples 


None 

Complex 

Input  signal  to  noise  ratio,  for  the  case  where  no  other  noise  module 
is  in  use 

Sampling  frequency 
Data  rate 

Number  of  samples  in  the  simulation 


LIBRARY  /  FUNCTION 
bklib/whi_gam_nse 


HETERODYNE 


DESCRIPTION 

Multiplies  the  input  by  a  complex  tine. 

MENU  PATH 

Utilities  >  Signal;  Processing  >  Heterodyne 
INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

het^freq  Heterodyne  frequency 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/heterodyne 


C-86 


HF  INTERFERENCE 


DESCRIPTION 

Adds  to  the  signal  an  interference  background  generated  from  a  stored  configuration  chosen 
to  be  characteristic  of  typical  HF  operating  environments.  The  user  selects  the  center  and 
width  of  the  operating  band,  and  the  m^ule  generates  the  interferers  which  lie  in  tiie 
selected  band. 

The  mathematical  form  of  the  interference  the  same  as  that  detailed  for  Random  Interfefence. 
MENU  PATH 

Background  >  Interference  >  HF 
INPUT  PORTS 
None 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

if_cnet_freq  Configuration  center  frequency 

if_bandwidtli  Configuration  one-sided  bandwidth 

simtime  Simulation  time 

sfreq  Sampling  frequency 

fees  Signal  center  frequency 

LIBRARY  /  FUNCTION 
bklib/ifhf 

Coded  in  C 


C-87 


HISTOGRAM 


DESCRIPTION 

Accumulates  a  histogram  of  the  incoming  signal.  The  user  specifies  the  lowest  and  highest 
values  captured  and  the  number  of  bins  .into  which  the  intervening  interval  is  divided.  The 
module  counts  the  number  of  input  samples  falling  into  each  bin,  as  well  as  the  number 
below  the  lowest  captured  value  and  the  number  above  the  highest  captured  value. 

MENU  PATH 


Interceptor  >  Metrics  >  Histogram 


INPUT  PORTS 
IN 

HOLD 

OUTPUT  PORTS 
OVER 
OUT 
UNDER 
PARAMETERS 
min 
max 


Complex 

Real 

Real 

Real  vector,  of  dimension  equal  to  the  number  of  bins 
Real 

Lowest  value  captured 
Highest  value  captured 


DEFAULT_VECLEN 

Number  of  bins 


LIBRARY  /  FUNCTION 
vulib/histogram 

Coded  in  C 


C-88 


HOLD  SWITCH 


DESCRIPTION 

If  the  selection  index  parameter  is  zero,  sets  all  ou^uts  to  TRUE,  Otherwise  sets  all  outputs 
to  TRUE  except  one,  which  is  identified  by  a  selection  index  parameter,  "^e  module  sets 
the  selected  port  to  FALSE.  Another  parameter  specifies  the  number  of  options,  that  is  the 
number  of  outputs  available  for  the  selection,  plus  one  for  the  case  where  none  is  selected. 
The  module  sets  the  value  of  the  both  parameters,  and  terminates  with  an  error  message  if 
either  is  out  of  range.  Setting  the  parent  module  parameter  permits  the  error  message  to  give 
the  containing  module  of  the  instance  that  has  an  error. 

MENU  PATH 


Utilities  >  Logic  >  Hold  Switch 
INPUT  PORT 
None 

OUTPUT  PORT 

OUT_l  Real 

OUT_2  Real 

OUT_3  Real 

OUT_4  Real 

OUT_5  Real 

PARAMETERS 

no_of_options  Number  of  options  (1-6) 

index  Selection  index,  between  0  and  (number  of  options  - 1) 

LIBRARY  /  FUNCTION 
utlib/switch  hold 


C-89 


INCOMPLETE  RUN 


DESCRIPTION 

Outputs  TRUE  until  the  number  of  executed  iterations  exceeds  the  total  number  of  iterations 
in  the  simulation  minus  2,  when  the  output  becomes  False. 

MENU  PATH 


Utilities  >  Incomplete  Run 
INPUT  PORT 

HOLD  Real 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 


sfreq 

collect_time 

Collects 


Simulation  sampling  frequency 
Detector  collect  time 
Number  of  collects 


LIBRARY  /  FUNCTION 


utlib/incomplete_run 


INTERCEPT  COUPLER 


DESCRIPTION 

Uses  the  transmitted  signal  and  background  noise  to  form  the  signal  plus  noise  (HI) 
and  noise  alone  ^0)  hypothetical  waveforms  for  input  to  the  intercept  detectors. 

MENU  PATH 

Channel  >  Couplers  >  Interceptor 
INPUT  PORT 

SIGNAL  Real 

BACKGROUND  Real 
OUTPUT  PORT 

HI  Real 

HO  Real 

PARAMETERS 
None 

USAGE  EXAMPLES 

SSVM>Demonstrations>Search  Chip  Rate  Detector 
SSVM>Denionstrations>Scanner 

LIBRARY  /  FUNCTION 

chlib/intcpt_coupler 


C-91 


INTERCEPTOR’S  ADAPTIVE  THRESHOLD  EXCISION  SYSTEM 
DESCRIPTION 

Applies  adaptive  threshold  excision  to  both  the  signal  present  Hj  and  signal  absent  Hq 
waveforms. 

MENU  PATH 

Interceptor  >  Preprocessor  >  Adaptive  Threshold  Excision 
INPUT  PORTS 

SIGNAL  Complex 

BACKGROUND  Complex 

OUTPUT  PORTS 

H 1  Complex 

H  0  Complex 

PARAMETERS 

exceedfactor  Exceedance  factor 

nexempt  Number  of  points  excluded  in  the  calculation  of  mean  square 

magnitude 

blocksize  Number  of  points  used  in  the  FFTs 

window_type  Type  of  window  used  in  the  forward  FFT 

sfreq  Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 
prlib/exadapt_int 


C-92 


INTERCEPTOR’S  FIXED  EXCISION  SYSTEM 


DESCRIPTION 

Applies  fixed  excision  to  both  the  signal  present  H  j  and  signal  absent  Hq  waveforms. 
MENU  PATH 

Interceptor  >  Preprocessor  >  Fixed  Excision 
INPUT  PORTS 

SIGNAL  Complex 

BACKGROUND  Complex 
OUTPUT  PORTS 

H 1  Complex 

HO  Complex 

PARAMETERS 

nlargest  Number  of  bins  to  be  excised 

blocksize  Number  of  points  used  in  the  FFTs 

window_type  Type  of  window  used  in  the  forward  FFT 
sfreq  Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 
prlib/exfixed_int 


C-93 


INTERCEPTOR’S  WHITENING  EXCISION  SYSTEM 


DESCRIPTION 

Applies  whitening  excision  to  both  the  signal  present  and  signal  absent  Hq  waveforms. 

MENU  PATH 

Interceptor  >  Preprocessor  >  Whitening  Excision 
INPUT  PORTS 

SIGNAL  Complex 

BACKGROUND  Complex 

OUTPUT  PORTS 

H 1  Complex 

HO  Complex 

PARAMETERS 

blocksize  Number  of  points  used  in  the  FFTs 

window_type  Type  of  window  used  in  the  forward  FFT 

sfreq  Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 
p  rli  b/exwhite_int 


C-94 


INTERCEPTOR’S  COMPOSITE  PREPROCESSOR 

DESCRIPTION 

Optionally  applies  one  of  three  frequency  based  excision  techniques:  fixed  adaptive 
threshold,  or  whitening. 

MENU  PATH 

Interceptor  >  Preprocessor  >  Composite 
INPUT  PORTS 

SIGNAL  Complex 

BACKGROUND  Complex 
OUTPUT  PORTS 

H 1  Complex 

H  0  Complex 

PARAMETERS 

index  Selection  index  (0-no  excision,  1-  Fixed,  2-Adaptive,  3- Whitening) 

sfreq  Simulation  sampling  frequency 

Additional  parameters  depending  in  index:  see  “INTERCEPT  FIXED  EXCISION”, 
“INTERCEPT  ADAPTIVE  EXCISION”,  and  “INTERCEPT  WHITENING  EXCISION” 

LIBRARY  /  FUNCTION 

prlib/intercept 


C-95 


INTERLEAVER 


DESCRIPTION 

Implements  a  block  interleaver,  which  permutes  data  by  sequentially  loading  an  NxM  matrix 
with  input  bits  row  by  row,  and  reading  the  bits  out  column  by  column. 

MENU  PATH 

Communicator  >  Modulator  >  Interleaver 
INPUT  PORTS 

IN  Real 

HOLD  Real 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 

ncols  Number  of  columns  in  interleaver  matrix 

n  r o  ws  Number  of  rows  in  interleaver  matrix 

fda  Data  rate 

sfreq  Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 
modlib/interleaver 


C-96 


TRANSPOSE 


DESCRIPTION 

Implements  the  core  a  block  interleaver,  permuting  data  by  sequentially  loading  an  NxM 
matrix  with  input  bits  row  by  row,  and  reading  the  bits  out  column  by  column.  This  inner 
block  is  called  by  both  the  top  level  interleaver  and  deinterleaver. 

MENU  PATH 

None 

INPUT  PORTS 

IN  Real 

HOLD  Real 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 

n  c  o  1  s  Number  of  columns  in  interleaver  matrix 

n  ro  ws  Number  of  rows  in  interleaver  matrix 

initial_value  Initial  value  in  interleaver  matrices 

LIBRARY  /  FUNCTION 
modlib/transpose 


C-97 


JITTER 


DES^JRIPTION 

Produces  trigger  pulses  at  randomly  occurring  time  steps.  The  random  timing  of  the  zeros  is 
constrained:  referring  to  the  intervd  between  successive  zeros  as  a  chip,  there  are  always  an 
integer  iiumber  of  chips  in  a  data  symbol,  and  the  chip  rate  only  changes  at  data  symbol 
boundaries.  The  instantaneous  chip  rate  is  generated  randomly  by  truncating  a  continuous 
uniform  random  variable  to  the  nearest  integer  multiple  of  the  ^ta  rate.  The  user  specifies 
the  mean  and  maximum  deviation  of  the  continuous  random  variable,  with  the  mean 
parametrized  as  a  multiple  of  the  data  rate,  and  the  size  of  the  chip  rate  fluctuations 
papnetrized  as  a  percentage  of  the  average  rate.  The  user  also  controls  the  interval  between 
chip  rate  changes,  parametrized  as  a  multiple  of  the  data  symbol  length,  and  the  user  may 
write  to  disk  the  instantaneous  chip  rate  and  the  module  output. 

The  terms  “jitter”  and  “dither”  are  synonymous. 

MENU  PATH 


Spread  Spectrum  >  Jitter 


INPUT  PORT 

None 

OUTPUT  PORT 

OUT 

Real 

PARAMETERS 

nchav 

Average  number  of  chips  per  symbol 

nsybj 

Number  of  symbols  between  chip  rate  changes 

hw 

Maximum  allowed  deviation  of  chip  rate  from  average  chip  rate 

writejitter 

1  to  write  instantaneous  chip  rate  and  chip  release  impulses  to  files,  0 
otherwise 

fda 

Data  rate 

sfreq 

Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 


walib/jitter 


C-98 


JITTERED  8ary  DS  SPREADER 


DESCRIPTION 

Multiplies  the  complex  input  by  a  8-ary  complex  pseudonoise  sequence  with  a  fluctiwtion 
chip  rate.  The  user  specifies  the  average  chip  rate  as  a  multiple  of  the  data  rate,  the  size  of 
the  chip  rate  fluctuations  as  a  percentage,  and  the  interval  between  chip  rate  changes  as  a 
multiple  of  the  data  symbol  length.  Optionally  stores  the  instantaneous  chip  rate  and  the 
sequence  of  impulses  that  trigger  each  chip.  If  the  input  signal  has  constant  unit  magnitude, 
the  output  power  is  1. 

MENU  PATH 

Spread  Spectrum  >  Jittered  Direct  Sequence  >  8ary 
INPUT  PORT 

IN  Complex 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

nchav  Average  number  of  chips  per  symbol 

nsybj  Number  of  symbols  between  chip  rate  changes 

h  w  Maximum  allowed  deviation  of  chip  rate  from  average  chip  rate 

write_jitter  1  to  write  instantaneous  chip  rate  and  chip  release  impulses  to  files,  0 

otherwise 

f  da  Data  rate 

sfreq  Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 
walib/j_8ary_ds 


C-99 


JITTERED  BINARY  DS  SPREADER 


DESCRIPTION 

Mdtiplies  the  complex  input  by  a  binary  complex  pseudonoise  sequence  with  a  fluctuation 
chip  rate.  The  user  specifies  the  average  chip  rate  as  a  multiple  of  the  data  rate,  the  size  of 
the  chip  rate  fluctuations  as  a  percentage,  and  the  interval  between  chip  rate  changes  as  a 
multiple  of  the  data  symbol  length.  Optionally  stores  the  instantaneous  chip  rate  and  the 
sequence  of  impulses  that  trigger  each  chip.  If  the  input  signal  has  constant  unit  magnitnHp 
the  output  power  is  1 . 

Delays  the  output  signal  by  the  input  delay  value  for  automatic  synchronization. 

Outputs  the  spreading  sequence  as  a  global  output  connector  called  ‘code’. 

MENU  PATH 

Spread  Spectrum  >  Jittered  Direct  Sequence  >  Binary 
INPUT  PORT 


IN 

Complex 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

nchav 

Average  number  of  chips  per  symbol 

nsybj 

Number  of  symbols  between  chip  rate  changes 

hw 

Maximum  allowed  deviation  of  chip  rate  from  average  chip  rate 

write_jitter 

1  to  write  instantaneous  chip  rate  and  chip  release  impulses  to  files,  0 
otherwise 

fda 

Data  rate 

sfreq 

Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 


walib/j_binary_ds 


C-100 


JITTERED  QUATERNARY  DS  SPREADER 


DESCRIPTION 

Multiplies  the  complex  input  by  a  quaternary  complex  pseudonoise  sequence  with  a 
fluctuation  chip  rate.  The  user  specifies  the  average  chip  rate  as  a  multiple  of  the  data  rate, 
the  size  of  the  chip  rate  fluctuations  as  a  percentage,  and  the  interval  between  chip  rate 
changes  as  a  multiple  of  the  data  symbol  length.  Optionally  stores  the  instantaneous  chip 
rate  and  the  sequence  of  impulses  that  trigger  each  chip.  If  the  input  signal  has  constant  unit 
magnitude,  the  output  power  is  1. 

Delays  the  output  signal  by  the  input  delay  value  for  automatic  synchronization. 

Outputs  the  spreading  sequence  as  a  global  output  connector  called  ‘code’. 

MENU  PATH 

Spread  Spectrum  >  Jittered  Direct  Sequence  >  Quaternary 
INPUT  PORT 


IN 

Complex 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

nchav 

Average  number  of  chips  per  symbol 

nsybj 

Number  of  symbols  between  chip  rate  changes 

hw 

Maximum  allowed  deviation  of  chip  rate  from  average  chip  rate 

write_jitter 

1  to  write  instantaneous  chip  rate  and  chip  release  impulses  to  files,  0 
otherwise 

fda 

Data  rate 

sfreq 

Simulation  sampling  fi-equency 

LIBRARY  /  FUNCTION 


walib/j_quatary_ds 


C-101 


LEADING  EDGE  MOD  N  IMPULSE  TRAIN 


DESCRIPTION 

Outputs  an  impulse  train  with  output  values  equal  to  zero,  except  for  one  sample  each  period 
when  the  output  is  one.  Differs  from  the  SPW  Mod  N  Impulse  Train  in  that  the  output  on 
the  first  iteration,  sample  number  zero,  is  one. 

MENU  PATH 

Utilities  >  Source  >  Lead  Mod  N  Imp  Train 
INPUT  PORT 
None 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 

period  Period  of  the  impulse  train  measured  in  samples.  If  the  parameter 

value  is  not  an  integer  the  period  used  is  the  nearest  integer. 

LIBRARY  /  FUNCTION 
utlib/Iead_train 


C-102 


LOWPASS/BANDPASS  FILTER 


DESCRIPTION 

Functions  as  either  a  lowpass  or  a  bandpass  filter,  with  the  user  making  the  choicfe  by 
setting  an  integer  flag  to  0  or  1  respectively. 

MENU  PATH 

Utilities  >  Signal  Processing  >  Low/Bandp  Filter 
INPUT  PORT 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

control  0  for  lowpass  or  1  for  bandpass 

1  p_edge  Lowpass  edge  frequency 

bp_  center  Bandpass  center  frequency 

b  p_b w  Bandpass  bandwidth  frequency 

sfreq  Simulation  sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/lp_bp_filter 


C-103 


LOWPASS  CYCLIC  FEATURE  DETECTOR  SYSTEM 


DESCRIPTION 

Uses  a  communication  signal  input  and  a  background  input  to  form,  a  signal  present  input 
to  one  detector  and  a  signal  absent  input  to  another  detector.  The  detectors  are  each  cyclic 
feature  detectors  for  lowpass  cyclic  frequencies.  The  result  of  processing  the  signal  present 
case  is  the  H  2  output  and  the  result  of  processing  the  signal  absent  case  is  the  Hq  output. 

Stores  detector  parameter  values  in  files. 

MENU  PATH 


Interceptor  >  Detector  >  Cyclic  Feature  >  Lowpass,  system 
INPUT  PORTS 


SIGNAL 

Complex 

BACKGROUND 

Complex 

OUTPUT  PORTS 

HI 

Complex 

H2 

Complex 

PARAMETERS 

cyclic_freq 

Cyclic  frequency 

filter_zone 

Spectral  frequency  zone:  0  for  lowpass,  1  for  bandpass 

lp_edge 

Lowpass  edge  frequency 

bp_center 

Bandpass  center  frequency 

bp_bw 

Bandpass  bandwidth  frequency 

delay_time 

Delay  in  delay,  conjugate  and  multiply  section  of  detector  systems 

coIIect_time 

Output  integration  time  in  detector  systems 

wswgn 

1  for  weak  signal,  white  Gaussian  noise  mode,  0  otherwise 

nwssb 

1  for  nonweak  signal,  simulated  background  mode,  0  otherwise 

C-104 


LOWPASS  CYCLIC  FEATURE  DETECTOR  SYSTEM,  continued 

sfreq  Sampling  frequency 

nsamples  Number  of  samples  in  simulation 

fda  data  rate 

LIBRARY  /  FUNCTION 
delib/cyc_featl 


C-105 


LOWPASS  CYCLIC  FEATURE  DETECTOR  CORE 


DESCRIPTION 

Simulates  the  action  of  a  cyclic  feature  detector  for  lowpass  cyclic  frequencies  envelope 
input.  Provides  optional  storage  of  intermediate  signals. 

MENU  PATH 

Interceptor  >  Detector  >  Cyclic  Feature  >  Lowpass,  core 
INPUT  PORTS 


IN 

Complex 

HOLD 

Real 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

cyclic_freq 

Cyclic  frequency 

filter^zone 

Spectral  frequency  zone:  0  for  lowpass,  1  for  bandpass 

lp_edge 

Lowpass  edge  frequency 

bp_center 

Bandpass  center  frequency 

bp_bw 

Bandpass  bandwidth  frequency 

delay_time 

Delay  in  delay,  conjugate  and  multiply  section  of  detector  systems 

collect_time 

Output  integration  time  in  detector  systems 

sfreq 

Sampling  frequency 

nsamples 

Number  of  samples  in  simulation 

LIBRARY  /  FUNCTION 


delib/cyc__featl_core 


C-106 


LOWPASS  CYCLIC  SPECTRUM 


DESCRIPTION 

Computes  the  lowpass  component  of  the  cyclic  spectmm  of  a  complex  signal  SDE.  Tlie 
lowpass  cyclic  spectmm  is  computed  at  a  user-specified  cyclic  frequency  and  for  spectral 
frequencies  between/g  -/sa/^  and/g  +/sa/2  where/g^  is  the  sampling  frequency  and/ 

g  is  the  signal  center  frequency. 

The  estimation  of  cyclic  spectra  combines  tiie  analysis  of  Appendix  E  with  functions  in 
COMDISCO’s  Signal  Display  Editor,  SDE.  The  final  implementation  in  SDE  is  a  maao 
which  has  performed  successfully  in  test,  although  errors  in  COMDISCO’s  documentation 
hinder  an  analytic  proof  of  correctness. 

MENU  PATH 

SDE:  SYSTEM/MACRO/EXECUTE/vmlib/CyclicJowpass 
INPUT  PORTS 
OUT 

OUTPUT  PORTS 
HI 
H2 

PARAMETERS 
numsigs 
alpha 
samp_freq 
numpts 
numffpts 

LIBRARY  /  FUNCTION 

spwdata/sdemacro/vmilb/Cyclic_low.mac 


Real 

Complex 

Complex 

Display  number  of  signal  to  analyze 
Cyclic  frequency  desired 
Sampling  frequency 

Number  of  samples  to  use  in  the  computation 
Number  of  points  in  the  FFT  (must  be  a  power  of  two) 


C-107 


LOWPASS  (KAISER)  FILTER 


MENU  PATH 

Spread  Spectrum  >  Lowpass  (Kaiser)  Filter 
INPUT  PORTS 

IN  Real 

DELAY_IN  Real 

OUTPUT  PORTS 

OUT  Complex 

DELAY_OUT  Real 

PARAMETERS 

f_c  Cutoff  frequency 

delta_f  Transition  region  width 

a_r  Stopband  attenuation 

s  f r e  q  Sampling  frequency 

LIBRARY  /  FUNCTION 
walib/kaiser_filter 
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MAGNITUDE  SQUARED 


DESCRIPTION 

Outputs  the  magnitude  squared  of  the  complex  input. 

MENU  PATH 

Utilities  >  Signal  Processing  >  Magnitude  Sqaured 


INPUT  PORTS 

IN 

Complex 

HOLD 

Real 

OUTPUT  PORTS 

OUT 

Real 

PARAMETERS 

None 

LIBRARY  /  FUNCTION 
utlib/mag_sqd 


ft 
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MAKE  ROC 


DESCRIPTION 


Treats  each  pair  of  inputs  as  the  outputs  of  the  exceedance  module,  that  is  as  a  vector  of 
exceedances  and  an  underflow,  and  uses  each  pair  of  inputs  to  estimate  the  probability 
distribution  function  of  the  corresponding  real  siganl.  There  are  two  pairs  of  inputs,  one  for 
the  signal  present  and  the  other  for  the  signal  absent  simulation,  and  the  modide  combines 
the  two  probability  distribution  estimates  to  form  probability  of  false  alarm,  probability  of 
detection,  pairs.  The  user  specifies  a  name,  filename,  in  the  parameter  table,  and  die 
module  stores  probability  paris  in  that  file,  in  the  directory  /pspwdata/vmsigs.  To 
control  statistical  signifigance,  the  user  specifies  a  minimum  exeedance  count  number  for 
the  signal  absent  input,  and  the  module  terminates  the  tail  of  the  ROC  at  the  bin  where  the 
number  of  counts  goes  below  the  minimum.  In  addition,  the  module  writes  the  complete  set 
of  probability  pairs  to  a  second  file  with  a  name  formed  by  adding  the  extension  .  all  to 
filename . 

MENU  PATH 
None 

INPUT  PORTS 


H1_UNDER 

H1_BINS 

H0_UNDER 

HO-BINS 

HOLD 


Real 

Real  vector 
Real 

Real  vector,  of  dimension  equal  to  dimension  of  H1_BINS 
Real 


OUTPUT  PORTS 


None 

PARAMETERS 

DEFAULT_VECLEN 

Number  of  components  in  the  input  vectors,  equal  to  the  number  of 
probability  pairs  stored  in  the  file  filename .  all . 
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MAKE  ROC,  continued 


HO_threshold  Minimum  number  of  signal  absent  exceedance  counts  for  bins  to  be 

included  in  the  file  filename 

filename  Name  of  file  where  the  module  stores  probability  pairs  calculated 

from  bins  with  signal  absent  counts  exceeding  the  minimum. 

LIBRARY  /  FUNCTION 
Yulib/store  roc 


Coded  in  C 


MMSE  WHITENING  FILTER 

DESCRIPTION 

Implements  an  adaptive  whitening  matched  filter,  which  acts  to  excise  interference  by 
estimating  the  signal  plus  noise  power  spectral  density  and  inverting  it  in  the  fi’equency 
domain.  Rather  than  actually  excising  interferers,  this  module  deemph^asizes  the  frequency 
bands  in  which  they  fall.  Blocking  effects  are  mitigated  by  overlapping  successive  blocks 
to  which  a  complementary  taper  has  been  applied. 

MENU  PATH 

Communicator  >  Preprocessor  >  MMSE 
INPUT  PORTS 

IN  Complex 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

blocksize  Length  of  FFT 

average_type  Periodogram  averaging  type,  either  ‘sliding_block’  or  ‘decaying’ 

averaging_block  Length  of  sliding  window  of  ’sliding_block’  averaging  type 

danielle_param  Width  of  frequency  domain  smoothing  window 

alpha  Decay  parameter  of  ‘decaying’  average  type 

window_type  Tapering  window  applied  to  lapped  reconstraction  blocks,  either 

‘banning’  or  ‘bartlett’  are  recommended 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
prlib/lapped_adapt 
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MODE  1 


DESCRIPTION 

Estimates  the  full  output  SNR,  defined  in  Appendix  F  as 


(S/N)o  = 


Varjz  |//i} 


and  writes  the  final  estimate  to  a  file.  The  user  specifies  the  filename  as  a  parameter. 

MENU  PATH 

Interceptor  >  Metrics  >  Mode  1 
INPUT  PORTS 


HI 

Complex 

HO 

Complex 

HOLD 

Real 

OUTPUT  PORTS 
None 

PARAMETERS 

filename  Filename  to  store  results 

LIBRARY  /  FUNCTION 
vulib/mode_l 
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MODE  2 


DESCRIPTION 

Estimates  the  full  output  SNR,  per  collected  bit,  and  writes  the  result  to  a  file.  The  full 
output  SNR  is  defined  as 


(S/iV)o  = 


E{z\H,]-E{z\Haf 
Var^z  |//i} 


and  the  number  of  collected  symbols  is  the  product  of  the  data  rate  and  the  collect  time. 


"^e  Vulnerability  Metric  extrapolates  the  result  to  a  range  of  collect  times.  After  the 
simulation,  a  window  appears  prompting  the  user  to  obtain  a  Gnuplot  screen  plot  by 
entering 


Plot2  filename 

in  an  operating  system  window,  where  filename  is  the  data  file  name  selected  by  the 
user.  The  plot  also  displays  the  name  of  the  underlying  data  file.  To  obtain  a  hard  copy  in 
addition  to  the  screen  plot,  the  user  enters 

Plot2  -h  filename 

If  the  default  printer  has  not  been  set  during  the  session,  the  user  may  have  to  specify  the 
printer  to  be  used  by  entering,  in  the  operating  system  window,  setenv  PRINTER 
followed  by  the  printer  name.  Quitting  the  graphics  window,  or  typing  a  keystroke  in  the 
operating  system  window  causes  the  graphical  output  to  disappear,  and  the  Gnuplot  session 
to  close. 

Plot  2  is  a  Bourne  shell  script  which  selects  and  edits  the  file 
Plot2 . hard .  template  if  the  flag  -h  is  present,  and 
Plot2  .  soft  .template  otherwise,  and  calls  Gnuplot  to  execute  the  resulting 
Gnuplot  script,  Plot2  .  gnu .  The  data  file  and  the  script  files  are  in  the  directory 
/spwdata/vmsign,  but  Plot 2  can  be  invoked  fi’om  any  directory  without 
specifying  a  path. 

MENU  PATH 

Interceptor  >  Metrics  >  Mode  2 
INPUT  PORTS 

H 1  Complex 
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MODE  2,  continued 


HO  Complex 

HOLD  Real 

OUTPUT  PORTS 
None 

PARAMETERS 

filename  Filename  to  store  results 

sfreq  Sampling  frequency 

fda  Data  rate 

collect_time  Collect  time 

collects  Number  of  collects 

LIBRARY  /  FUNCTION 
vulib>mode  2 
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MODE  3 


DESCRIPTION 


Estimates  the  full  output  signal  noise  ratio  defined  as 

Var{z  |//o} 


and  stores  the  result  in  a  file  as  a  coefficient  which  permits  extrapolation  of  the  result  to  a 
range  of  collect  times  and  a  range  of  interceptor’s  input  signal  power.  The  extrapolation  to  a 
range  of  signal  powers  is  accurate  only  for  low  signal  powers,  and  the  module  ensures  that 
results  are  presented  over  a  range  where  they  are  accurate. 

Plots  filename 

in  an  operating  system  window,  where  filename  is  the  data  file  name  selected  by  the 
user.  The  plot  also  displays  the  name  of  the  underlying  data  file.  To  obtain  a  hard  copy  in 
addition  to  the  screen  plot,  the  user  enters 

Plot  3  -h  filename 

If  the  default  printer  has  not  been  set  during  the  session,  the  user  may  have  to  specify  the 
printer  to  be  used  by  entering,  in  the  operating  system  window,  setenv  PRINTER 
followed  by  the  printer  name.  Quitting  the  graphics  window,  or  typing  a  keystroke  in  the 
operating  system  window  causes  the  graphical  output  to  disappear,  and  the  Gnuplot  session 
to  close. 

Plots  is  a  Bourne  shell  script  which  selects  and  edits  the  file 
Plot3  .hard. template  if  the  flag  -h  is  present,  and 
Plots .  soft .  template  otherwise,  and  calls  Gnuplot  to  execute  the  resulting 
Gnuplot  script.  Plot 2  .gnu .  The  data  file  and  the  script  files  are  in  the  directory 
/spwdata/vmsign,  but  Plots  can  be  invoked  from  any  directory  without 
specifying  a  path. 

MENU  PATH 

Interceptor  >  Metrics  >  Mode  3 
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MODE  3,  continued 


INPUT  PORTS 
SIG 
HO 

HOLD 

OUTPUT  PORTS 
None 

PARAMETERS 
filename 
store_signals 
parent 
fda 

colIect_time 
sfreq 
collects 
LIBRARY  /  FUNCTION 
vulib>mode  3 


Complex 

Complex 

Real 


Filename  to  store  results 
Store  signals? 

Parent  module 
Data  rate 
Collect  time 
Sampling  frequency 
Number  of  collects 
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MODE  4 


DESCRIPTION 

Estimates  the  output  signal  noise  ratio  defined  as 


Var{z  |//„} 


and  stores  the  result  in  a  file  as  a  coefficient  which  permits  extrapolation  of  the  result  to  a 
range  of  collect  times  and  a  range  of  interceptor’s  input  signal  power.  The  extrapolation  to  a 
range  of  signal  powers  is  accurate  only  for  low  signal  powers,  and  the  module  ensures  that 
results  are  presented  over  a  range  where  they  are  accurate. 

The  module  uses  the  results  of  analytic  calculations  to  determine  the  denominator,  thereby 
shortening  the  computation  in  two  ways.  First  the  signal  absent  case  does  not  need  to  be 
simulated,  and  second,  estimates  of  the  mean  converge  faster  than  estimates  of  the  variance, 
so  the  simulation  does  not  have  to  be  as  long  to  achieve  a  given  accuracy. 

The  extrapolation  to  a  range  of  signal  powers  is  accurate  only  for  low  signal  powers,  and 
the  module  ensures  that  results  are  presented  over  a  range  where  they  are  accurate. 

The  user  also  has  the  option  of  storing  intermediate  statistics  on  disk.  The  user  selects  the 
filenames. 

After  the  simulation,  a  window  appears  prompting  the  user  to  obtain  a  Gnuplot  screen  plot 
by  entering 


Plot4  filename 

in  an  operating  system  window,  where  filename  is  the  data  file  name  selected  by  the 
user.  The  plot  ^so  displays  the  name  of  the  underlying  data  file.  To  obtain  a  hard  copy 
in  addition  to  the  screen  plot,  the  user  enters 

Plot4  -h  filename 

If  the  default  printer  has  not  been  set  during  the  session,  the  user  may  have  to  specify  the 
printer  to  be  used  by  entering,  in  the  operating  system  window,  setenv  PRINTER 
followed  by  the  printer  name.  Quitting  the  graphics  window,  or  typing  a  keystroke  in  the 
operating  system  window  causes  the  graphical  output  to  disappear,  and  die  Gnuplot  session 
to  close.  Plots  is  a  Bourne  shell  script  which  selects  and  edits  the  file 
Plots  .hard. template  if  the  flag  -h  is  present,  and 
Plots  .  soft  .template  otherwise,  and  calls  Gnuplot  to  execute  the 
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MODE  4,  continued 


resulting  Gnuplot  script.  Plot 2  .gnu.  The  data  file  and  the  script  files  are  in  the 
directory  /spwdata/vmsign,  but  Plot3  can  be  invoked  from  any  directory 
without  specifying  a  path. 

MENU  PATH 

Interceptor  >  Metrics  >  Mode  4 
INPUT  PORTS 


SIG 

Complex 

NOISE_BW 

Real 

HOLD 

Real 

OUTPUT  PORTS 

None 

PARAMETERS 

com_ran_pow 

Communicator’s  range  input  power.  Fixed  at  1.0  in  SSVM 

filename 

Filename  to  store  results 

fda 

Data  rate 

com_dB_in 

Communicator’s  input  SNR,  dB 

int_ran_adv 

Interceptor’s  noise  advantage,  dB  of  power 

int_nse_adv 

Interceptor’s  noise  advantage,  dB  of  power 

colIect_time 

Collect  time 

collects 

Number  of  collects 

sfreq 

Sampling  frequency 

LIBRARY  /  FUNCTION 
vulib>niode  4 
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MODE  SELECTION 


DESCRIPTION 

Logically  processes  the  parameter  values,  and  selects  one  of  four  ou^uts  to  set  to  low,  and 
sets  the  rest  to  high. 

MENU  PATH 

Interceptor  >  Metrics  >  Mode  Selection 
INPUT  PORT 
None 

OUTPUT  PORTS 

HOLD_l  Real 

HOLD_2  Real 

HOLD_3  Real 

HOLD_4  Real 

PARAMETERS 

long_stat_coll  Long  stationary  collect? 
quad_weak  Quadratic  processing,  with  weak  signal? 

ideal_bkgd  Ideal  background? 

LIBRARY  /  FUNCTION 
vulib  >  select  mode 
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MSK  DEMODULATOR 


DESCRIPTION 

Demodulates  an  MDSK  waveform,  producing  a  binary  sequence. 
MENU  PATH 

Communicator  >  Demodulator  >  MSK 
INPUT  PORT 


IN 

Complex 

OUTPUT  PORT 

OUT 

Real 

PARAMETERS 

fda 

Bit  rate 

sfreq 

Sampling  rate 

LIBRARY  /  FUNCTION 
modli  b/msk_demod 
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MSK  MODULATOR 


DESCRIPTION 

A  minimum  shift  keyed  modulator,  implemented  as  shaped  offset  QPSK. 
MENU  PATH 

Communicator  >  Modulator  >  MSK 


INPUT  PORT 

IN 

Real 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

fda 

Bit  rate 

sfreq 

Sampling  rate 

LIBRARY  /  FUNCTION 
modHb/msk_mod 
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MSK  SOURCE 


DESCRIPTION 

Outputs  a  narrowband  complex  envelope  resulting  from  interleaving  a  pseudorandom  bit 
source,  and  using  the  result  as  the  data  stream  for  a  binary  phase  shift  keying  modulation. 
The  user  sets  the  numbers  of  columns  and  rows  to  define  the  interleaving,  and  accepting  the 
default  value  of  1  for  both  produces  a  non-interleaved  data  stream. 

MENU  PATH 

Communicator  >  Modulated  Source  >  MSK 
INPUT  PORT 
None 

OUTPUT  PORT 

OUT  Complex 

DELAY_OUT  Real 

PARAMETERS 

n  c  o  1  s  Interleaver  columns 

n  r  o  ws  Interleaver  rows 

fda  Bit  rate 

sfreq  Sampling  rate 

LIBRARY  /  FUNCTION 

modulib/msk_source 
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NYQUIST  FILTERED  WAVEFORM 


DESCRIPTION 

Implements  realizable  finite  impulse  response  low  pass  filtering  by  convolving  its  input  with 
a  truncated  sini'jfj  lx 

MENU  PATH 

Utilities  >  Nyquist  Filter 

INPUT  PORT 

IN  Real 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 

b  w  Bandwidth  of  Nyquist  filter 

impulse_length  Length  in  seconds  of  impulse  response 
sfreq  Sampling  rate 

LIBRARY  /  FUNCTION 
utlib/nyq_filter 

Codedin  C 
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COMPOSITE  OUTPUT  SNR  METRICS 


DESCRIPTION 

Logically  processes  the  parameter  values  describing  the  interception  scenario,  and  selects 
one  of  four  ou^ut  SNR  modules,  MODE  1,  MODE  2,  MODE  3,  or  MODE  4,  to  process 
the  incoming  signals. 

MENU  PATH 

Interceptor  >  Metrics  >  Composite  Output  SNR 
INPUT  PORTS 


HI 

Complex 

HO 

Complex 

OUTPUT  PORTS 

None 

PARAMETERS 

filename 

Filename  to  store  results 

com_ran_pow 

Communicator’s  range  input  power.  Fixed  at  1.0  in  SSVM 

int_ran_adv 

Interceptor’s  noise  advantage,  dB  of  power 

long_stat_coU 

Long  stationary  collect? 

quad  _weak 

Quadratic  processing  with  weak  signal? 

ideal_bkgd 

Ideal  background 

colIect_time 

Collect  time 

collects 

Number  of  collects 

sfreq 

Sampling  frequency 

fda 

Data  rate 

com_dB_in 

Communicator’s  input  SNR,  dB 
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OUTPUT  SNR  METRICS,  continued 


USAGE  SNR  METRICS 

SSVM>Demonstrations>CorreIation  Radiometer 
SSVM>Demonstrations>TotaI  Power  Radiometer 

LIBRARY  /  FUNCTION 

vulib>snro  metrics 
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OVERFLOW  PROTECT 


DESCRIPTION 

Replaces  the  input  with  1  during  an  initial  transient  of  length  specified  by  the  user.  After  the 
transient  the  ouq)ut  equals  the  input.  Avoids  transient  overflow  errors  for  example  if  the 
input  is  a  denominator  and  equal  to  zero  at  the  beginning  of  the  simulation. 

MENU  PATH 

Utilities  >  Logic  >  Overflow  Protect 
INPUT  PORTS 

IN  Complex 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

init_protd_tm  Length  if  transient 

sf  re  q  Sampling  rate 

LIBRARY  /  FUNCTION 
utIib/ovrflw_protect 
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PHASOR 


DESCRIPTION 

Outputs  a  phasor  with  frequency  and  phase  exactly  as  given  by  the  respective 
inputs/^  and  (pjj 

—  e 

where is  the  simulation  sampling  frequency.  The  frequency  is  not  adjusted  to  be  the 
inverse  of  an  integer  number  of  simulation  sampling  intervals. 

MENU  PATH 

Utilities  >  Sources  >  Phasor 
INPUT  PORT 

FREQ  Real 

PHASE  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/phasor 
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POWER  COMPENSATION 


DESCRIPTION 

Adaptively  multiplies  the  signal  by  a  scalCT  to  maintain  a  desired  power.  The  user  sets  the 
desired  power  as  a  parameter,  and  the  module  estimates  the  incoming  power.  The  user  can 
also  specify  a  start  time  to  permit  transients  to  decay  before  power  compensation  is 
performed. 

MENU  PATH 

Spread  Spectrum  >  Power  Compensation 
INPUT  PORT 

IN  Complex 

DELAY_IN  Real 

OUTPUT  PORT 

OUT  Complex 

DELAY_OUT  Real 

PARAMETERS 

time_on  Time  on 

ave_pow  Desired  average  power 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 

walib/power_compense 
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QUARTERNARY  DS  SPREADER 


DESCRIPTION 

Multiplies  the  complex  input  by  a  quaternary  complex  pseudonoise  sequence  specified  by 
the  user  as  a  multiple  of  the  data  rate.  If  the  input  signal  has  constant  unit  magnitude,  the 
output  power  is  1. 

Delays  the  output  signal  by  the  input  delay  value  for  automatic  synchronization. 

Ouq)uts  the  spreading  sequence  as  a  global  output  connector  called  ‘code’. 

MENU  PATH 

Spread  Spectrum  >  Direct  Sequence  >  Quaternary 
INPUT  PORT 

IN  Complex 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

n  c  h  Number  of  chips  per  symbol 

fda  Data  rate 

sfreq  Simulation  sampling  frequency 

nsamples  Number  of  samples  in  simulation 

LIBRARY  /  FUNCTION 
walib/quatary_ds 
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QPSK  DEMODULATOR 


DESCRIPTION 

Demodulates  an  quaternary  PSK  waveform,  producing  a  binary  sequence. 

MENU  PATH 

Communicator  >  Demodulator  >  QPSK 
INPUT  PORT 


IN 

Complex 

OUTPUT  PORT 

OUT 

Real 

PARAMETERS 

fda 

Bit  rate 

sfreq 

Sampling  rate 

LIBRARY  /  FUNCTION 

modlib/qpsk_demod 
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QPSK  MODULATOR 


DESCRIPTION 

A  quaternary  PSK  modulator. 

MENU  PATH 

Communicator  >  Modulator  >  QPSK 
INPUT  PORT 

IN  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

fda  Bit  rate 

sfreq  Sampling  rate 

LIBRARY  /  FUNCTION 
modIib/qpsk_mod 
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QPSK  SOURCE 


DESCRIPTION 

Outputs  a  narrowband  complex  envelope  resulting  from  interleaving  a  pseudorandom  bit 
source,  and  using  the  result  as  the  data  stream  for  a  binary  phase  shift  keying  modulation. 
The  user  sets  the  numbers  of  columns  and  rows  to  define  the  interleaving,  and  accepting  the 
default  value  of  1  for  both  produces  a  non-interleaved  data  stream. 

MENU  PATH 

Communicator  >  Modulated  Source  >  QPSK 
INPUT  PORT 
None 

OUTPUT  PORT 

OUT  Complex 

DELAY_OUT  Real 

PARAMETERS 

ncols  Interleaver  columns 

n  r  o  ws  Interleaver  rows 

f  da  Bit  rate 

sfreq  Sampling  rate 

LIBRARY  /  FUNCTION 

moduIib/qpsk_source 
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RANDOM  DIGITS 


DESCRIPTION 

Produces  a  pseudorandom  sequence  of  nonbinary  numbers,  using  a  linear  congruential 
generator.  Output  range  may  be  as  large  as  2  ''24-1. 

MENU  PATH 

Utilities  >  Source  >  Random  Digits 
INPUT  PORT 
None 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 

symbol_size  Range  of  possible  outputs.  One  more  than  the  maximum  possible 

value. 

starting_state  Initial  value. 

LIBRARY  /  FUNCTION 
utIib/random_digits 

Coded  in  C 
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RANDOM  INTERFERENCE 


DESCRIPTION 

Randomly  generates  an  interference  configuration  and  adds  the  resulting  interference  to  the 
signal.  The  user  specifies  the  number  of  interferers,  and  the  maxima  and  minima  of  the 
interferer  width  and  carrier  frequency  distributions,  and  the  mean  and  standard  deviation  of 
the  interferer  power  distribution.  A  configuration  of  interferers  is  then  randomly  generated 
from  the  specified  distributions  and  the  resulting  interference  background  is  added  to  the 
signal.  The  frequency  and  bandwidth  of  each  interferer  are  generated  from  uniform 
distributions  with  upper  and  lower  limits  specified  by  the  user.  The  power  of  the  interferers 
is  generated  from  a  power  law  probability  density,  with  the  mean  power  and  the  variance 
specified  by  the  user.  The  bandwidth  is  implemented  by  BPSK  modulating  each 
interference  carrier  at  a  modulation  frequency  equal  to  the  bandwidth.  TTie  inteiforence 
configuration  is  stored  in  a  file.  The  user  may  specify  the  filename,  which  defaults  to 
rand.conf.. 


The  interference  is  given  by 


N 


=1 


J2n(frfo)n/ 


I  Floori 


V  V 


nW^ 


\ 

I 


where  n  is  the  sample  number,  N  is  the  number  of  interferers,  /  is  the  simulation 
sampling  frequency,  /q  is  the  center  frequency  of  the  communication  signal,  and  ,/,• , 
Wi  and  are  respectively  the  power,  carrier  frequency,  one-sided  bandwidth  and  a 
pseudonoise  binary  data  sequence  of  the  i  th  interferer.  The  data  sequence  has  values  ±  1 . 


MENU  PATH 


Background  >  Interference  >  Random 
INPUT  PORT 
None 

OUTPUT  PORT 

OUT  Complex 

OUTPUT  FILE 

data/rand.conf 

PARAMETERS 

if_num  Interferers 

if_cent_f  req  Configuration  center  frequency 
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RANDOM  INTERFERENCE,  continued 


if_bandwidth 

if_res 

if_niin 

if_niax 

if_mean 

if_std 

simtime 

sfreq 

fees 

LIBRARY  /  FUNCTION 
bklib/ifrand 

Coded  in  C 


Configuration  one-sided  bandwidth 

Frequency  resolution.  Carrier  frequencies  are  selected  at  random 
from  a  gridwidth  spacing  equal  to  the  frequency  resolution 

Minimum  interferer  width 

Maximum  interferer  width 

Mean  interferer  power,  in  dB  signal.  If  the  mean  interferer  power  is  a 
factor  of  is  E  above  the  signal  power,  if_mean  is  equal  to 
lOlog  10(E). 

The  standard  deviation,  in  dB,  of  the  distribution  of  interferer  power. 
If  the  standard  deviation  of  the  power  distribution  is  O,  if_std  is 
equal  to  lOlog  jq  (<^)- 

Simulation  time 
Sampling  rate 
Signal  center  frequency 


C-136 


READ  FROM  FILE 


DESCRIPTION 

Successfully  outputs  the  numbers  form  a  column  of  real  numbers  in  file  named  by  the  user. 
If  the  hold  port  is  high  the  progression  through  the  file  stops  until  the  hold  goes  low  again. 

MENU  PATH 

Utilities  >  Input/Output  >  Read  From  File 
INPUT  PORTS 

HOLD  Real 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 

filename  Name  of  file  containing  column  of  numbers 

LIBRARY  /  FUNCTION 
utlib/read_from_file 

Coded  in  C 
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REAL  INPUT  SWITCH 


DESCRIPTION 

Selects  the  output  from  one  of  the  real  inputs.  Parameters  specify  both  the  number  of 
permitted  inputs,  and  which  line  is  used  for  the  output.  TTie  module  tests  the  value  of  the 
both  parameters,  and  terminates  with  an  error  message  if  either  is  out  of  range.  Setting  the 
parent  module  parameter  permits  the  error  message  to  give  the  containing  module  of  the 
instance  that  has  an  error. 

Coded  in  C. 

MENU  PATH 

Utilities  >  Logic  >  Switch,  Real,  Parameter 
INPUT  PORT 


IN_0 

Real 

IN_1 

Real 

IN_2 

Real 

IN_3 

Real 

IN_4 

Real 

IN_5 

Real 

OUTPUT  PORT 

OUT 

Real 

PARAMETERS 

no_of_inputs  Number  of  inputs  (1-6) 

index  Selection  index  (0-  no_of_inputs-l) 

LIBRARY  /  FUNCTION 
utlib/switch  in  real 
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RECEIVER 


DESCRIPTION 

Simulates  a  real  receiver  or  tuner,  which  filters  the  incoming  real  signal  to  reduce  noise 
contamination  of  the  detection  processing,  and  converts  the  result  to  a  complex  envelope 
form  centered  on  the  center  frequency  of  the  input  passband.  Since  the  entire  simulation 
uses  complex  envelope  representation,  the  simulation  model  of  the  real  bandpass  filtering  is 
complex  lowpass  filtering,  The  Receiver  module,  by  incoiporating  the  Center  Frequency 
between  the  receiver  center  frequency  and  the  sign^  center  frequency.  The  filter  model 
used  is  a  Butterworth  HR  filter  with  an  adjustable  order  set  at  20. 

MENU  PATH 

Utilities  >  Signal  Processing  >  Receiver 
INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

b  w  i  f  Sampling  frequency 

fees  Signal  center  frequency 

feed  Receiver  center  frequency 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/receiver 
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RICIAN  ROC 


DESCRIPTION 

Samples  the  output  of  a  detector  at  intervals  equal  to  the  detector’s  collect  time;  forms 
estimates  of  the  mean  and  variance  of  the  incoming  signal;  treats  the  magnitude  of  the  mean 
as  the  amplitude  of  a  tone,  and  the  variance  as  the  variance  of  Gaussian  noise;  and  computes 
the  receiver  operating  characteristic  (ROC)  of  the  detection  of  the  tone  in  the  Gaussian  noise 
using  an  envelope  centered  on  the  tone  frequency.  The  resulting  ROC  is  an  approximation 
to  the  tme  ROC  because  the  envelope  used  may  not  be  exactiy  centered  on  the  feature 
frequency,  and  the  detector  output  noise  will  in  general  not  be  exactly  Gaussian.  By  setting 
the  flag  in  the  parameter  table,  the  user  may  also  store  the  progressive  estimates  of  tiie 
variance  and  the  amplitude  of  the  mean  for  subsequent  viewing  in  SDE.  The  user  specifies 
the  name,  filename  ,  of  the  file  to  which  the  results  are  written. 

After  the  simulation,  the  user  enters 

Rician  filename 

in  an  operating  system  window,  where  filename  is  the  data  file  name  selected  by  the 
user  in  the  spectrum  module.  The  result  is  a  Gnuplot  screen  plot  showing  spectral  strength 
as  a  function  of  frequency.  The  plot  also  displays  the  name  of  the  underlying  data  file.  To 
obtain  a  hard  copy  in  addition  to  the  screen  plot,  the  user  enters 

Rician  -h  filename 

If  the  default  printer  has  not  been  set  during  the  session,  the  user  may  have  to  specify  the 
printer  to  be  used  by  entering,  in  the  operating  system  window,  setenv  PRINTER 
followed  by  the  printer  name.  Quitting  the  graphics  window,  or  typing  a  keystroke  in  the 
operating  system  window  causes  the  graphical  output  to  disappear,  and  tfie  Gnuplot  session 
to  close. 

The  module  writes  the  mean  and  variance  to  files  called  rician__roc_mean  .dat 
and  rician__roc  mean.dat  respectively.  When  the  simulation  is  over,  the 
custom  coded  module~STORE  RICIAN  ROC  RESULTS  reads  both  files  and  writes  the 
variance  and  the  amplitude  of  the  mean  to  a  single  file  with  name,  filename,  specified 
by  the  user  in  the  parameter  table  of  RICIAN  ROC.  Rician  is  a  Bourne  shell  script 
which  selects  and  edits  the  file  Rician  .hard. template  if  the  flag  "h  is 
present,  and  Rician  .  soft  .template  otherwise,  and  calls  Gnuplot  to  execute  the 
resulting  Gnuplot  script,Rician  .  gnu .  The  data  file  and  the  script  files  are  in  the 
directory  /spwdata/vmsign,  but  Rician  can  be  invoked  from  any  directory 
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RICIAN  ROC,  continued 


A  second  Bourne  shell  script  also  permits  the  user  to  integrate  communication  performance 
results  with  the  Rician  ROC.  After  the  simulation,  the  user  obtains  a  screen  plot  by 
entering 

Rician  filenamel  filename2 

where  the  additional  filename,  f  ilenaine2,  is  the  name  specified  by  the  user  in  the  BIT 
ERROR  RATE  module  to  store  the  bit  error  rate  estimate.  The  functionality  of  BerRic  is 
similar  to  the  functionality  of  Rician  ,  with  the  addition  that  the  bit  error  rate  estimates 
and  the  name  of  the  file  used  to  store  it  are  both  shown  in  the  title  of  the  plot.  As  with 
Rician,  using  the  flag  -h  produces  a  hard  copy  as  well  as  the  screen  plot. 

MENU  PATH 

Interceptor  >  Metric  >  Rician  ROC 

INPUT  PORT 
HI 

OUTPUT  PORTS 
None 

PARAMETERS 

write_moments 
filename 
coIIect_time 
sfreq 
nsamples 

LIBRARY  /  FUNCTION 
vulib/roc  rician 


Complex 


Store  progressive  moment  estimates  if  equal  to  1,  do  not  if  equal  to  0 

Name  of  file  to  store  the  moment  estimates 

Collect  time  used  in  the  detector  module 

Sampling  frequency 

Number  of  iterations  in  the  simulation 
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ROC 


DESCRIPTION 

Forms  histograms  of  the  magnitudes  of  samples  of  two  input,  one  for  signal  present  and  the 
other  for  signal  absent,  and  uses  the  results  to  estimate  the  receiver  operating  characteristic 
(ROC).  The  sampling  interval  is  equal  to  the  collect  time  used  in  the  detector  module.  The 
user  specifies  the  lowest  and  highest  values  of  the  detection  threshold  to  be  considered,  the 
number  of  bins  into  which  the  intervening  range  is  divided,  and  the  name,  filename,  of 
the  file  to  which  the  results  are  written. . 

After  the  simulation,  the  user  enters 

ROC  filename 

in  an  operating  system  window,  where  filename  is  the  data  file  name  selected  by  the 
user  in  the  spectrum  module.  The  result  is  a  Gnuplot  screen  plot  showing  spectral  strength 
as  a  function  of  frequency.  The  plot  also  displays  the  name  of  the  underlying  data  file.  To 
obtain  a  hard  copy  in  adchtion  to  the  screen  plot,  the  user  enters 

ROC  -h  filename 

If  the  default  printer  has  not  been  set  during  the  session,  the  user  may  have  to  specify  the 
printer  to  be  used  by  entering,  in  the  operating  system  window,  setenv  PRINTER 
followed  by  the  printer  name.  Quitting  the  graphics  window,  or  typing  a  keystroke  in  the 
operating  system  window  causes  the  graphical  output  to  disappear,  and  die  Gnuplot  session 
to  close. 

ROC  is  a  Bourne  shell  script  which  selects  and  edits  the  file  ROC .  hard .  template 
if  the  flag  -h  is  present,  and  ROC .  soft .  template  otherwise,  and  calls  Gnuplot  to 
execute  the  resulting  Gnuplot  script,  ROC .  gnu .  The  data  file  and  the  script  files  are  in 
the  directory /spwdata/vmsign,  but  ROC  can  be  invoked  from  any  directory 
without  specifying  a  path. 

MENU  PATH 

Interceptor  >  Metrics  >  Rician  ROC 
INPUT  PORT 


HI 

Complex 

HO 

Complex 

HOLD 

Real 
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ROC,  continued 


OUTPUT  PORTS 
None 

PARAMETERS 


filename 

Name  of  file  to  store  the  moment  estimates 

no_of_bins 

Number  of  bins  dividing  region  between  minimum  and  maximum 
threshold 

min 

Lowest  considered  threshold 

max 

Highest  considered  threshold 

collect_time 

Collect  time  used  in  the  detector  module 

sfreq 

Sampling  frequency 

nsamples 

Number  of  iterations  in  the  simulation 

LIBRARY  /  FUNCTION 
vulib/roc 
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SCANNER 


DESCRIPTION 

Cyclically  scans  a  set  of  frequency  channels  and  stops  when  the  average  power  detected 
exceeds  a  threshold.  The  scanner  frequency  then  remains  fixed  until  Ae  average  power 
drops  below  the  threshold,  when  the  scanning  continues.  The  channel  center  frequencies 
are  either  uniformly  spaced  in  a  band  specified  by  the  user,  or  listed  in  a  file  supplied  by  the 
user. 


The  file  must  exist  in  /caedata/data.  The  power  average  is  a  moving  average  over  an  interval 
equal  to  the  channel  dwell  time.  The  threshold  is  a  user-specified  multiplier  times  die 
average  power  obtained  from  a  reference  signal  that  is  typically  known  to  contain  no  signal. 

The  user  has  the  options  to  store  intermediate  signals  inside  the  scanning  radiometer  and 
inside  the  reference  radiometer.  The  Vulnerability  Metric  stores  the  signals  in 
/spwdata/vmsigs,  with  filenames  based  on  the  stem  supplied  by  the  user. 

MENU  PATH 


Interceptor  >  Detector  >  Radiometer>  Scanner 
INPUT  PORTS 


HI 

Complex 

RFNCE 

Complex 

HOLD 

Real 

OUTPUT  PORTS 

FREQ 

Real 

PARAMETERS 


filename 

file_or_auto 

auto  scan  bw 


Name  of  file  containing  channel  center  frequencies 

1  to  scan  a  band  at  regular  intervals,  0  to  read  center  frequencies 
from  the  file 

Auto  scan  bandwidth 

scan_cent_freq  Center  frequency  of  the  scan 
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SCANNER,  continued 


channel_bw 

Full  bandwidth  of  each  channel 

dwell_tinie 

Time  spent  at  each  frequency  channel 

threshold 

Number  of  dB  above  the  noise  reference  to  halt  scanning 

store_scan_rad 

Store  scanning  radiometer  signals?  (l-yes,  0-no) 

store_ref_rad 

Store  reference  radiometer  signals?  (1-yes,  0-no) 

stem 

Filename  stem  for  storing  signals  (no) 

sfreq 

Sampling  frequency 

fees 

Signal  center  frequency 

USAGE  EXAMPLE 

SSVM>Denionstrations>Scanner 
LIBRARY  /  FUNCTION 
delib/scanner 
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SINE  WAVE 


DESCRIPTION 

Outputs  a  sine  wave  with  frequency  and  phase  determined  by  the  respective  input  signals. 
Unlike  the  SPW  Function  Generator,  the  Vulnerability  Metric  Sine  Wave  module  does  not 
change  the  frequency  to  make  the  period  an  integer  of  sampling  intervals. 

MENU  PATH 

Utilities  >  Sources  >  Sine  Wave 
INPUT  PORTS 


FREQ 

Real 

PHASE 

Real 

OUTPUT  PORT 

OUT 

Real 

PARAMETERS 

sfreq 

Sampling  frequency 

LIBRARY  /  FUNCTION 

utlib/sine_wave 
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SPECTRA 


DESCRIPTION 

Filters  and  subsamples  the  input,  then  forms  the  magnitude  squared  of  the  Fourier 
Transform  of  the  subsampled  signal,  and  finally  performs  a  sliding  block  average  of  the 
resulting  spectrum.  The  process  repeats  and  die  averaged  spectra  are  stored  in  ASCII 
format  in  a  data  file  with  name  specified  by  the  user  as  a  parameter.  The  user  also  specifies 
the  value  of  a  spectrum  width  parameter  and  the  actual  spectrum  width  is  rounded  up  to  a 
submultiple  of  the  sampling  frequency. 


Actual  Spectmm  Width  = 


Sampling  Frequency 


Roor  (Sampling  Frequency/  Spectram  Width  Parameter) 


The  bandwidth  of  the  input  filter  is  set  equal  to  the  Actual  Spectrum  Width  to  eliminate 
aliasing,  resulting  in  spectrum  edges  contaminated  by  the  filter  band  edges.  As  m  aid  to  the 
user,  updating  an  additional  parameter  shows  the  duration  of  signal  input  required  to  form 
each  averaged  spectrum. 


After  the  simulation,  the  user  enters 

Waterfall  filename 


in  an  operating  system  window,  where  filename  is  the  data  file  name  selected  by  the 
user  in  the  spectrum  module.  The  result  is  a  Gnuplot  screen  plot  showing  spectral  strength 
as  a  function  of  firequency.  The  plot  also  displays  the  name  of  the  underlying  data  file.  To 
obtain  a  hard  copy  in  adchtion  to  the  screen  plot,  the  user  enters 

Waterfall  -h  filename 

If  the  default  printer  has  not  been  set  during  the  session,  the  user  may  have  to  specify  the 
printer  to  be  used  by  entering,  in  the  operating  system  window,  setenv  PRINTER 
followed  by  the  printer  name.  Quitting  the  graphics  window,  or  typing  a  keystroke  in  the 
operating  system  window  causes  the  graphical  output  to  disappear,  and  tiie  Gnuplot  session 
to  close. 

Waterfall  is  a  Bourne  shell  script  which  selects  and  edits  the  file 
Waterfall  .hard,  template  if  the  flag  -h  is  present,  and 
Waterfall .  soft  .template  otherwise,  and  calls  Gnuplot  to  execute  the  resulting 
Gnuplot  script,Waterf  all .  gnu .  The  data  file  and  the  script  files  are  in  the  directory 
/spwdata/vmsign,  but  Waterfall  can  be  invoked  from  any  directory  without 
specifying  a  path. 


MENU  PATH 

Interceptor  >  Detector  >  Spectral  >  Spectra 
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SPECTRA,  continued 


INPUT  PORT 


IN 

Complex 

HOLD 

Real 

OUTPUT  PORTS 

None 

PARAMETERS 

spectrum_width 

Approximate  width  of  final  spectra 

fft_size 

Number  of  points  in  spectrum,  power  if  2 

window_tupe 

Windowing  method  used  in  FFT 

spec_avg_block 

Size  of  sliding  block  used  to  average  spectra 

feed 

Input  filter  center  frequency 

filename 

Name  of  file  to  store  the  moment  estimates 

sfreq 

Sampling  frequency 

fees 

Signal  center  frequency 

LIBRARY  /  FUNCTION 
delib/spectra 
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SPECTRAL  ESTIMATE  AND  MATCH 


DESCRIPTION 

Custom  coded  module  preforais  the  frequency  domain  operations  necessary  to  iii5)lement 
the  adaptive  whitening  filter.  It  accepts  a  frequency  domain  block,  updates  the  running 
spectral  estimate,  and  divides  componentwise  by  that  estimated  to  form  a  frequency  domain 
output. 

MENU  PATH 
None 

INPUT  PORTS 

X_re  Real  vector 

X_im  Real  vector 

HOLD  Real 

OUTPUT  PORTS 

Y_re  Real 

Y_im  Real 

PARAMETERS 

blocksize  Number  of  points  used  in  the  FFTs 

average_type  Periodogram  averaging  type,  either  ‘sliding_block’  or  ‘decaying’ 

average_block  Length  of  sliding  window  for  ‘sliding_block’  averaging  type 

daniell_param  Width  of  frequency  domain  smoothing  window 

nexempt  Number  of  points  excluded  in  the  calculation  of  mean  square 

magnitude 

alpha  Decay  parameter  for  ‘decaying’  averaging  type 

sfreq  Simulation  sampling  frequency 
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SPECTRAL  ESTIMATE  AND  MATCH,  continued 


LIBRARY  /  FUNCTION 
prlib/spectral_est 
CODED  INC 
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SPECTRUM 


DESCRIPTION 

Starts  at  a  time  specified  by  the  user,  filters  and  subsamples  the  input,  then  forms  the 
magnitude  squared  of  the  Fourier  Transform  of  the  subsampled  signal,  and  finally  performs 
a  sliding  block  average  of  the  resulting  spectmm.  The  averaged  spectra  are  stored  in  ASCII 
format  in  a  data  file  with  name  specified  by  the  user  as  a  parameter.  The  user  also  specifies 
the  value  of  a  spectmm  width  parameter  and  the  actual  spectmm  width  is  rounded  up  to  a 
submultiple  of  the  sampling  frequency. 


Actual  Spectmm  Width  = 


Sampling  Frequency 


Floor  (Sampling  Frequency/  Spectmm  Width  Parameter) 


The  bandwidth  of  the  input  filter  is  set  equal  to  the  Actual  Spectmm  Width  to  eliminate 
aliasing,  resulting  in  spectrum  edges  contaminated  by  the  filter  band  edges.  As  an  aid  to  die 
user,  updating  an  additional  parameter  shows  the  duration  of  signal  input  required  to  form 
each  averaged  spectmm. 


After  the  simulation,  the  user  enters 


Spectrum  filename 

in  an  operating  system  window,  where  filename  is  the  data  file  name  selected  by  the 
user  in  the  spectrum  module.  TTie  result  is  a  Gnuplot  screen  plot  showing  spectral  strength 
as  a  function  of  frequency.  The  plot  also  displays  the  name  of  the  underlying  data  file.  To 
obtain  a  hard  copy  in  adition  to  the  screen  plot,  the  user  enters 

Spectrum  -h  filename 

If  the  default  printer  has  not  been  set  during  the  session,  the  user  may  have  to  specify  the 
printer  to  be  used  by  entering,  in  the  operating  system  window,  setenv  PRINTER 
followed  by  the  printer  name.  Quitting  the  graphics  window,  or  typing  a  keystroke  in  the 
operating  system  window  causes  the  graphical  output  to  disappear,  and  die  Gnuplot  session 
to  close. 


Spectrum  is  a  Bourne  shell  script  which  selects  and  edits  the  file 
Spectrum,  hard,  template  if  the  flag  -h  is  present,  and 
Spectrum  .soft  .template  otherwise,  and  calls  Gnuplot  to  execute  the  resulting 
Gnuplot  script, Spectrum .  gnu .  The  data  file  and  the  script  files  are  in  the  directory 
/spwdata/vmsign,  but  Spectrum  can  be  invoked  from  any  directory  without 
specifying  a  path. 

MENU  PATH 

Interceptor  >  Detector  >  Spectral  >  Spectra 
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SPECTRUM,  continued 


INPUT  PORTS 


IN 

Complex 

HOLD 

Real 

OUTPUT  PORTS 

None 

PARAMETERS 

start_time 

Time  to  begin  collecting  data  to  determine  averaged  spectram 

spectruni_width 

Approximate  width  of  final  spectra 

fft_size 

Number  of  points  in  spectrum,  power  if  2 

window_tupe 

Windowing  method  used  in  FFT 

spec_avg_bIock 

Size  of  sliding  block  used  to  average  spectra 

feed 

Input  filter  center  frequency 

filename 

Name  of  file  to  store  the  moment  estimates 

sfreq 

Sampling  frequency 

fees 

Signal  center  frequency 

LIBRARY  /  FUNCTION 
delib/spectrum 
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SQUARE  WAVE 


DESCRIPTION 

Outputs  a  square  wave  with  frequency  exactly  as  given  by  the  frequency  input/„ 

Zf^  =  SquareWave (27t/„  ^Ifsa^ 

where  is  the  simulation  sampling  frequency.  In  contrast  to  the  SPW  implementation 
of  function  generators,  the  frequency  is  not  rounded  to  the  inverse  of  an  integer  number  of 
simulation  sampling  intervals. 

MENU  PATH 

Utilities  >  Sources  >  Square  Wave 
INPUT  PORT 

FREQ  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/square_wave 
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STEPPER 


DESCRIPTION 

Produces  a  stepped  output  which  cycles  through  a  set  of  values  specified  by  the  user.  In  the 
pto  mode,  the  values  are  regularly  spaced  and  specified  by  the  starting  value,  the 
increment,  and  the  number  of  values.  Alternatively  the  user  may  store  arbitrary  values  ina 
file  in  the  directory  / caedata/data/,  and  specify  a  scaling  factor  and  a  offset  to  be 
applied  to  the  values  in  the  file.  The  user  controls  the  name  of  the  file,  and  the  dwell  time, 
which  is  the  time  interval  between  changes  in  the  ou^ut. 

A  file  with  the  chosen  name  must  be  present  in  /caedata/data/  ,  even  if  it  is  not 
used  because  the  user  chooses  the  auto  mode. 

MENU  PATH 

Utilities  >  Sources  >  Stepper 
INPUT  PORT 


HOLD 

Real 

OUTPUT  PORTS 

OUT 

Real 

PARAMETERS 

file_or_auto 

1  to  scan  a  band  at  regular  intervals,  0  to  read  center  frequencies 
firom  the  file 

dwell_time 

Time  spent  at  each  fi-equency  channel 

start_freq 

Initial  output  value  for  auto  mode 

step_freq 

Step  in  output  value  for  auto  mode 

no_of_vaIues 

Number  of  values  for  auto  mode 

filename 

Name  of  file  containing  channel  center  frequencies 

scale_factor 

Scaling  applied  to  values  in  file 

offset 

Offset  applied  to  scaled  values  from  file 

sfreq 

Sampling  frequency 

LIBRARY  /  FUNCTION 


utlib/stepper 
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STORE  COMPLEX  VECTORS 


DESCRIPTION 

Provides  optional  storage  of  intermediate  complex  vectors  for  subsequant  analyisis 

MENU  PATH 

Utilities  >  Input/Output  >  Store  Complex  Vectors 
INPUT  PORTS 


INI 

Complex 

IN2 

Complex 

INS 

Complex 

HOLD 

Real 

OUTPUT  PORT 
None 

PARAMETERS 

blocksize  Vector  length 

no  of  vectors  Number  of  vectors  to  store 

filename_stem  Filename  stem 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 

utlib/sto_cmplx_vecs 
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STORE  DETECTOR  DATA 


DESCRIPTION 

Stores  detector  parameter  settings  to  files  on  disk. 

MENU  PATH 

Interceptor  >  Detector  >  Utilities  >  Store  Detector  Data 
INPUT  PORT 


None 

OUTPUT  PORT 


None 

PARAMETERS 


bwif 

fchd 

nsamples 

fda 

colIect_time 


Sampling  frequency 
Interceptor’s  chip  rate  estimate 
Number  of  samples  in  simulation 
Data  rate 

Output  integration  time  in  detector  systems 


LIBRARY  /  FUNCTION 


deIib/store_det_data 
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STORE  ON  TERMINATING 


DESCRIPTION 

Writes  the  signal  value  at  the  IN  port  at  the  end  of  the  simulation  to  a  file.  The  user 
specifies  the  name  of  the  file  as  a  parameter,  and  the  file  is  written  to  the  Graphics 
Directory.  The  Graphics  Data  directory  is  specified  in  the  file  /spwsysA''minclude/all.h,  and 
is  the  same  for  all  graphics  output  data. 

MENU  PATH 

Interceptor  >  Metrics  >  Store  on  Terminating 
INPUT  PORTS 

IN  Real 

HOLD  Real 

OUTPUT  PORT 
None 

PARAMETERS 

filename  Filename 

LIBRARY  /  FUNCTION 
vulib/sto_on_termng 
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STORE  POWER 


DESCRIPTION 

Computes  the  power  of  the  input  signal,  averaged  over  the  whole  simulation  result  in  a  file 
on  disk. 

MENU  PATH 

Utilities  >  Input/Output  >  Store  Power 
INPUT  PORTS 

IN  Complex 

OUTPUT  PORT 
None 

PARAMETERS 

nsamples  Number  of  iterations  in  simulation 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/store_power 
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STORE  REAL  SIGNALS 


DESCRIPTION 

Provides  optional  storage  of  intermediate  real  signals  for  subsequent  analysis 


MENU  PATH 

Utilities  >  Input/Output  >  Store  Real  Signals 

INPUT  PORTS 

INI 

Real 

IN2 

Real 

IN3 

Real 

HOLD 

Real 

OUTPUT  PORT 

None 

PARAMETERS 

no_of_vecs 

Number  of  signals  to  store  (3,  or  less) 

stem 

Filename  stem 

sfreq 

Sampling  frequency 

LIBRARY  /  FUNCTION 

utlib/store  real 
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STORE  RICIAN  ROC 


DESCRIPTION 

After  the  simulation  completes,  writes  the  final  signal  value  at  each  of  the  two  inputs  to  a 
file,  with  the  value  at  the  ABSMEAN  port  in  the  first  line  of  the  file,  and  the  value  at  tiie 
VAR  port  on  the  second  line.  The  user  specifies  the  name  of  the  file  as  a  parameter,  and  the 
file  is  written  to  the  Graphics  Directory.  The  Graphics  Data  directory  is  specified  in  the  file 
/spwsys/  vminclude/all.h,  and  is  the  same  for  all  graphics  ou^ut  data. 

Codedin  C 

MENU  PATH 

Interceptor  >  Metrics  >  Store  Rician  ROC 
INPUT  PORTS 

ABSMEAN  Real 

VAR  Real 

HOLD 

OUTPUT  PORT 
None 

PARAMETERS 

filename  Name  of  file  to  store  amplitude  of  mean  and  variance 

LIBRARY  /  FUNCTION 
vulib/sto  ric  roc 
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STORE  (COMPLEX)  SIGNALS 

DESCRIPTION 

Provides  optional  storage  of  intermediate  (complex)  signals  for  subsequent  analysis 

MENU  PATH 

Utilities  >  Input/Output  >  Store  Signals 
INPUT  PORTS 


INI 

Complex 

IN2 

Complex 

IN3 

Complex 

HOLD 

Real 

OUTPUT  PORT 
None 

PARAMETERS 

no_of_vecs  Number  of  signals  to  store  (3,  or  less) 

stem  Filename  stem 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/store 
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STORE  SPECTRA 


DESCRIPTION 

Stores  the  input  vectors  and  a  vector  of  frequencies  in  columns  in  an  ASCII  file 
/spwdata/vmsigs/.  At  each  iteration  the  elements  of  the  input  vector  at  SPECTRA 
are  reordered  to  place  the  zero  frequency  component  in  the  center  to  the  first  column  in  the 
file.  Parameters  define  the  width  of  the  spectrum  and  frequencies,  and  the  module 
computes  a  vector  of  frequencies  and  appends  the  column  in  the  &e.  The  user  may  zoom  in 
on  the  input  spectrum  by  only  plot  the  spectrum,  centered  at  zero  frequency,  with  the  ratio 
of  the  total  spectrum,  size  being  a  power  of  two  specified  by  the  user.  Cod^  in  C. 

MENU  PATH 

Interceptor  >  Detector  >  Utilities  >  Store  Spectra 
INPUT  PORTS 

SPECTRA  Real  Vector 

TIMES  Real  Vector 

HOLD  Real 

OUTPUT  PORT 
None 

PARAMETERS 

spectrum_width  Width  of  spectra 
freq_zoom  Zoom  factor,  power  of  2 

DEFAULT„VECLEN 

Number  of  points  in  spectrum,  power  of  2 
Input  filter  center  frequency 
Name  of  file  to  store  spectrum 

delib/store_spectra 


feed 

filename 

LIBRARY  /  FUNCTION 
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STORE  THREE  DIMENSIONAL  DATA 


DESCRIPTION 

Treats  each  component  of  the  input  vector  as  the  value  of  a  function  of  two  variables.  The 
value  of  the  first  variable  increments  from  one  iteration  to  the  next,  while  the  value  of  the 
second  variable  increments  from  one  component  of  the  input  vector  to  the  next  The  module 
stores  the  data  in  a  file  as  triplets,  with  one  triplet  per  line,  and  each  triplet  consisting  of  die 
values  of  the  first  variable,  the  second  variable  and  the  vector  component.  Tabs  separate 
numbers  on  a  line.  Coded  in  C. 

MENU  PATH 

Utilities  >  Input/Output  >  Store  3d  Data 
INPUT  PORTS 

IN  Real  vector 

HOLD  Real 

OUTPUT  PORT 
None 

PARAMETERS 

start  1  First  variable  starting  value 

step  1  First  variable  increment  value 

start_2  Second  variable  starting  value 

step_2  Second  variable  increment  value 

DEFAULT_VECLEN 

Number  of  different  values  assumed  by  the  second  variable,  equals 
dimension  of  input  vector 

filename  Name  of  file  to  store  spectrum 

LIBRARY  /  FUNCTION 
utlib/store_3d_data 
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SWITCHED  RADIOMETER  SYSTEM 


DESCRIPTION 

Uses  a  communication  signal  input  and  a  background  input  to  form,  a  signal  present  input 
to  one  detector  and  a  signal  absent  input  to  another  detector.  The  detectors  are  each 
radiometer  detectors,  and  an  additional  input  port  permits  the  user  to  supply  the  switching 
frequency  which  both  adds  to  the  incoming  signd  and  serves  as  the  switched  reference 
source  processing  the  signal  present  case  is  the  Hj  output  and  the  result  of  processing  the 

signal  absent  case  is  the  Hq  output.  Stores  detector  parameter  values  in  files. 

MENU  PATH 


Interceptor  >  Detector  >  Radiometer  >  Switched,  system 
INPUT  PORTS 


SIGNAL 

Complex 

BACKGROUND 

Complex 

RECVR_NOISE 

Complex 

OUTPUT  PORTS 

HI 

Complex 

HO 

Complex 

PARAMETERS 

feed 

Receiver  center  frequency 

bwif 

Receiver 

sw_freq 

Switching  frequency 

wswgn 

1  for  weak  signal,  white  Gaussian  noise  mode,  0  otherwise 

nwssb 

1  for  nonweak  signal,  simulated  background  mode,  0  otherwise 

fees 

Signal  center  frequency 
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SWITCHED  RADIOMETER  SYSTEM,  continued 

sfreq  Sampling  frequency 

nsamples  Number  of  samples  in  simulation 

fda  data  rate 

collect  time  Output  integration  time  in  detector  systems 

LIBRARY  /  FUNCTION 
delib/sw_radiom 
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SWITCHED  RADIOMETER  CORE 


DESCRIPTION 

Stimulates  the  action  of  a  switched  radiometer  detector, with  complex  envelope  additional 
input  permits  the  user  to  supply  receiver  noise,  which  lx)th  adds  signal  and  serves  as  the 
switched  reference. 

MENU  PATH 

Interceptor  >  Detector  >  Radiometer  >  Switched,  core 
INPUT  PORTS 


IN 

Complex 

RECVR_NOISE 

Complex 

OUTPUT  PORT 

OUT 

Complex 

PARAMETERS 

write_signals 

1  to  store  intermediate  signals,  0  otherwise 

feed 

Receiver  center  frequency 

sw_freq 

Switching  frequency 

bwif 

Receiver 

collect_time 

Output  integration  time  in  detector  systems 

sfreq 

Sampling  frequency 

fees 

Signal  center  frequency 

nsamples 

Number  of  samples  in  simulation 

LIBRARY  /  FUNCTION 

delib/sw  radiom  core 
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TAPPED  DELAY  LINE 


DESCRIPTION 

Outputs  the  weighted  sum  of  the  input  signal  and  N  -1  delayed  samples  of  the  input  signal. 
The  number  of  taps,  N,  equals  the  integer  truncation  of  the  product  of  the  spread  ^  sp  and 

the  simulation  sampling  frequency  The  weights  are  read  from  an  input  file  built  by  the 
user.  The  number  of  taps  must  be  less  than  or  equal  to  512.  The  output  signal  yfnj  is 
given  by 


N-1 

y(n)  =ic}^(n-k) 
k  =0 

where  x(n)  is  the  input  signal,  n  is  the  sample  number,  are  the  weights  input  by  the  User, 
and 


N  =  Floor  ) 

MENU  PATH 

Channel  >  Tapped  Delay 
INPUT  PORT 

IN  Complex 

OUTPUT  PORTS 

OUT  Complex 

INPUT  FILE 

data/weights_512 

Each  tap  weight  is  listed  on  a  separate  line  in  the  format  magnitude  phase 

PARAMETERS 

spread  Largest  delay  time,  in  seconds 

sf re  q  Sampling  frequency 

LIBRARY  /  FUNCTION 
chlib/tapdelay 
Codedin  C 
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TOWER  POWER  RADIOMETER  SYSTEM 


DESCRIPTION 

Uses  a  communication  signal  input  and  a  background  input  to  form,  a  signal  present  input 
to  one  detector  and  a  signal  absent  input  to  another  detector.  The  detectors  are  each  total 
power  radiometer  detectors.  The  result  of  processing  the  signal  absent  case  is  the  Hq 

output  Stores  detector  parameter  values  in  files. 

MENU  PATH 

Interceptor  >  Detector  >  Radiometer  >  Total  Power,  system 
INPUT  PORTS 


SIGNAL 

Complex 

BACKGROUND 

Complex 

OUTPUT  PORTS 

HI 

Complex 

HO 

Complex 

PARAMETERS 

feed 

Receiver  center  frequency 

bwif 

Receiver 

collect_time 

Output  integration  time  indetector  systems 

wswgn 

1  for  weak  signal,  white  Gaussian  noise  mode,  0  otherwise 

nwssb 

1  for  nonweak  signal,  simulated  background  mode,  0  otherwise 

fees 

Signal  center  frequency 

sfreq 

Sampling  frequency 

nsamples 

Number  of  samples  in  simulation 

fda 

data  rate 

LIBRARY  /  FUNCTION 
delib/radiometer 
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TOTAL  POWER  RADIOMETER  CORE 


DESCRIPTION 

Simulates  the  action  of  a  radiometer  detector  with  complex  envelope  input  Provides 
optional  storage  of  intermediate  signals. 

MENU  PATH 

Detector  >  Radiometer  >  Total  Power,  core 
INPUT  PORTS 

IN  Complex 

HOLD  Real 

OUTPUT  PORTS 

OUT  Complex 

PARAMETERS 

write  signals  1  to  store  intermediate  signals,  0  otherwise 
f  ce  d  Receiver  center  frequency 

bwif  Receiver 

collect  time  Output  integration  time  indetector  systems 

sfreq  Sampling  frequency 

fees  Signal  center  frequency 

nsamples  Number  of  samples  in  simulation 

LIBRARY  /  FUNCTION 

delib/radiometer_cor 
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TRIGGERED  FREQUENCY  HOPPING 

DESCRIPTION 

Performs  incoherent  frequency  hopping,  shifting  the  incoming  signal’s  frequency  and  phase 
by  random  offsets  which  change  each  time  the  trigger  input  goes  high,  the  input  signal 
has  constant  unit  magnitude,  the  output  power  is  1. 

MENU  PATH 

Spread  Spectrum  >  Trig  Freq  Hopping 
INPUT  PORTS 

IN  Complex 

TRIGGER  Real 

HOLD  Real 

OUTPUT  PORT 

OUT  Complex 

PARAMETERS 

nfreq  Number  of  available  frequency  offsets 

maxfreqs  Maximum  frequency  deviation 

write_f reqs  1  to  store  frequency  offsets 

sfreq  Sampling  frequency 

f  da  Data  rate 

LIBRARY  /  FUNCTION 
walib/trig_inc_fh 
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TRUNCATED  NYQUIST  FILTER 


DESCRIPTION 
Coded  in  C. 

MENU  PATH 

Utilities  >  Trunc  Nyquist  Filter 
INPUT  PORT 

IN  Real 

OUTPUT  PORT 

OUT  Real 

PARAMETERS 

f  c  Cutoff  frequency 

delta_f  Width  of  transition  region 

a  r  Stopband  attenuation 

sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
utlib/kaiser_lpf 
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TWO  DIMENSIONAL  STEPPER 

DESCRIPTION 

Outputs  two  incrementing  reals,  with  the  second  output  cycling  through  its  entire  range  each 
time  the  first  output  changes  once.  The  first  output  does  not  cycle,  but  keeps  incrementing 
indefinitely. 

MENU  PATH 

Utilities  >  Sources  >  2D  Stepper 
INPUT  PORT 
None 

OUTPUT  PORTS 

OUTl  Real 

OUT2  Real 

PARAMETERS 

start_l  First  value  output  at  OUTl 

step_l  Increments  of  the  OUTl  value 

start_2  First  value  output  at  OUT2 

step_2  Increments  of  the  OUT2  value 

no_of_vaIues_2  Number  of  distinct  values  taken  at  OUT2 
LIBRARY  /  FUNCTION 
utIib/two_d_stepper 
Coded  in  C 
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VECTOR  LOG  (BASE  10) 


DESCRIPTION 

Outputs  a  vector  z  which  is  the  pointwise  log  (base  10)  of  the  input  vectors  x. 

H  =logiQ(Xi} 

MENU  PATH 

Utilities  >  Vector  >  Vector  Log  (Base  10) 

INPUT  PORTS 

IN  Real  Vector 

HOLD  Real 

OUTPUT  PORT 

OUT  Real  vector 

PARAMETERS 

DEFAULT_VECLEN 

Length  of  input  and  output  vectors 

LIBRARY  /  FUNCTION 
utlib/vector_loglO 
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VECTOR  MAGNITUDE  SQUARED 


DESCRIPTION 

Outputs  a  vector  z  which  is  the  pointwise  sum  of  squares  of  the  input  vectors  x  and  y. 

2  ,  2 
Zi  =Xi  +yi 

where  the  subscript  i  denotes  the  vector  index. 

MENU  PATH 

Utilities  >  Vector  >  Vector  Magnitude  Sqd 
INPUT  PORTS 

REAL  Real  vector 

REAL  Real  vector 

HOLD  Real 

OUTPUT  PORT 

OUT  Real  vector 

PARAMETERS 

vector_length  Length  of  input  and  output  vectors 
LIBRARY  /  FUNCTION 
utlib/vec_niag_sqd 
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VECTOR  MEAN 


DESCRIPTION 

Outputs  a  vector  z  which  is  the  sliding  block  average  of  the  input  vector  .  If  the  length  of 
the  sliding  block  is  N  and  z^  *  denotes  the  i  th  component  of  the  input  vector  at  the  /i  th 
simulation  iteration,  then 

N-1 

^n,i  ~  J  ^  n-m,i 
m  =0 

MENU  PATH 

Utilities  >  Vector  >  Vector  Mean 
INPUT  PORTS 

REAL  Real  vector 

IMAG  Real  vector 

HOLD  Real 

OUTPUT  PORT 

OUT  Real  vector 

PARAMETERS 

blocksize  Number  of  input  vectors  used  in  sliding  block  average 

DEFAULT_VECLEN 

Length  of  input  and  ou^ut  vectors 

LIBRARY  /  FUNCTION 
utlib/vec_mean 
Coded  in  C 
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VECTOR  SUBSAMPLER 


DESCRIPTION 

Creates  an  output  vector  by  subsampling  the  input  vector.  The  subsampling  ratio  is  the 
integer  truncation  of  the  ratio  between  the  input  and  output  vector  lengths.  Using  vector 
indices  starting  at  1,  the  output  vector  contains  exactly  those  elements  of  the  input  vector 
with  indices  that  are  an  integer  multiple  of  the  subsampling  ratio. 

MENU  PATH 

Utilities  >  Vector  >  Vector  Subsampler 
INPUT  PORTS 

IN  Real  Vector 

HOLD  Real 

OUTPUT  PORT 

OUT  Real  vector 

PARAMETERS 

IN_IOVEC_LEN  Output  vector  length 
LIBRARY  /  FUNCTION 

utlib/vec_subsampler 
Codedin  C 
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VHF  INTERFERENCE 


n 


DESCRIPTION 

Adds  to  the  signal  a  weak  medium  or  strong  interference  background  generated  from  a 
stored  configuration  chosen  to  be  characteristic  of  VHF  operating  environments.  The  user 
selects  the  center  and  width  of  the  operating  band,  and  the  module  generates  the  interferers 
which  lie  in  the  selected  band. 

The  mathematical  form  of  the  interference  the  same  as  that  detailed  for  Random  Interference. 


MENU  PATH 

Background  >  Interference  >  VHF  -  Weak 
Background  >  Interference  >  VHF  -  Medium 
Background  >  Interference  >  VHF  -  Strong 
INPUT  PORTS 


None 

OUTPUT  PORT 

OUT  Complex 

INPUT  FILES 

data/wvhf.conf 

data/mvhf.conf 

data/svhf.conf 


PARAMETERS 

if_cent_freq 

if_bandwidth 

simtime 

sfreq 

fees 


Configuration  center  frequency 
Configuration  one-sided  bandwidth 
Simulation  time 
Sampling  frequency 
Signal  center  frequency 
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VHF  INTERFERENCE,  continued 


LIBRARY  /  FUNCTION 
bklib/ifwvhf 
bklib/ifmvhf 
bklib/ifsvhf 
Coded  in  C 
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WHITENING  EXCISION  CORE 


DESCRIPTION 

Combines  the  two  input  vectors  to  form  a  complex  vector,  and  divides  each  resulting 
component  by  its  magnitude.  The  user  specifies  the  length  of  the  vectors.  The  sp^tral 
resolution,  or  bin  size,  is  the  sampling  frequency  divided  by  the  FFT  size.  If  the  input 
magnitude  of  any  component  is  sm^er  Aan  the  machine’s  arithmatic  resolution,  thd 
magnitude  is  set  to  1,  and  the  phase  is  arbritrarily  set  to  0. 

MENU  PATH 

None 

INPUT  PORTS 

X-re  Real  Vector 

X-im  Real  Vector 

OUTPUT  PORTS 

Y-re  Real  Vector 

Y-im  Real  Vector 

PARAMETERS 

blocksize  Number  of  points  used  in  the  FFTs 

window_type  Type  of  window  used  in  the  forward  FFT 
sfreq  Sampling  frequency 

LIBRARY  /  FUNCTION 
p  rli  b/exwh  ite_cmp 

Coded  in  C 
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WHITE  NOISE 


DESCRIPTION 

Produces  a  discrete  complex  envelope  model  of  continuous  bandpass  white  noise.  The 
noise  can  be  pure  Gaussian,  pure  Gamma,  or  any  mixture  of  the  two,  as  specified  by  the 
user.  The  noise  power  in  both  the  Gaussian  and  Gamma  Noise  modules  is  parameterized 
by  the  comrnunicator’s  input  signal  to  noise  ratio,  which  is  computed  for  the  case  where 
only  one  noise  module  is  present.  Using  more  than  one  noise  module  at  a  time  requires 
manually  determining  the  input  signal  to  noise  ratio. 

MENU  PATH 

Background  >  Noise  >  White 
INPUT  PORT 
None 

OUTPUT  PORTS 

OUT  Complex 

HOLD  Real 

PARAMETERS 

gauss__propn  Gaussian  proportion  of  variance 

Communicator’s  range  input  power,  fixed  at  1 
Communicator’s  input  signal-to-noise  ratio,  dB 
Sampling  frequency 
Data  rate 

LIBRARY  /  FUNCTION 
bklib/white  noise 


com_ran_pow 

snr_in 

sfreq 

fda 
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MISSION 


ROME  LABORATORY 


Mission.  The  mission  of  Rome  Laboratory  is  to  advance  the  science  and 
technologies  of  command,  control,  communications  and  intelligence  and  to 
transition  them  into  systems  to  meet  customer  needs.  To  achieve  this, 
Rome  Lab: 


a.  Conducts  vigorous  research,  development  and  test  programs  in  all 
applicable  technologies; 

b.  Transitions  technology  to  current  and  future  systems  to  improve 
operational  capability,  readiness,  and  supportability; 

c.  Provides  a  full  range  of  technical  support  to  Air  Force  Materiel 
Command  product  centers  and  other  Air  Force  organizations; 

d.  Promotes  transfer  of  technology  to  the  private  sector; 

e.  Maintains  leading  edge  technological  expertise  in  the  areas  of 
surveillance,  communications,  command  and  control,  intelligence,  reliability 
science,  electro-magnetic  technology,  photonics,  signal  processing,  and 
computational  science. 


The  thrust  areas  of  technical  competence  include:  Surveillance, 
Communications,  Command  and  Control,  Intelligence,  Signal  Processing, 
Computer  Science  and  Technology,  Electromagnetic  Technology, 
Photonics  and  Reliability  Sciences. 


